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A. S. Mildvan and C. M. Grisham 
I. Introduction 


The enzyme-catalyzed transfer of simple and substituted phosphoryl 
groups, such as nucleotidyl groups, are among the most ubiquitous and 
fundamental processes in biochemistry. Thus organisms conserve and 
expend energy by phosphoryl transfer to ADP and from ATP. They 
synthesize, repair and degrade nucleic acids by nucleotidyl transfer to 
polynucleotides or to water. Almost all of the enzymes which catalyze 
these reactions require divalent cations as added cofactors or contain 
them as tightly bound components which purify with the enzyme. 
In animal cells the predominant divalent cation available for enzyme 
activation is Mg?+. For magnetic resonance studies with pure enzymes, 
Mg2+ may generally by replaced by the paramagnetic ions Mn2+ or Co2+ 
which provide comparable activation and which serve as probes for 
investigating structural and kinetic properties of the active site (7, 2). 

This article will consider recent advances in the mechanism of model 
phosphoryl transfer reactions, magnetic resonance studies of phosphoryl 
and nucleotidyl transerring enzymes in solution and Xray diffraction 
studies in the crystalline state which have clarified or at least have 
suggested the role of the essential divalent cation. 


II. Chemical Mechanisms of Phosphoryl Transfer 


Phosphoryl and nucleotidyl transfer reactions are nucleophilic displace- 
ments on phosphorus (3—6), and like analogous displacements on carbon 
(7) or on metal ions (5), have been found to take place by mechanisms 
varying between two extreme or limiting cases: 1. In the dissociative 
or Syl mechanism the initial departure of the leaving nucleophile, 
yields the planar triply coordinate metaphosphate anion as a reactive 
chemical intermediate (3, 8), which then combines with the entering 
ligand on either face of the metaphosphate plane 2. In the associative 


! F 5 Vi E ` i 


AN 
X — PAn X: — AR — P— Y (1) 
q — th 


Ol 


or Sy2 mechanism, the entering nucleophile first binds to phosphorus 
along the axis of the trigonal bipyramid resulting in the formation of a 
pentacoordinate phosphorane either as a transition state with a lifetime 
t = 10-13 sec or as a reactive intermediate, (z > 10-13 sec). The leaving 
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nucleophile then dissociates along the opposite axis of the trigonal 
bypyramid. As described by Equation 2, the stereochemistry of the 
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Sy2 displacement is “in line”, that is, the bond angle between entering 
and leaving groups is 180°. 

A stereochemical possibility unique to the associative mechanism is 
“adjacent” displacement in which the bond angle between entering and 
leaving groups is 90°. This possibility results from the ability of the pen- 
tacoordinate phosphorane intermediate to undergo pseudorotation (9, 70) 
which, with certain restrictions (77) interchanges axial and equatorial 
ligands. Nucleophilic displacements on stable phosphoranes appear to 
proceed via an associative mechanism with hexacoordinated phosphorus 
intermediates (72, 73). The biochemical relevance of such displacements 
has not been established. 

As in carbon chemistry (6) and in coordination chemistry (5) mechan- 
isms intermediate between the limiting dissociative and associative 
cases are found. Model studies (3, 4) indicate that both of the limiting 
mechanisms of phosphoryl transfer are substantially accelerated by 
electron withdrawal from the leaving group, lowering its pKa. Such 
electron withdrawal may be accomplished by selective protonation, 
metal coordination or by appropriate covalent substituents. 

Protonation, metal coordination, or further esterification of the trans- 
ferrable phosphoryl! group itself, profoundly retards the Syl mechanisms 
(3), presumably because of the decreased availability of lone pair electrons 
of oxygen for x bonding to phosphorus as metaphosphate is expelled, 
but accelerates the Sy2 mechanism by charge neutralization (6), and 
possibly by electron withdrawal from phosphorus. In addition, metal 
coordination could accelerate Sy2 displacements by five additional 
mechanisms (5). Chelation of the phosphoryl group could 


a) induce strain in the ground state, 


b) stabilize the phosphorane transition state, 
c) or serve as a template for pseudorotation. 
Metal ions could also 
d) coordinate and promote the nucleophilicity of the entering ligand or 


e) simultaneously coordinate the entering ligand and the phosphoryl 
group undergoing nucleophilic attack transforming a bimolecular to a 
unimolecular reaction (3, 6). 
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Since more modes of accelerating Sy2 displacements have been de- 
fined, and since greater stereochemical control is possible with them, it 
might be anticipated that more enzymes would utilize associative rather 
than dissociative pathways. 


III. Metal-Phosphorus Distances 


Almost all enzyme catalyzed nucleophilic displacements on phosphorus 
require divalent metal ions either as added cofactors or in metalloenzymes. 
From the above mechanistic considerations, direct coordination by the 
metal ion of the phosphoryl group undergoing transfer or of the attacking 
nucleophile or both would strongly suggest an Sy2 mechanism and would 
argue against an Syl mechanism, while coordination of the leaving 
group would facilitate either an Syl or Sy2 mechanism. Hence, a deter- 
mination of the Jocation of the metal ion with respect to the phosphoryl 
group undergoing transfer in the active complex would yield insight into 
the mechanism. As pointed out in detail elsewhere (7, 2), nuclear relaxa- 
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Fig. 1. Manganese to phosphorus distances in inner sphere, distorted inner sphere, 
adjacent pyrophosphate and second sphere complexes from crystallographic and 
model building studies. Original references are given in (74—76). Corresponding 
distances in Ca2+ complexes would be 0.2 A greater and in Mg2+ complexes 0.1 A 
smaller 
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tion induced by paramagnetic metal ions such as Mn?+ can be used to 
determine distances between the paramagnetic metal ion and magnetic 
nuclei (19F, 1H, 31 P, 13C) of substrates and substrate analogs in enzyme 
complexes. When such studies are carried out at several magnetic fields 
to yield an independent evaluation of the correlation time, the precision 
attainable in such distance calculations, (< 0.1 A) exceeds that of protein 
crystallography. This results from the fact that the nuclear relaxation 
rate obtained from the NMR data is inversely related to the sixth power 
of the distance while the electron density map derived from the crystallo- 
graphic data is related to the first power of the distance. Small molecule 
crystallography is a higher resolution technique with a precision better 
than + 0.01 A. From the results of small molecule crystallography and 
molecular model studies summarized in Fig. 1 (74—76) the Mn?* to phos- 
phorus distance is ~ 3 A for an inner sphere monodentate complex of 
tetrahedral phosphate; ~ 4 A for an axial oxygen complex of a penta- 
coordinate phosphorane i.e., a “distorted” inner sphere complex; ~ 5 A 
for Mn2+ chelated by an adjacent pyrophosphate in a triphosphate 
chain; and ~ 6 A for a second sphere complex in which an inner sphere 
water or a ligand of comparable size intervenes between the Mn?* and 
phosphate. Corresponding distances in Ca2+ complexes would be 0.2 A 
greater and in Mg2+ complexes 0.1 A smaller. 


IV. Nucleotidyl Transfer Enzymes 


A. Staphylococcal Nuclease 


As a result of extensive studies of its primary, secondary and tertiary 
structure by chemical degradation (77), and synthesis (78), crystallo- 
graphic (79, 20), and magnetic resonance techniques (27), a great deal of 
information is available on the active site of this Ca2+ activated nucleo- 
tidyl transfer enzyme. Comprehensive reviews of each aspect of this work 
have been published (77, 22, 23), and various mechanisms of this reac- 
tion have been considered in detail (76). A reasonable mechanism, 
consistent with all of the structural and kinetic data is an in-line, associa- 
tive (Sw2) nucleophilic displacement on the phosphodiester phosphorus 
by a hydroxyl ligand coordinated to the enzyme bound Cašr (Fig. 2). 
An Syl mechanism is deemed unlikely because of the strong hydrogen 
bonding of the phosphate oxygens by two arginine-guanidinium cations 
from the enzyme ((22), Fig. 2)) which would inhibit an Sy1 mechanism, 
and by the lack of a catalytic role for the nearby and essential Ca2+ 
in such a mechanism. The enzyme bound Ca?+, which is coordinated by 
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Fig. 2. Suggested mechanism for Staphyloccal nuclease based on crystallography 
(19, 20, 22), chemical degradation (77), synthesis (78), nuclear magnetic resonance 
(27, 23, 24), and kinetic studies (25) 


3 or 4 carboxylate ligands, is at a distance from the phosphorus (4.7 + 
0.2 A) which could permit an intervening water or hydroxyl ligand. 
Alternatively the value of 4.7 A could result from a time or space average 
of 50% inner sphere complex (3.1 A) and 50% second sphere complex 
(6.3 A). 


B. DNA Polymerases 
The reaction catalyzed by this class of enzymes is (26, 27): 


Initiator 
3' 


TTT 3 OH . dATP 24 
AGCAACT + {dcp Mg 
TCGTTGACATG eee. ern (3) 
As we vee ce Pe er, 
Template Elongated chain aon 
AGCAACTG + PP; 


TCGTTGACATG 
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The templates which determine the base sequence of the product are 
either polydeoxyribonucleotides or polyribonucleotides, the initiator is 
either an oligodeoxyribonucleotide or an oligoribonucleotide, and the 
substrates are deoxyribonucleoside triphosphates. A proposed distinction 
between enzymes which utilized only polydeoxyribonucleotide templates 
from those which used only polyribonucleotide templates, and the intro- 
duction of the term “reverse transcriptase” for the latter are probably 
not justified on chemical grounds since DNA polymerases and “reverse 
transcriptases’’ can, with varying efficiencies, copy either template (28— 
30). 


It is now clear that tightly bound, stoichiometric Zn?+ is present in 
homogeneous DNA polymereases purified from four diverse sources: 
a bacterial enzyme (E. coli. polymerase I (37, 32)), two viral enzymes 
(T4 bacteriophage (32) and avian myeloblastosis virus (Povesz, B., 
Battula, N. and Loeb, L. A., unpublished results), and an animal enzyme 
(sea urchin polymerase (37)). Kinetic analysis of the inhibition of sea 
urchin polymerase and E. coli DNA polymerase I by o-phenanthroline 
reveals that this inhibition is wholly or in part competitive with respect 
to DNA (27, 37) suggesting that the enzyme bound Zn?+ interacts with 
DNA. This suggestion is supported by the observation that DNA 
polymerase enhances the effect of Zn?* on the nuclear quadrupolar 
relaxation rate of Br- ions indicating that the enzyme bound Zn2+ 
can coordinate Br- ions (32). This enhanced effect is progressively 
reduced on titration with an oligonucleotide, (dT}e-9, to 1/3 of its 
original value with an end point suggesting that 1 molecule of 
oligonucleotide binds to 1 molecule of polymerase and displaces 2/3 
of the Br- from the enzyme bound Zn®+. The presence of Zn?+ at the 
DNA binding site is further supported by the detection of 2 Zn2+ per 
mole of E. colt RNA polymerase (33) which also must bind a DNA 
template. A mechanistic role for Zn2+ rather than simple DNA binding 
is suggested by the specific inhibition by o-phenanthroline of RNA 
chain initiation with RNA polymerase (33), and by the negative finding 
by B. M. Alberts, C. Springgate, L. A. Loeb, and A. S. Mildvan that the 
gene 32 protein, a non-basic and non-catalytic protein which binds 
tightly to single stranded DNA, does not contain significant amounts of 
Zn2+ (less than 0.1 g atoms/mole). 


A reasonable mechanistic role for Zn2+ at the DNA site consistent 
with the above results is to coordinate and promote the nucleophilicity 
of the 3’-OH group of deoxyribose at the growing point of the initiator 
strand ((32, 34), Fig. 3)). Model reactions of this type (35) and their 
microscopic reverse (36, 37) have been described. This role for Zn2+ 
would be assisted by a nearby general base to deprotonate the 3’ hydroxyl 
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group. In this regard it is of interest that several, but not all DNA 
polymerases appear to be sulfhydryl enzymes (26, 27). 
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Fig. 3. Proposed mechanism of E. coli DNA polymerase I showing the role of the 
metal ions (34). Steps I>II>III represent primer elongation and III>IV>V 
represent primer translocation as discussed in Ref. (34) 


Unlike DNA analogs, deoxynucleoside triphosphate substrates do 
not displace Br- ions from the enzyme-bound Zn®+ of E. coli polymerase 
I (32). In an important series of papers reviewed elsewhere (26) Korn- 
berg and his colleagues established that this enzyme binds template 
in absence of the added divalent cation activator, Mg2+ (38), but re- 
quires Mg?+ for the tight binding of the substrates (39). A quantitative 
study of the binding of the divalent activator, Mn2+, and nucleotides 
to E. colt polymerase I using EPR and water proton relaxation (34) 
revealed that the enzyme binds Mn?+ at three types of sites: one very 
tight site (Kp < 10-6 M), five intermediate sites (Kp = 10-5 M), and ~ 12 
weak sites (Kp =10-3 M). From kinetic studies, the tight and inter- 
mediate sites may function catalytically but the weakest sites are 
probably inhibitory. The tight site and intermediate sites all interact 
with deoxynucleoside triphosphates or monophosphates, suggesting 
that the enzyme pre-aligns substrate molecules prior to polymerization 
(34). Moreover, the relaxation rates of water in solutions of ternary com- 
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plexes of deoxynucleoside triphosphates and monophosphates are the 
same (34) and the affinities for triphosphates and monophosphates are 
similar (34, 38, 39) suggesting that both compounds interact with the 
enzyme bound Mn?+ in the same way, presumably through their æ- 
phosphoryl groups. All of these findings are consistent with an associative 
nucleophilic displacement mechanism for DNA polymerase ((34) Fig. 3)). 

The geometry of this displacement, in line or adjacent, is unknown. 


V. Phosphoryl Transfer Enzymes 


A. Pyruvate Kinase 


From 31P-NMR studies at one frequency of the ternary complex of 
muscle pyruvate kinase, Mn, and phosphoglycolate, a competitive analog 
of the substrate phosphoenolpyruvate, a Mn to P distance of 3.5 +0.3 A 
was obtained suggesting a distorted inner sphere complex ((75), Fig. 4)). 
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Fig. 4. Proposed mechanism of pyruvate kinase in A, the phosphoryl transfer 
reaction and B, the enolization of pyruvate, based on nuclear relaxation (75, 43, 
45—48) and chemical modification studies (49—571) 


An alternative explanation of this distance, suggested by the EPR 
spectrum of Mn2+ (40), is that it represents a time average of two 
complexes: 32% inner sphere (2.9 A) and 68% second sphere (6.1 A). 
This possibility arises because of the inherently higher frequency or 
“shutter speed” of the EPR experiment (2 x 101° sec-1) as compared 
with that of the nuclear relaxation experiment (5 x 106 sec-1). Hence 
exchange between the two complexes with a rate constant k such that 
5 x 106 sec-1< k <2 x 1010 sec-1 would result in an average structure 
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in the NMR experiment but not in the EPR experiment. The addition 
of the activator K+ which may form a metal bridge between the enzyme 
and the carboxyl group of phosphoenolpyruvate and its analogs (75) 
causes a reorientation of the inner sphere complex as indicated by proton 
NMR and an increase in the amount of inner sphere complex as suggested 
by EPR (40). A pyruvate kinase-Mn-OgPF bridge complex has pre- 
viously been found by 19F NMR (47). A reevaluation of the Mn to F 
distance in this complex (74) using better estimates of the correlation 
time at one frequency (3.7 x 10-9 sec) and of the extent of formation of 
the ternary complex in the NMR experiment (4%) from computer 
analyzed titration curves (75) yields a Mn to F distance of 4.3 + 1.5 A 
in agreement with the Ca to F distance of 4.2 A in the monodentate 
crystalline Ca—OgPF complex (42). Hence the phosphoenolpyruvate 
analogs phosphoglycolate and FPO3 (which is also a product of the 
fluorokinase reaction) form inner sphere E-Mn-phosphoryl bridge com- 
plexes consistent with an associative or Sy2 nucleophilic displacement 
on the phosphorus of phosphoenolpyruvate by ADP ((43), Fig. 4A)). 
Although phosphoenolpyruvate itself could not be studied by NMR 
because of slow exchange it competes with these analogs in the NMR 
and in kinetic studies (75, 47, 44) and it yields a very similar EPR spec- 
trum for enzyme-bound Mn as that observed in the quaternary pyruvate 
kinase-Mn-K-phosphoglycolate complex (40). 

The geometry of a kinetically active complex of pyruvate kinase, 
Mn?+, K+, phosphate, and pyruvate which catalyzes the enolization 
of pyruvate has been elucidated by 13C and 3!1P NMR ((45) Fig. 4B)). 
The distances from Mn to the carboxyl and carbonyl carbon atoms of 
pyruvate (7.3 + 0.1 A) establish the substrate to be in the second coordi- 
nation sphere, while the lower limit Mn2* to phosphorus distance (4.5 + 
0.4 A) suggests either a distorted inner sphere complex or the rapid 
averaging of 6% inner sphere complex with 94% second sphere complex. 
Mechanisms for pyruvate kinase, consistent with the geometries are 
given in Fig. 4 (75, 43, 45). 

The distance between the bound monovalent cation and bound 
paramagnetic divalent cation, originally determined by Tl* NMR to be 
5 A and 8 A in the presence and absence of phosphoenolpyruvate (46, 47) 
has been confirmed by Na+ NMR (Gupta, R. K., and Mildvan, A. S., 
unpublished), and by frequency dependent proton NMR using mono- 
methylammonium as the activator (48). The possible widespread 
applicability of the use of monomethylammonium rather than Na+ or 
Tl+ would render such distance determinations much easier because of 
the general availability of proton NMR spectrometers (48). 

With respect to the ADP binding site of pyruvate kinase, a functional 
lysyl ammonium group and a thiol group have been detected by a com- 


10 


Metal Ions in Enzyme Catalyzed Phosphoryl Transfer 


bination of methods: chemical modification (49, 50), combined with 
binding studies using EPR and nuclear relaxation (57) and kinetic 
studies of the detritiation of pyruvate, a partial reaction (57). The modi- 
fication of the thiol group consisted of converting it to a disulfide, a 
reaction mediated by 5,5’-dithio-bis-2-nitrobenzoate. While the mod- 
ifications of both groups abolish ADP or ATP binding and preserve the 
binding and detritiation of pyruvate, the thiol group is closer to the 
Mn2+ and phosphoenolpyruvate binding sites as indicated by large 
decreases in the number and exchange rate of rapidly exchanging water 
ligands on Mn?+ and in the affinity for phosphoenolpyruvate in the 
disulfide enzyme (57). 


B. Phosphoenolpyruvate Carboxylase 


PEP Carboxylase catalyzes a reaction homologous with that of pyruvate 
kinase (52) in which an oxygen of bicarbonate attacks the leaving phos- 


Ooc (0. pt05 coo” 
N M2+ 
a to C=O + PO,0H7 
HO—tc + cl = 3 4 
`o" Ken | 4) 
CH, 
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phoryl group of phosphoenolpyruvate and the COe derived from bicar- 
bonate adds to the Cs position of the pyruvate enolate to form oxaloace- 
tate. Mn2+ and substrate binding studies measuring the proton relaxation 
rate of water suggested the formation of an enzyme-Mn-PEP complex 
with carboxylase from peanut cotyledons (52). More recent kinetic data 
with the enzyme purified from spinach provides additional indirect 
evidence for a metal bridge complex to the phosphoryl group (Mzztorko, 
H., Nowak, T., and Mildvan, A. S., unpublished). Thus, changing the 
divalent cation from Mg2+ to Mn?+ produces a 49-fold decrease in the 
dissociation constant of the competitive inhibitor L-phospholactate. 
From the dissociation constants of the inhibitors L-lactate and methyl 
phosphate the factor of 49 is ascribed to a 9-fold tightening of phosphate 
binding but only a 5-fold tightening of carboxyl binding of L-phospho- 
lactate. 


C. Fructose Diphosphatase 


From 31P nuclear relaxation studies (53), an Sy2 mechanism seems likely 
for the reaction catalyzed by the alkaline form of fructose diphosphatase 
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since an enzyme-Mn-fructose-1-phosphate bridge complex with a Mn 
to P distance of 3.3+0.3A has been detected. 


D. Creatine Kinase 


As implied above, a dissociative or Syl mechanism for phosphoryl 
transfer is more difficult to establish than an Sy2 mechanism. An Syl 
mechanism may be operative in the reaction catalyzed by creatine kinase 
as suggested by the observation that planar trigonal monoanions such 
as nitrate and formate (54) which are analogs of metaphosphate, (Fig. 5) 





ADP Nitrate Creatine 


Fig. 5. Proposed structure of the abortive quaternary nitrate complex of creatine 
kinase (55) which may mimic the transition state of the enzymatic reaction, based 
on kinetic (54), EPR (55) and nuclear relaxation studies (88) 


protect the enzyme from modification of its essential thiol. However, the 
finding that M2+-ADP and creatine raise the affinity of the enzyme for 
nitrate as detected kinetically (54) and by proton relaxation (55, 56), 
and, conversely, that nitrate raises the affinity of the enzyme for creatine 
and ADP as detected by kinetics (54) and proton relaxation (53, 57) 
suggest the formation of a complex in which there is significant 
ADP---NOz bonding and creatine~--NOg3 bonding. Such a complex 
would mimic a transition state structure for an in line Sy2 mechanism 
((55) (Fig. 5)). A crucial question in this connection is the location of the 
bound divalent cation. 

A model of the substrates at the active site, based on the super- 
position of the geometry of several inactive and active “ground state” 
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complexes determined by proton NMR suggests that the Mn coordinates 
the leaving ADP group (58). Moreover, the long Mn?+ creatine-CHa- 
distance (9.8 A) did not differ in the enzyme-creatine-Mn-ADP complex 
and in the equilibrium mixture, suggesting that Mn2+ remains on the 
a and ñ phosphorus atoms in the quaternary ATP complex. X-ray 
diffraction studies of the phosphoglycerate kinase-Mn-ADP system 
reveals an enzyme-metal-substrate bridge complex to both phosphoryl 
groups of ADP (59). However, more recent data on the ground state and 
“transition state” complexes of creatine kinase, substrates, and formate, 
indicate the Mn2+ to creatine distance to be <8.5 A (56). The Mn to 
formate proton distance (<5.1 A) is compatible with but does not 
establish direct coordination of formate (56, 60). EPR and proton re- 
laxation data reveal a progressive decrease in the symmetry of the ligand 
field at Mn2+ as the ternary and quaternary complexes are formed 
(54—57, 60), consistent with the entrance of ligands on Mn2* from the 
protein and possibly from the formate. Nucleotides or creatine analogs 
with a higher Vmax appear to form larger fractions of the more asymmetric 
complex (56, 57, 60). Hence while the Mn2+ appears not to coordinate the 
transferable phosphoryl group in the ground state complexes. it may 
well do so in the transition state and in the ternary-enzyme-ATP-Mn 
complex. If this were the case, an Sy2 mechanism would be indicated, 
as with pyruvate kinase, and the metal would have to migrate during the 
course of the reaction. The failure of Cri! ADP (67-63) and Cri ATP 
(63) to serve as substrates for various kinases may be due in part to the 
inability of this trivalent metal to migrate in these stable complexes. 


E. Phosphoglucomutase 


Phosphoglucomutase catalyzes the reversible transfer of a phosphoryl 
group from the 1 or 6 position of glucose 1,6 diphosphate to a serine 
hydroxyl group on the enzyme yielding a phosphoenzyme intermediate 
(64). Direct observation of this enzyme-bound group by 31P-NMR (65) 
established it to be a phosphoserine from its chemical shift and that of 
the protons to which it is coupled and shows it to be highly immobilized 
or “frozen to the enzyme” from its relaxation rates suggesting direct 
activation of the phosphoryl group (W. J. Ray, Jr., and A. S. Mildvan, 
unpublished). The metal ion required for this reaction is bound at a site 
6 A from the phosphoserine P suggestive of a second sphere complex, 
(Fig. 6) although only 20% of the phosphorus was visualized in this 
relaxation experiment, leaving open the possibility of an extensively 
broadened additional signal due to an inner sphere component. A pre- 
dominantly (96%) second sphere complex with only a 4% inner sphere 
contribution is supported by the detection of a 4.8+0.1 A Mn-to-P 
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distance in the ternary complex of the dephospho-enzyme with Mn and 
methyl phosphonate, an analog of phosphate (66). The methyl phospho- 
nate is very weakly bound to the enzyme at this site with a dissociation 
constant of 0.1 M. This site is absent on the phosphoenzyme. On both 
forms of the enzyme, a second binding site for methyl phosphonate is 
found 10.5 +0.5 A from the metal with a lower dissociation constant of 
1 mM which apparently serves to position the binding of either of the 
phosphates of the intermediate glucose 1,6 diphosphate (66, 67). Hence a 
“one dimensional model” of the active site may be constructed (Fig. 6) 


PGM 


——$10.5+0.5A——> 


° 
«—5-6A—>| 


— a 
CH3PO3 
Mn 10PO3 i i 
Kp 10°M OM IO °M 


Fig. 6. One dimensional model of the active site of phosphoglucomutase from 31P 
and proton relaxation studies (65, 66, 68) and Mn?* and CH3PO0$ binding studies 
(67, 68). Kp represents the dissociation constants of Mn2+ and CHsPOs from the 
indicated sites. The weaker CH3PO3 binding site represents the phosphoryl 
transfer site and the stronger binding site positions the substrate 


with the metal ion as the reference point, the weak phosphate binding 
site (corresponding to the phosphoserine phosphate hence to the phos- 
phoryl transfer site) at ~5 A from the metal and the tighter binding 
phosphate positioning site 10—11 A away. The glucose 1,6 diphosphate 
which interchanges its phosphoryl group 100 times faster than it dis- 
sociates (64) must therefore alternatively shift the binding of its two 
phosphates between the two sites without dissociating from the enzyme. 
The calculated distances, together with water proton relaxation (68) 
and molecular model studies, suggest two possible roles for the metal ion 
(Fig. 7). By coordinating the 3-OH group of the sugar the metal-sugar 
bond could serve as the pivot for the phosphate shift. Alternatively, or 
additionally, by coordinating the serine hydroxyl, the metal could 
increase its nucleophilicity on phosphorus and its effectiveness as a leav- 
ing group from phosphoserine (Fig. 7). The second role is suggested by 
the decrease in water relaxation caused by the phosphoenzyme-Mn 
complex as compared with the dephosphoenzyme-Mn complex (67). 
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Similarly the possibility of an Sy2 mechanism is suggested by the im- 
mobilization of phosphorus and by the possibility of an inner sphere 
contribution to the Mn-P distance. 





Fig. 7. Modes of binding of the intermediate glucose 1,6 diphosphate to phospho- 
glucomutase consistent with the geometry of Fig. 6. The enzyme bound Mn?+ 
promotes the nucleophilicity of the serine hydroxyl group and provides a pivot 
for the interchange of the 1 and 6 phosphoryl groups between the phosphoryl 
transfer and phosphoryl positioning sites 


F. Alkaline Phosphatase 


The presence of a metal on a phosphoryl transferring-enzyme provides 
no assurance that the metal is directly involved in phosphoryl transfer. 
Thus with alkaline phosphatase, no direct interactions of Cl- with enzyme 
bound Zn?+ (69) or water with enzyme-bound Mn2* (70) were detected 
by nuclear relaxation. Similarly no direct interaction of phosphate with 
enzyme bound Co2+ (77) or Mn?+ (77, 72) was detected by š1P nuclear 
relaxation. A Mn2+ to phosphate distance of 7.3 A was calculated from 
NMR data on the inactive Mn®*-enzyme (73) indicative of a second 
sphere complex. These results are in accord with crystallographic data 
on the enzyme which at 7.7 A resolution indicate that substrates cannot 
easily gain direct access to the metal site (74). More recent proton 
relaxation studies with the Cu2+ enzyme, which retains 5% of the ac- 
tivity, indicate the presence of a rapidly exchanging axial hydroxyl 
ligand on Cu2+ suggesting that the active metals may promote the 
nucleophilicity of the water molecule which is to attack the phosphorus 
(75). 


G. (Nat + K+)-ATPase 


The membrane-bound Na+ and K+ activated ATPase partially purified 
from kidney microsomes is believed to function in the active transport 
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of Nat and K+ (76). A phosphoenzyme intermediate has been detected 
(77), in which the phosphoryl group is covalently bound in a mixed 
anhydride linkage to the carboxyl group of aspartic acid (78). Inter- 
estingly, this high energy intermediate appears to form from ATP in the 
presence of Nat (77) or from P; in the presence of K+ (79) suggesting 
that the nature of the monovalent cation influences the mode of 
phosphorylation of the enzyme. In addition to Nat and K+ the ATPase 
requires a divalent cation such as Mg2* or Mn?* for activity (76). 

The interactions of Mn2+ with the membrane-bound (Nat +Kt)- 
ATPase from sheep kidney medulla have been examined by kinetic and 
magnetic resonance techniques (80, 87). EPR and water proton relaxation 
rate studies show that the enzyme binds Mn?+ at one tight binding site 
(Kp =0.88 uM). Kinetic studies yield an activator constant for Mn?* 
of 0.88 uM, identifying the one tight Mn?*+ binding site as the active site 
of the ATPase. 

Competition studies indicate that Mg2+ binds at this site with a 
dissociation constant of 1 mM (80) in agreement with its activator 
constant (82). Inorganic phosphate and methylphosphonate bind to the 
enzyme-Mn2+ complex with similar, extremely high affinities suggesting 
covalent bonding, and decrease the enhancement, «*, of 1/Tı of water 
protons. Analysis of the frequency dependence of 1/T 1 of water indicates 
that this change in enhancement is due to a decrease from 4 to 3 in the 
number of rapidly exchanging water protons in the coordination sphere 
of enzyme-bound Mn2+. The relative effectiveness of Na+ and K+ in 
facilitating ternary complex formation with P; was found to be a func- 
tion of $H (80). The fraction of the total de-enhancement produced by 
P; in the presence of Na+ and Mn?+ closely follows the titration curve 
for the second ionization of P;. Analogous behavior at a higher $H was 
observed for methylphosphonate, which has a pKg for its second ioniza- 
tion which is 0.9 H units higher than that of P;. These results indicate 
that Na+ induces the phosphate monoanion to coordinate to the enzyme- 
bound Mn2+, while K+ causes the phosphate dianion to coordinate to 
the enzyme-bound Mn?2+, Thus protonation of an enzyme-bound phos- 
phoryl group would convert a K+ binding site to a Nat binding site. 
Dissociation constants for K+ and Nat estimated from the NMR titra- 
tions (80) agreed with kinetically determined activator constants of 
these ions (83, 84) consistent with binding to the active site. All of these 
results are consistent with the mechanism of Fig. 8 for the ATPase 
and for monovalent cation transport (80). Because of direct coordination 
to Mn?+ (80) the phosphoryl transfer reactions (85) are depicted as Sx2 
processes, In this mechanism the ionization state of phosphate at a single 
enzyme active site controls the binding and transport of Nat or Kt. 
This model also accounts for the order of magnitude weaker binding of 
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Inside || Outside 





Fig. 8. Proposed mechanism for the (Nat + K*)-ATPase and for cation transport 
(80), based on EPR (80), nuclear relaxation (80, 87) and kinetic studies (80, 86), 
the detection of phosphoryl transfer (85) to an aspartyl residue (78), yielding a 
covalent phosphoenzyme intermediate (77, 79). Steps (7—3) show the binding of 
Mn*+, Nat and the phosphate monoanion to the enzyme. Steps (4—77) represent 
the ATPase activity. Active cation transport occurs via steps (4—6) 


Na+ compared to K+, since Na+ binds to a phosphate monoanion and 
K+ binds to a dianion. Recent T1+-NMR data indicates this monovalent 
cation binds 4.0 A from the enzyme bound Mn?* and that this distance 
increases to 5,4 A on binding of phosphate consistent with the mechanism 
of Fig. 8 (87). However, the possibility of T1+ binding at an inhibitory site 
rather than its activating site has not been excluded. 
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VI. Summary and Conclusions 


1. The detection of metal to phosphorus distances intermediate between 
those of inner sphere and second sphere complexes on Staphylococcal 
nuclease, pyruvate kinase, and phosphoglucomutase may well reflect 
the averaging of two or more complexes exchanging rapidly within the 
sampling time defined by the X-ray and NMR methods. Multiple com- 
plexes could result from the rapid exchange of inner and second sphere 
complexes or from the rapid migration of a metal along a polyphosphate 
chain. Alternatively unusual distances may reflect distortion of the metal- 
phosphoryl complexes on the enzyme. Such distortion may be a part of 
the catalytic process itself (87). 

2. Since metal coordination or immobilization of the transferred 
phosphoryl group by multiple hydrogen bonds would inhibit the forma- 
tion of a metaphosphate intermediate in an Syl mechanism and would 
facilitate nucleophilic attack in an Sy2 mechanism, the latter process 
seems likely for the reactions catalyzed by staphylococcal nuclease, 
DNA polymerase, pyruvate kinase, fructose diphosphatase, phospho- 
glucomutase, (Na +K) ATPase and possibly PEP carboxylase. In 
creatine kinase where an Syl mechanism is possible, the enzyme would 
have to prevent access of nucleophiles other than ADP and creatine to 
the reactive metaphosphate intermediate. 

3. In four of the enzymes considered here, Staphylococcal nuclease, 
DNA polymerase (Zn site), phosphoglucomutase, and possibly alkaline 
phosphatase, a reasonable role for the enzyme-bound divalent cation is 
to promote the deprotonation and thereby to increase the nucleophilicity 
of water or of the hydroxyl group of the substrate which is to attack the 
phosphoryl group undergoing transfer. In the microscopic reverse of this 
mechanism, the role of the metal would be to promote the electrophilicity 
of the oxyanion which is leaving phosphorus. This role may also be oper- 
ative in the creatine kinase reaction in which Mn?+ may coordinate the 
leaving ADP portion of ATP in the phosphoryl transfer. In the reactions 
catalyzed by pyruvate kinase, DNA polymerase (Mg site), and (Nat + K+) 
ATPase, the role of the divalent cation appears to be directly at the trans- 
ferred phosphoryl group itself — either to shield charge, to withdraw elec- 
trons, or to stabilize the pentacovalent intermediate or transition state. 
Analogous roles may well be operative with fructose diphosphatase and 
PEP carboxylase in which no ambiguities are introduced by a polyphos- 
phate chain in the substrate. Metals thus appear to play either of two gen- 
eral roles in phosphoryl and nucleotidyl transfer reactions. They may 
coordinate the atom or group which attacks (or leaves) the transferred 
phosphoryl group, or they may coordinate an oxygen atom of the trans- 
ferred phosphoryl group itself. 
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H. P. C. Hogenkamp and G. N. Sando 
I. Introduction 


The ribonucleotide reductase system, consisting of ribonucleotide 
reductase, thioredoxin and thioredoxin reductase, catalyzes the irrevers- 
ible reduction of the four common ribonucleoside-5’-phosphates to the 
corresponding 2’-deoxyribonucleoside-5’-phosphates. This reaction pro- 
vides the deoxy ribonucleotide precursors necessary for DNA synthesis. 


‘ : , ; Ribonucleotide 
Ribonucleotide + Thioredoxin-(SH)¢ 


Reductase 


Deoxyribonucleotide + Thioredoxin-Sa 


Thioredoxin 


Thioredoxin-S, + NADPH + H+ ———— 


Reductase 


Thioredoxin-(SH)a + NADP+ 


The observation that these deoxyribonucleotide precursors exist in 
very low concentration in cells suggests that ribonucleotide reduction 
is a crucial and rate-controlling step in DNA synthesis and cell division. 
Indeed the level of reductase activity has been correlated with DNA 
replication (7), cell proliferation (2) and tumor growth rate (3, 4). 

Ribonucleotide reductase has been studied extensively in two bac- 
terial systems: Escherichia colt (5) and Lactobacillus leichmannii (6). 
These two systems differ in their cofactor requirements and in the level 
of phosphorylation of the substrate. 

Studies on the reductases from animal and plant origin indicate that 
they resemble the enzyme from E. coli, while the adenosylcobalamin- 
dependent reductase, first described in L. leichmannii has been so far 
detected only in procaryotic organisms, in the Euglenophyta, Euglena 
gracilis and Astasia longa (7—10) and in the fungus Pithomyces chartarum 
(77). 

Ribonucleotide reductase activity in non-dividing cells is usually low; 
optimum enzymatic activity In microorganisms is attained during the 
mid-to-late-logarithmic growth phase (72—74). Striking increases in 
reductase activity occur following infection of bacteria or animal cells 
with viruses (75—79). The best characterized ribonucleotide reductases 
from mammalian tissues are those from Novikoff ascites hepatoma 
(20, 27) and regenerating rat liver (22), but reductases have also been 
reported in many rapidly proliferating tissues such as embryonic tissues 
from the chick (23) and the rat (24), leukemic mouse spleen (79), 
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regenerating liver (25), bone marrow (26—28), developing brain (29), 
leukocytes from blood of patients with acute and chronic myelocytic 
leukemia (27), rat hepatomas (3) and Yaba pox virus tumor (78). Ribo- 
nucleotide reductase has also been partially purified from Saccharomyces 
cerevisiae (30) and has been detected in extracts of germinating wheat 
embryo (37) and of root tips of the broad bean, Victa faba (32), in 
preparations of developing adults of the silkworm, Hyalophora cecropia 


(33) and of the sea urchin, Arbacia punctulato, during early embryogenesis 
(34). 


The reduction of ribonucleotides to deoxyribonucleotides is linked 
to NADPH by thioredoxin and thioredoxin reductase. Thioredoxin is a 
small protein which is oxidized from a dithiol to a disulfide form during 
ribonucleotide reduction. The dithiol form of thioredoxin is regenerated 
by NADPH and a specific flavoprotein, thioredoxin reductase. The 
“thioredoxin system” consisting of thioredoxin and thioredoxin reduc- 
tase was first identified as the reducing system from E. coli by Reichard 
and co-workers (35, 36) and both proteins have since been purified to 
homogeneity (37, 38). 


Thioredoxin systems have also been demonstrated in L. leschmannu 
(39), rat Novikoff hepatoma (40), T4 phage infected E. coli (47, 42), 
yeast (43), rat liver (44), calf liver (45) and germinating wheat embryo 
(32). 

Very recently ribonucleotide reduction has been reviewed by Hender- 
son and Paterson (46) and by Follmann (47) ; earlier reviews are those by 
Reichard and co-workers (48—51) and by Blakley (52). Ribonucleotide 


reduction has also been covered in several issues of Annual Review of 
Biochemistry (33—56). 


II. Non-Heme Iron Containing Ribonucleotide Reductases 


The non-heme iron containing ribonucleotide reductases catalyze the 
reduction of ADP, GDP, CDP and UDP to the corresponding 2’-deoxy- 
ribonucleotides. This type of reductase has been purified to homogeneity 
from E. cols B using conventional enzyme purification techniques (57— 
39), while a reductase has been isolated in essentially pure form from T4 
phage infected E. coli B using affinity chromatography on dATP-Se- 
pharose (60). Ribonucleotide reductase has also been purified to homo- 
geneity from extracts of Novikoff rat tumor (67) and from bone marrow 
(28) and more than 150-fold from regenerating rat liver (22). 
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A. Ribonucleoside Diphosphate Reductase from E. col 


Ribonucleotide reductase of E. coli consists of two non-identical subunits, 
which are purified together during the initial steps of the purification 
procedure; however upon further purification the reductase separates 
into two subunits, proteins Bl and B2 (58). The overall yield of the two 
proteins is low, in particular that of protein B1, which readily decom- 
poses into smaller subunits. This decomposition can be counteracted by 
the addition of dithiols. An improved procedure for the purification of 
protein B1 using affinity chromatography on dATP-Sepharose has been 
described by Thelander (59). 


Separately, protein B1 (mol. weight 160,000) and protein B2 (mol. 
weight 78,000) are inactive and do not catalyze any known partial reac- 
tion. However, in the presence of TTP, magnesium ions and dithio- 
threitol proteins B1 and B2 form a 1:1 complex (mol. weight 245,000) 
which is catalytically active (58, 59, 62, 63). 


Molecular weight determinations by ultracentrifugation in guani- 
dine hydrochloride and by dodecylsulfate (SDS) polyacrylamide gel 
electrophoresis showed that protein B 1 consists of two polypeptide chains 
of similar or identical size. Both chains have isoleucine as the carboxyl 
terminus while the amino terminal residues are different: glutamic acid 
and aspartic acid. In the absence of magnesium ions or upon oxidation 
in air protein B1 dissociates into its subunits; the latter dissociation 
involves the oxidation of cysteine residues and leads to inactivation. 
Both dissociation and inactivation can be reversed by incubation in 
presence of dithiothreitol. The allosteric effectors, ATP, dATP, dGTP 
and TTP, which regulate both the level of activity as well as the sub- 
strate specificity of the reductase bind only to protein B1. Subunit 
B2 does not bind any of the nucleotides. In the absence of the effectors 
the enzyme reduces the four ribonucleoside diphosphates at a slow rate. 
Addition of the positive effector ATP stimulates the reduction of the 
two pyrimidine nucleotides CDP and UDP, while the reduction of GDP 
is stimulated by TTP and that of ADP by dGTP. In contrast dATP acts 
as a negative effector and inhibits the reduction of all four substrates. 
The inhibition by dATP is counteracted by ATP and thus the ratio 
between these two nucleotides determines the activity of the enzyme 
(57). The results of binding studies suggest that protein B1 contains 
four binding sites. Two sites with high affinity for dATP and with the 
capacity to bind the other effectors were designated h-sites; while the 
other two sites with a lower affinity for dATP and with no detectable 
capacity to bind TTP or dGTP were designated l-sites. Brown and 
Reichard (64) have suggested that the l-sites regulate the level of activity 
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of the enzyme while the h-sites are involved in the regulation of the 
substrate specificity. 

Protein B2 also consists of two identical or very similar polypeptide 
chains, each with a molecular weight of about 40,000 and contains two 
atoms of inorganic iron (65). The two polypeptide chains do not dissociate 
except under denaturing conditions. The metal can be removed from the 
protein by prolonged dialysis against 8-hydroxyquinoline or by precipita- 
tion with acid ammonium sulfate. Removal of the metal does not change 
the hydrodynamic properties of the protein and indeed “apoprotein B2” 
retains its ability to form a 1:1 complex with protein B1. However 
removal of iron results in loss of enzyme activity which can only be 
restored by treating apoprotein B2 with ferrous ion or less efficiently with 
ferric ion. Protein B2 has a characteristic absorption spectrum with 
maxima at 280 nm (e = 11.2 x 104), 365 nm (6.0 x 103), 410 nm (4.1 x 
103), 480 nm (0.5 x 103), 600 nm (0.3 x 108) and a steep shoulder at 
325 nm (7.6 x 103). The iron content of the protein is closely connected 
with this spectrum, the peaks at 410, 365, 480 and 600 nm disappear 
on removal of the metal and reappear on reactivation with Fet?. Treat- 
ment of protein B2 with hydroxyurea or hydroxylamine, which both 
inhibit ribonucleotide reduction, does not result in removal of the iron. 
However when protein B2 is treated with hydroxyurea or hydroxy- 
lamine the peak at 410 nm, but not the broad band at 365 nm, diminishes 
concomitant with enzyme activity. Protein B2 does not contain labile 
sulfide. The electron spin resonance (ESR) spectrum of protein B2 
recorded from 14 to 293 °K shows a slightly asymmetric doublet centered 
at g=2.0047 with a splitting of 16.5 gauss (66). This paramagnetic 
species is a stable organic free radical because Mössbauer spectra of 
[°’Fe] protein B2 as well as magnetic susceptibility measurements have 
shown that the iron atoms are diamagnetic. The paramagnetic species 
is, however, dependent on the presence of iron; apoprotein B2 does not 
give an ESR signal. The paramagnetic species is destroyed by treatment 
with hydroxylamine or hydroxyurea. On the other hand hydroxylamine 
treatment of [5?Fe] protein B2 does not affect the Mössbauer spectrum. 
These results clearly demonstrate that the enzymatic activity of protein 
B2 depends on the presence of the organic radical and that the radical 
in turn depends on the presence of iron (67). Aged preparations of protein 
B2, which show decreased enzyme activity, with simultaneous loss in 
radical concentration, but without changes in iron concentration, can be 
reactivated by removal of the metal followed by reconstitution with 
ferrous ions. The Mossbauer spectrum of [5?Fe] protein B2 as well as 
the electronic spectrum of protein B2 after treatment with hydroxyla- 
mine are very similar to those of methydroxohemerythrin and oxy- 
hemerythrin, suggesting that protein B2 contains two non-identical 
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high spin Fet3 ions in an antiferromagnetically coupled binuclear com- 
plex. Since protein B2 consists of two apparently identical polypeptide 
chains Atkin and co-workers (67) proposed that each atom of iron is 
bound to a separate polypeptide chain. Furthermore since the Méssbauer 
spectra show the presence of two equal populations of iron, they also 
suggested that by analogy with the proposed oxyhemerythrin structure 
(68) each binuclear iron complex in protein B2 is intrinsically asymmetric 
and contains one iron atom of each class. 

When protein Bl and B2 are mixed in the presence of magnesium 
ions and dithiothreitol active ribonucleotide reductase with an Soo, w 
of 9.7S is formed. Stimulatory effectors, such as ATP and TTP, do not 
effect complex formation. In contrast, in the presence of the negative 
effector, dATP, at concentrations which inhibit enzyme activity a larger 
complex is formed with an S20, w of 15.55. Both complexes contain equi- 
molar amounts of each subunit. A heavy complex is also formed in the 
presence of mixtures of other nucleoside triphosphates which inhibit 
enzyme activity. On the other hand the formation of this heavy inactive 
complex is prevented by ATP at concentrations which reverse the inhibi- 
tion by dATP (63). More recent experiments (59) have shown that the 
interaction between proteins Bl and B2 in the presence of dATP is 
strongly influenced by the presence of sucrose, and indeed in the absence 
of sucrose subunits Bl and B2 with dATP form a complex with an S20, w 
of 22.15. 


B. Ribonucleoside Diphosphate Reductase Induced in E. coli by Infection 
with Bacteriophage T4 


The T4 ribonucleotide reductase is similar to the enzyme from E. colt 
(60). The phage enzyme has a molecular weight of approximately 
225,000 and polyacrylamide gel electrophoresis in sodium dodecylsulfate 
indicated the presence of equimolar amounts of 2 polypeptide chains 
a and f with molecular weights of 85,000 and 35,000 respectively. These 
observations suggest that the native phage reductase is a tetramer, 
woo, with æa corresponding to protein Bl and fe corresponding to 
protein B2 of the E. coli reductase. In contrast to the E. coh reductase 
the phage enzyme does not dissociate into «, and fa during the purifi- 
cation procedure, and thus the phage enzyme does not require magnesium 
ions for activity. The bacterial enzyme requires Mgt? for the formation 
of the BI—B2 complex. The electronic spectrum, ESR spectrum and iron 
content of the phage enzyme are similar to that of protein B2. Competi- 
tion experiments with nucleotide substrates and effectors indicate 
that the phage enzyme contains only one catalytic site and one allo- 
steric site (69). The allosteric properties of T4 ribonucleotide reductase 
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are very similar to those of the enzyme from E. coli. However, the 
action of dATP is opposite for the two reductases: dATP is a prime 
effector for the reduction of the pyrimidine ribonucleotides in the T4 
system, while in the E. coli reductase reaction dATP inhibits the reduc- 
tion of all ribonucleotides. 


C. Ribonucleoside Diphosphate Reductase from Mammalian Sources 


Ribonucleotide reductase from Novikoff ascites hepatoma resembles the 
enzyme from E. colt in its requirement for ATP, magnesium ions and a 
dithiol for optimal activity. The mammalian enzyme is stimulated by 
ferric ions and inhibited by hydroxyurea (27). During chromatography 
in the absence of magnesium ions the reductase separates into two 
components Pl and P2. P1 is retained by dATP-Sepharose and thus 
resembles the B-1 fraction of the E. coli system (67). The reductase from 
regenerating rat liver has been purified extensively and its requirements 
are also reminiscent of those of the bacterial enzyme (22). The purified 
enzyme requires the addition of ATP, Mgt?, Fet® or Fet? and a dithiol 
reducing agent for optimal activity. Incubation of the rat liver enzyme 
in the presence of 8-hydroxyquinoline, hydroxylamine or 0-phenanthro- 
line resulted in an approximately 85% inhibition of reductase activity. 
With some enzyme preparations only marginal stimulation with Fet3 
or Fe+2 ions was obtained, suggesting that, like the bacterial enzyme, 
different enzyme preparations vary with respect to iron content and 
probably organic radical concentration. Recently ribonucleotide reduc- 
tase from rabbit bone marrow has been purified approximately 3600 fold 
(28); chromatography on dATP-Sepharose separates this reductase into 
two inactive protein fractions. The bound fraction (S1) was eluted with 
dATP and corresponds to the B1 protein of the E. coli reductase. The 
unbound fraction (S52) was further purified by chromatography on DEAE- 
Sephadex. During the purification ribonucleotide reductase activity 
became dependent on the addition of ferrous ions, suggesting that the 
iron cofactor is loosely bound to the bone marrow reductase. 


III. Adenosylcobalamin-Dependent Ribonucleotide Reductases 


Ribonucleotide reductases with an absolute requirement for adenosyl- 
cobalamin (Fig. 1) as a coenzyme have been demonstrated only in micro- 
organisms (7). These reductases may be classified into two groups based 
on the nature of the nucleotide substrate utilized: ribonucleoside tri- 
phosphate reductases, which act upon the ribonucleoside triphosphates 
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R = CHsCONH3 
R’ = CHeCHeCONHe 


Fig. 1, Structure of adenosylcobalamin 


and ribonucleoside diphosphate reductases which prefer the ribonucleoside 
diphosphates as substrates. 


A. Ribonucleoside Triphosphate Reductases 


1. Ribonucleoside Triphosphate Reductase from L. leichmannii 


The reductase from L. leichmanni acts only upon the nucleoside tri- 
phosphates and has an absolute requirement for adenosylcobalamin as 
a coenzyme. This reductase has been purified to homogeneity using con- 
ventional purification techniques (70, 77). The enzyme is also retained 
on an adenosylcobalamin-Sepharose column in the presence of a dithiol 
and dGTP, a positive effector, and can be removed from the affinity 
column when dGTP is omitted from the elution buffer. Unfortunately 
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under these conditions the enzyme is slowly inactivated (72). More 
recently affinity chromatography on dGTP-Sepharose has been used for 
the purification of the enzyme. Surprisingly this affinity column 
separated a homogeneous preparation of the enzyme, obtained by con- 
ventional purification techniques into four active fractions. The nature 
of the differences between these four forms has not yet been established. 
(73). 

The molecular weight of the enzyme purified by conventional techni- 
ques is 76,000. Equilibrium sedimentation under various denaturing 
conditions as well as polyacrylamide gel electrophoresis in the presence 
of SDS showed no evidence of subunit structure. The reductase contains 
a single N-terminal amino acid, L-serine, and a single C-terminal amino 
acid, L-lysine, while cyanogen bromide treatment gave nine different 
peptide fragments expected for a protein containing eight methionine 
residues. The results firmly establish that the reductase from L. leich- 
mannii is a Monomeric enzyme with a molecular weight of 76,000. The 
presence of substrates or the modifier, dGTP, did not affect the hydro- 
dynamic properties of the enzyme under the limited conditions investi- 
gated. 

Like the reductase from E. colt, the enzyme from L. leichmanmit is 
subject to allosteric control by the deoxyribonucleotide products. In 
the absence of effectors the enzyme reduces the nucleotide substrates 
at markedly different rates; however these rates are selectively modified 
by deoxyribonucleotide triphosphates which function as specific activa- 
tors. Stimulation of ATP reduction is caused by dGTP, of CTP by dATP 
and of UTP by dCTP (32, 74). The reduction of GTP, the best substrate, 
is not affected significantly by any of the common deoxyribonucleotides. 
Kinetic studies have shown that activation is due to binding at an allo- 
steric site and that such a site may accommodate not only the deoxy- 
ribonucleotide effectors but also the ribonucleotide substrates. Inter- 
action of a deoxyribonucleotide effector with this allosteric site results 
in a marked decrease of the Km for adenosylcobalamin suggesting that 
the enzyme undergoes a conformational change with a consequent 
increase in the affinity for the coenzyme (75). Indeed the observations 
that deoxyribonucleoside triphosphates are able to promote the hydro- 
gen exchange reaction (76), the formation of a radical pair (77) as well 
as the degradation of the coenzyme (78), indicates that these allosteric 
effectors have a pronounced effect on enzyme-coenzyme interaction. 
Results from a study in which analogs of ribonucleotides were tested 
for their ability to act as substrates, activators or inhibitors of ribonucleo- 
tide reduction suggested that the reductase has two effector sites with 
different requirements for the sugar moiety of the nucleoside triphospha- 
tes (79). 
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In contrast, direct binding studies by equilibrium dialysis and differ- 
ential fluorescence methods indicated the presence of only one binding 
site for all the deoxyribonucleotide effectors. Competition experiments 
demonstrated that ribonucleotides are weakly bound to the effector site, 
suggesting that at high concentration the ribonucleotide substrates can 
also modulate enzyme activity. On the other hand in the same study 
binding of ribonucleotide to a substrate site could not be demonstrated 
in the necessarily incomplete reaction systems. The authors suggest 
that their inability to detect enzyme-substrate interaction at a substrate 
site is a consequence of the fact that the ribonucleotide substrate is the 
last reactant to bind in an ordered reaction sequence (77). It should be 
pointed out that these binding studies, showing the presence of only one 
effector site, are done in incomplete reaction mixtures and thus no 
catalysis takes place. 

In contrast with the reductase from E. coli, the ribonucleotide reduc- 
tase system from L. leichmannii does not have an absolute requirement 
for Mgt? (52, 80). 


2. Ribonucleoside Triphosphate Reductase from Euglena gracilis 


The adenosylcobalamin-dependent reductase from this eucaryotic or- 
ganism has recently been purified to homogeneity (87). The molecular 
weight determined by gel filtration was approximately 440,000, while 
polyacrylamide gel electrophoresis in the presence of SDS showed only a 
single band with a molecular weight of approximately 100,000, suggesting 
that the native enzyme is composed of four subunits. The native enzyme 
readily dissociates into smaller units, this dissociation can be counteracted 
by dithiothreitol and dATP. Earlier the molecular weight of a partially 
purified preparation of the reductase was estimated to be in the range of 
140,000 to 150,000 (73). Like the reductase from L. leichmannii the 
enzyme from E. gracilis is a nucleoside triphosphate reductase with an 
absolute requirement for a dithiol and adenosylcobalamin. However, 
unlike the bacterial enzyme the reductase from E. gracilis reduces 
GTP and ATP at approximately the same rate and CTP and UTP at 
approximately half this rate. Furthermore the rate of reduction is not 
stimulated by a specific deoxyribonucleotide but rather the reduction of 
a particular purine ribonucleotide is inhibited by all deoxyribonucleoti- 
des except its own product. TTP appears to be the key modulator of 
reductase activity in E. gracilis, in that this deoxyribonucleotide inhibits 
the reduction of all ribonucleotides. Hamilton (87) has found that CTP 
reduction is stimulated by dATP and ATP; he points out that the homo- 
geneous reductase preparation may have lost a regulatory subunit 
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during the purification procedure and thus may have different kinetic 
properties from the partially purified reductase. 


3, Ribonucleoside Triphosphate Reductase from Thermus X-1 


A deoxyadenosylcobalamin-dependent ribonucleoside triphosphate re- 
ductase has been partially purified from cell free extracts of the extreme 
thermophile, Thermus X-1 (14). The enzyme preparation catalyzed the 
reduction of GIP and CTP at comparable rates, while UTP and ATP 
were reduced at only one-tenth the rate of GTP reduction. Only the 
dithiols could serve as reducing substrates. The enzyme has a temperature 
optimum of 70°, and the allosteric regulation of the enzyme activity is 
also temperature-dependent. The reduction of ATP 1s specifically 
stimulated by dGTP only at a higher temperature. Maximum stimulation 
of ATP reduction is observed at approximately 75°, while no stimulation 
can be detected at 37°. The molecular weight determined by gel filtration 
was approximately 80,000 but no information about the subunit structure 
is yet available. 


4. Ribonucleoside Triphosphate Reductase from Pithomyces chartarum 


The reductase detected in crude extracts of this fungus showed barely 
detectable activity with purine nucleoside diphosphates and was most 
active with ATP and GTP (77). 


B. Ribonucleoside Diphosphate Reductases 


Adenosylcobalamin-dependent reductases which utilize ribonucleoside 
diphosphates as substrates have been detected in three different micro- 
organisms: Rhizobium meliloti (82), Bacillus megaterium (10) and Coryne- 
bacterium nephridu (83). These enzymes require a dithiol as the reducing 
substrate but do not have an absolute requirement for a divalent metal 
ion. The reductase from C. nephridii has been purified to homogeneity 
and preliminary experiments indicate a molecular weight of approxi- 
mately 150,000 and the presence of two subunits. When R. meliloti is 
grown in a cobalt deficient medium the synthesis of ribonucleotide re- 
ductase apoenzyme is derepressed. Furthermore the cells become ab- 
normally elongated and contain subnormal levels of DNA (9). The 
reductase from R. meliloti has been partially purified, while the reductase 
from B. megaterium has only been detected in crude extracts, thus firm 
conclusions concerning their substrate requirements and allosteric control 
must await further studies with more purified preparations. 
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IV. The Catalytic Mechanism 


Although the reduction of ribonucleotides to the corresponding 2’- 
deoxyribonucleotides is catalyzed by enzyme systems differing in their 
cofactor requirements and/or in the level of phosphorylation of the sub- 
strates, the overall reduction process is very similar in all systems. 
For all the systems NADPH is the ultimate reductant, the hydrogen 
is transferred by thioredoxin reductase to thioredoxin, which in turn is 
oxidized by ribonucleotide reductase with the concomitant production 
of a 2’-deoxyribonucleotide. In these systems the thioredoxin-thiore- 
doxin reductase reducing system can be replaced by dithiols such as 
dihydrolipoate or dithiothreitol. 


A. The Stereochemistry of the Reduction 


Initial studies with the enzyme systems from E. coli (84) and from L. 
leichmannii (85, 86) demonstrated that both reductases catalyze the 
transfer of hydrogen from water, via the hydrogens of the dithiol, to the 
deoxyribonucleotide. Chemical and enzymatic degradation studies 
of the 2’-deoxyribonucleotide product showed that tritium from tritiated 
water was associated exclusively with C-2’ of the deoxyribosyl moiety. 
In further studies (86—89) ribonucleotide reduction was carried out in 
deuterium oxide and the deuterated deoxyribonucleotides analyzed by 
nuclear magnetic resonance (NMR) spectroscopy. The NMR spectrum 
of deoxyadenosine, derived from deoxyadenosine triphosphate formed 
by the action of the reductase from L. leichmannii on ATP in deuterium 
oxide showed the simplification of the splitting patterns of the resonances 
of H-1’, H-2’ and H-3’, consistent with the stereospecific introduction 
of one deuterium atom (87). Based on an extension of the Karplus 
equation the downfield octet in the 2’ multiplet was assigned to H-2’ 8 
(Hs) in deoxyadenosine. Since in the NMR spectrum of the enzymatically 
synthesized deuterated deoxyadenosine the upfield octet of the 2’ multi- 
plet was absent, the authors concluded that the deuterium atom was 
introduced in the H-2’« (Hx) position with retention of configuration. 
An independent study of the CDP + dCDP conversion by the reductase 
of E. coli in deuterium oxide (88) led to the same conclusions. However, 
Durham and co-workers pointed out that the Karplus equation when 
applied to furanoses may not be valid; thus these conclusions about the 
stereochemistry of the reductase reaction are not unequivocal. Abrams 
and co-workers (89) using the L. leichmannii reductase and CTP as sub- 
strate based the assignments of the chemical shifts and coupling constants 
of monodeuteriodeoxycytidine on the spectral properties of model 
compounds. Their results also suggested that the 2’-hydroxyl of CTP 
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had been replaced by deuterium with retention of configuration. The 
stereochemistry of the reductase reactions was unequivocally established 
by Fraser-Reid and Radatus (90) who prepared to two specifically labeled 
2’-monodeuteriodeoxycytidines (97) and compared their NMR spectra 
with the spectrum of monodeuteriodeoxycytidine reported by Reichard 
and co-workers (65). This comparison proved conclusively that mono- 
deuteriodeoxycytidine formed enzymatically was (2’R)-[2’-2H]-deoxy- 
cytidine. Independently David and Eustache (92) showed that the deuter- 
ated product of the reductase reaction was not the 2’S-isomer. David and 
Rouzeau (93) used [2’-2H]CDP as a substrate for ribonucleotide reductase 
from E. coli and showed that the label was retained in the product 
without inversion. The 2’-deoxycytidine diphosphate product was trans- 
formed by a series of reactions into 5’-trityl-O2, 3’-cyclouridine, which in 
turn in the presence of potassium tert-butoxide underwent an elimination 
to give radioactive 2’,3’-dideoxy-2’,3’-didehydrouridine. Since most of 
these E2 reactions involve trans eliminations, the presence of the label 
in the product indicated that the tritium label was czs to the leaving group 
thus establishing that the enzymatic reduction had occurred with reten- 
tion of configuration. 

The results from these experiments clearly exclude reaction mecha- 
nisms involving first the elimination of water with the formation of 2’,3’ 
or 1’,2’ double bond in the ribofuranose ring followed by hydrogenation 
of the double bond. They also exclude a classical bimolecular nucleo- 
philic substitution of the hydroxyl group at C-2’ by hydride ion, which 
would proceed with inversion of configuration at C-2’. Because the co- 
factor requirements for the ribonucleotide reductase system from E. colt 
(ATP and Mgt+) resemble those for enzymatic phosphorylation reactions, 
Verheyden and Moffatt (94) suggested that the displacement of the 2’- 
hydroxyl group may be facilitated through the intermediate formation 
of a 2’-phosphate- or pyrophosphate derivative. However cytidine 2’- 
phosphate 5’-diphosphate was not converted to dCDP by the reductase 
from E. colt. A different mechanism was suggested by the similarity 
between the adenosylcobalamin-dependent ribonucleotide reductase 
and dioldehydrase reactions. In the dioldehydrase reaction the partici- 
pation of both hydroxyl functions in the intramolecular redox reaction 
was demonstrated using !8O-labeled propanediol (95). However, studies 
in which 180-labeled adenosine triphosphates were used as substrates for 
the reductase of L. leichmannii indicated that the 3’-hydroxy] function 
does not migrate during the reaction. With [2’-180)ATP as substrate the 
dATP product did not contain any label, while with [3’-18O]ATP as 
substrate all the label could be recovered in the product (96, 97). Further- 
more when [2’-18O|]ATP was incubated with the reductase and all the 
other reaction components except the reducing agent, no label was lost 
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to the solvent, suggesting that in the absence of the reducing agent no 
reversible dissociation of the carbon-oxygen bond occurs. Indirect 
evidence indicates that in the E. colt reductase system also no oxygen 
migration takes place. Adenosine derived from E. colt grown on a medium 
containing [1-180]glucose or [2-180]fructose as the sole carbon source 
was labeled with isotopic oxygen in carbons 2’, 4’ and 5’, whereas deoxy- 
adenosine isolated from DNA contained corresponding amounts of label 
only in carbons 4’ and 5’ (98, 99). 


B. The Nature of the Hydrogen Transfer 


Both ribonucleotide reductases catalyze the transfer of hydrogen from 
the dithiol form of thioredoxin to the 2’-position of the deoxynbonucleo- 
tide product. In the adenosylcobalamin-dependent reductase system the 
mechanism of this hydrogen transfer has been studied extensively. Con- 
siderable evidence has been accumulated which indicates that in all 
adenosylcobalamin-dependent reactions the coenzyme functions as an 
intermediate hydrogen carrier. Thus when ribonucleotide reductase from 
L. leichmannit is incubated with a ribonucleoside triphosphate and a 
dithiol in tritiated water, tritium label is not only incorporated into the 
2’-position of the deoxyribonucleotide but also into the coenzyme. Chem- 
ical degradation of the coenzyme demonstrated that during the reductase 
reaction both hydrogens of the 5’-methylene group bound to cobalt ex- 
change with solvent hydrogens (76, 700, 707, 702). Since these two 
C-5’hydrogens of adenosylcobalamin are not stereochemically equi- 
valent, this observation suggests that during the reductase reaction these 
two hydrogens become stereochemically equivalent. In contrast to the 
other adenosylcobalamin-dependent reactions tritium from [5’-3Ha]- 
adenosylcobalamin is not transferred to the product but only to the sol- 
vent. This is a consequence of the fact that in the reductase reaction 
hydrogen is transferred from a thiol, via adenosylcobalamin to the nucleo- 
tide. Because solvent hydrogens readily exchange with sulfur-bound 
hydrogens, solvent hydrogen is incorporated into the deoxyribonucleo- 
tide product and equilibration between solvent hydrogens and the C-5’ 
methylene hydrogens of the coenzyme occurs. This exchange reaction 
between labeled adenosylcobalamin and the solvent is faster than the 
reduction and thus no label from the tritiated coenzyme is found in the 
deoxyribonucleotide. Furthermore hydrogen exchange between solvent 
and adenosylcobalamin is not dependent on reduction of the substrate. 
Hydrogen transfer from coenzyme to water and from water to coenzyme 
is also promoted by deoxyribonucleotide triphosphates in the absence of 
substrates (76). In spite of this the exchange reaction has been considered 
an integral part of the reductase reaction because both exchange and 
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reduction have an absolute requirement for a dithiol, both reactions are 
catalyzed by the same homogeneous enzyme preparation and treatment 
of the enzyme with protein reagents, such as N-bromosuccinimide, 
affects the rate of reduction and exchange to the same extent. 

Based on isotope exchange experiments with dioldehydrase. Abeles 
(703) postulated the following general mechanism for the hydrogen 
transfer reaction. 


H 
PH x 8 
c + Hp 
H He ~~ H 
SH, + X ° — I: X H, —— PH, + X `° 
Hy Hp 
Hyp N 
H 
PH > e 
b + H 


SH and PH represent substrate and product, XHH is enzyme bound 
adenosylcobalamin and I is an enzyme bound intermediate. The proposed 
reaction involves a two step sequence in which an intermediate is formed 
containing three equivalent hydrogens, two derived from the 5’-methy- 
lene carbon of the coenzyme and one from propanediol. 

A similar intermediate with three equivalent hydrogens has been 
proposed for the ribonucleotide reductase system (76). However, the 
ribonucleotide reductase system is unique in that it catalyzes the hydro- 
gen transfer between two different substrates: a dithiol and a ribonucleo- 
tide. 


He. _ < He 
RSH, + X 55° RSX Hs 
b 


He 
i 7 H 
RSX Ha +RF'OH — R'He + X r, + RS 


Hp NN l n 
R'H, + X i + RS 


R'Hp + X He + RS 
Ha 


RSH and R’OH represent the dithiol and ribonucleotide substrates, 
XHH is the enzyme bound adenosylcobalamin, RSXHHH is an enzyme 
bound intermediate and RS represents the oxidized dithiol. The isotope 
exchange experiments suggest that the enzyme bound intermediate can 
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be generated in presence of only the reducing substrate. Thus the first 
half reaction describes the hydrogen exchange between solvent hydrogens 
and the 5’-methylene hydrogens of the coenzyme, while the second half 
reaction depicts the reduction of the ribonucleotide by the intermediate. 


C. The Nature of the Intermediate(s) 


One of the subunits (protein B2) of ribonucleotide reductase from E. 
colt contains two iron atoms and a variable amount of a paramagnetic 
species. This paramagnetic species, which has been identified as an or- 
ganic radical, is directly correlated with enzyme activity (67). The iron 
of protein B2 does not appear to participate as an electron carrier because 
the part of the electronic spectrum of protein B2 due to iron is not al- 
fected by a variety of oxidizing and reducing agents. Furthermore the 
electronic spectrum is not affected by the addition of thioredoxin or any 
other component — alone or in combination — of the ribonucleotide re- 
ductase system. Atkin and co-workers (67) have suggested that the two 
iron atoms function in the generation as well as the stabilization of the 
organic free radical. Radical formation may involve an iron-catalyzed 
one electron aerobic oxidation of the radical precursor, or a one electron 
reduction of the precursor by Fet? to yield Fet3 in protein B2. According 
to these authors the organic radical is the functional contribution of 
protein B2 to ribonucleotide reductase and they thus suggest that the 
reduction of ribonucleoside diphosphates by the E. coli reductase in- 
volves the participation of radical intermediates. Evidence for the pres- 
ence of paramagnetic intermediates in the ribonucleotide reductase sys- 
tem of L. leichmannii has been obtained from spectrophotometric stopped 
flow studies (704) and ESR measurements (77) by Blakley and co-workers. 
In the presence of dGTP, an allosteric effector, and excess dihydrolipoate 
adenosylcobalamin reacts rapidly (¢1;2 = 17 ms) with equimolar amounts 
of reductase. The changes in the ultraviolet and visible spectrum cor- 
respond to those predicted for the partial conversion of adenosylcobala- 
min to cob(II)alamin, suggesting that the carbon-cobalt bond of the 
coenzyme is cleaved homolytically. The other product of the homolytic 
cleavage, 5’-deoxyadenosyl radical, cannot be trapped in a pool of 5- 
deoxyadenosine. When the temperature of the reaction mixture is lower- 
ed from 37° to 5°, the spectrum of adenosylcobalamin is regenerated, 
indicating that this homolytic cleavage is a reversible process. The homo- 
lytic cleavage of the carbon-cobalt bond of adenosylcobalamin can also 
be followed by ESR spectroscopy if the reaction mixture is freeze- 
quenched in isopentane at 130 °K. The ESR signal is centered on a mean 
effective g value of 2.119. Computer simulation of the spectrum is neces- 
sary for analysis of the hyperfine structure. The authors point out that 
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the principal features of the spectrum are due to a single paramagnetic 
species. On the other hand estimates of the yield of the paramagnetic 
species suggest the production of two unpaired electrons per molecule of 
intermediate. Evidently one paramagnetic species is responsible for the 
major features of the ESR signal while the other gives rise to a very 
broadened signal underlying much of the spectrum. The ESR signal of a 
reaction mixture containing adenosylcobalamin, selectively enriched 
with 90% carbon-13 in the 5’-position of the adenosyl moiety (705), 
does not show any additional hyperfine splitting due to the !%C-nucleus. 
These results suggest that the major features of the ESR spectrum are 
due to an unpaired electron associated with Cot? and that the broadened 
signal is due to the 5’-deoxyadenosyl radical. Thus although the elec- 
tronic spectrum of the intermediate is very similar to that of cob(II)- 
alamin, the ESR spectrum is not at all like that of cob(IT)alamin, either 
free in solution or enzyme-bound. These paramagnetic intermediates de- 
rived from adenosylcobalamin are formed not only in the presence of 
allosteric effectors such as dGTP but also in the presence of the substrates 
GTP and ATP. For instance in the presence of GTP the paramagnetic 
intermediates first accumulate, reach a maximum at 15 to 20 ms and then 
diminish until a low steady state concentration is reached. 

The spectral studies with the E. coli and the L. leichmannii ribonucleo- 
tide reductase system indicate that in both systems one of the ‘‘active 
intermediates’ is an organic free radical. In the E. coli system this or- 
ganic free radical is generated by an unknown interaction between two 
iron atoms and an unknown radical precursor while in the L. leichmannii 
system the organic free radical is generated by a homolytic cleavage of 
the carbon-cobalt bond of adenosylcobalamin. In a subsequent reaction 
either organic radical reacts with a dithiol, such as reduced thioredoxin, 
to yield a reducing agent capable of the hydrogenolysis of the secondary 
alcohol at carbon 2’ of the substrate. 

Two additional paramagnetic species have been detected in the L. 
leichmannii ribonucleotide reductase system. When adenosylcobalamin 
is incubated with a molar excess of reductase, a dithiol and either a ribo- 
nucleoside triphosphate or a deoxyribonucleoside triphosphate, the coen- 
zyme is slowly degraded (#1;2 =30 min) to cob(I])alamin and 5’-deoxy- 
adenosine (78, 706). The ESR spectrum of the paramagnetic cobamide 
shows unique hyperfine and superhyperfine splitting not observed for 
cob(II)alamin in solution or as a solid. 

In the presence of a complete ribonucleotide reductase system: 
reductase, adenosylcobalamin, ribonucleoside triphosphate, appropriate 
allosteric effector and thiol, another type of ESR signal is produced. 
This signal consists of two major lines with mean g values of 2.032 and 
1.965 and has been referred to as the “doublet” spectrum (707). In reac- 
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tion mixtures containing adenosylcobalamin, a molar excess of reductase, 
dGTP, ATP and glutathione this doublet reaches a maximum after 20 
minutes at 37° and then slowly declines; during the same time the highly 
resolved spectrum of cob(II)alamin steadily increases. After 60 minutes 
incubation most of the cobalamin is degraded to cob(II)alamin and 5’- 
deoxyadenosine. Almost identical “doublet” ESR signals have been 
observed in some of the other adenosylcobalamindependent reactions 
such as glyceroldehydrase (708), dioldehydrase (709), and ethanolamine 
deaminase (770), suggesting that in all these reactions a similar para- 
magnetic species is generated. It is noteworthy that the ESR spectrum of 
protein B2 of the E. cols reductase also shows a doublet structure centered 
at g =2.0047, however in the contrast with the doublet of the cobala- 
mindependent reaction the doublet of protein B2 is quite narrow. Of 
course the ESR spectrum of protein B2 does not show the broad re- 
sonance at low field, due presumably to cob(II)alamin. In the L. leich- 
mannii reductase system the most pronounced doublet signal is obtained 
with a poor reducing substrate, glutathione, or with a ribonucleotide 
analog, ara-ATP, that is not a substrate. These results resemble those 
obtained by Abeles and co-workers with dioldehydrase (777). Using a 
substrate analog, 2-chloroacetaldehyde, labeled with 2H or 13C they 
were able to demonstrate that the doublet signal represents a radical 
derived from the substrate analog. Although chloroacetaldehyde is 
not a substrate for dioldehydrase it is nevertheless able to exchange 
hydrogen with the 5’-methylene hydrogens of adenosylcobalamin, 
suggesting that the reaction of chloracetaldehyde with the dioldehydrase 
coenzyme complex produces an intermediate(s) similar to that (those) 
found in the normal catalytic process. Indeed the reaction with chloro- 
acetaldehyde can be regarded as the first half of the hydrogen transfer 
reaction. In the normal catalytic process the steady state level of the 
intermediates is probably too low to be detected by ESR spectroscopy. 
This paramagnetic intermediate in the ribonucleotide reductase system 
has not yet been identified, however, by analogy with the dioldehydrase 
system, the doublet signal may represent a radical derived from the 
reducing substrate or from the ribonucleotide substrate. Hamilton et al. 
(707) have pointed out that the doublet spectrum does not resemble 
published spectra of sulfur radicals. On the other hand the doublet 
spectrum probably does not represent a 2’-deoxyribonucleotide radical, 
because if such a radical were formed with ara-ATP as the nucleotide, 
hydrogen exchange between the 5’-methylene hydrogens of adenosyl- 
cobalamin and the 2’-hydrogen of ava-ATP would be expected. Since 
in the ribonucleotide reductase reaction these hydrogens of the coenzyme 
exchange with solvent hydrogens, solvent hydrogen would be incor- 
porated in the 2’-position of ara-ATP. However, no incorporation of 
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tritium could be demonstrated when this nucleotide was incubated in 
tritiated water with the reaction components required for the genera- 
tion of the doublet ESR signal (772). 

The results described above suggest the presence of at least two inter- 
mediates in the reductase reaction of L. leichmannii and are consistent 
with the following reaction scheme: 
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The first step in this reaction scheme involves the homolytic cleavage 
of the carbon-cobalt bond of the coenzyme to cob(Il}alamin and a 5- 
deoxyadenosyl radical. This intermediate (I) is responsible for the spectral 
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changes observed in the stopped-flow studies and can be detected by 
ESR spectroscopy if the reaction mixture is freeze-quenched. In the 
second step the 5’-deoxyadenosyl radical abstracts a hydrogen atom from 
a dithiol, provided either by the enzyme or by the reducing substrate, 
to yield 5’-deoxyadenosine, and a sulfur radical. 

Both reactions 1 and 2 are reversible and account for the hydrogen 
exchange reaction. In the third reaction the dithiol radical abstracts 
a hydroxyl radical from carbon-2’ of the nucleotide and generates a 2’- 
deoxyribonucleotide radical, water and a disulfide. In the last reaction 
the 2’-deoxyribonucleotide radical abstracts a hydrogen atom from 
5’-deoxyadenosine to yield the 2’-deoxyribonucleotide and the coenzyme 
radical pair. The doublet spectrum is generated by either a sulfur 
radical (II) or a 2’-deoxyribonucleotide radical (III); the unusual line 
width and splitting of this signal are probably due to the interaction of 
the sulfur- or carbon radical with the unpaired electron of cob(I]}alamin 
(113). Slow degradation of all or some of the three intermediates to 
5’-deoxyadenosine and cob(II)alamin probably occurs in vitro. 


V. The Thioredoxin System 


The thioredoxin system, consisting of thioredoxin and thioredoxin 
reductase, was originally discovered as the hydrogen carrier system, 
which provides, with NADPH, the reducing potential for the reduction 
of ribonucleotides (5, 35). Since then considerable evidence has been 
accumulated to indicate that this or a closely related system also parti- 
cipates in a variety of other enzymatic reductions. For instance thiore- 
doxin can function as an electron carrier between NADPH and several 
disulfides, such as insulin, lipoate and oxidized glutathione. Furthermore 
Porqué et al. (114) have shown that thioredoxin and thioredoxin reduc- 
tase from yeast can function as hydrogen carriers in the reduction of 
methionine sulfoxide and sulfate. 


A. Thioredoxin 


Thioredoxin from E. coli is a heat stable small acidic protein (pI 4.4—4.5) 
with a molecular weight of about 12,000 (37, 115—118). It consists of a 
single polypeptide chain with 108 amino acid residues; the complete 
amino acid sequence has been determined (Fig. 2). Thioredoxin can occur 
in the oxidized or reduced state. Reduced thioredoxin contains two 
cysteine residues in positions 32 and 35 which can be oxidized enzy- 
matically or chemically to the disulfide bridge of cystine. Only two 
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1 5 10 15 
H- Ser Asp Lys Ile Ile His Leu Thr Asp Asp Ser Phe Asp Thr Asp 


Leu 
Val 
30 25 20 
Sa Cys Trp Glu Ala Trp Phe Asp Val Leu Ile Ala Gly Asp Ala Lys 
Gly 
Pro 


Cys Lys Met Ile Ala Pro Ile Leu Asp Glu Ile Ala Asp Glu Tyr 


35 40 45 
50 Gin 
Gly 
Thr Gly Pro Asn Gin Asp Ile Asn Leu Lys Ala Val Thr Leu Lys 
65 60 55 
Ala 
Pro 


Lys Thr lle Gly Arg Gly Ile Pro Thr Leu Leu Leu Phe Lys Asn 


70 75 80 

Gly 

85 Glu 
100 

Lys Leu Gln Gly Lys Ser Leu Ala Gly Val Lys Thr Ala Ala Val 

95 90 
Glu 
Phe 
105 


Leu Asp Ala Asn Leu Ala -OH 


Fig. 2. Structure of E. coli thioredoxin (37) 


amino acids — glycine and proline — are spaced between the half cystines; 
the proline residue holds the disulfide bridge in a fixed and exposed posi- 
tion. The N-terminal half of thioredoxin contains twelve acidic amino 
acids and four basic residues; on the other hand the C-terminal region 
contains many hydrophobic residues. This distribution of amino acid 
residues suggests that the N-terminal part of thioredoxin, including 
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the disulfide bridge of the active center is exposed to the solvent, while 
the C-terminal region is buried m the inside of the molecule by its 
hydrophobic residues. The redox potential of the couple thioredoxin 
Sa/thioredoxin (SHg) is —0.26 V, pH 7.0, 20° indicating that reduced 
thioredoxin is a better reducing agent than glutathione and cysteine. 
The two tryptophan residues of thioredoxin (Trp-28 and Trp-31) are 
located close to the disulfide bridge and seem to play a major role in 
maintaining an enzymatically active conformation of the protein. 
Chemical modification of both tryptophans to I-formyl-tryptophanyl 
residues (779) or to oxindolyl residues (720) resulted in complete inactiva- 
tion. Tryptophan fluorescence studies of the oxidized and reduced form 
of thioredoxin suggested that one or both of the residues are involved 
in a conformational change that occurs when the disulfide bond is reduced 
(127, 122). This conformational change of the protein must be localized 
in nature because the sedimentation behavior, the fluorescence polariza- 
tion, the ORD and the CD do not change on reduction of the disulfide 
bridge. Reichard (5) suggests that this conformational change is involved 
in the binding of the oxidized and reduced forms to their respective 
proteins in the ribonucleotide reductase system. The reduced form of 
thioredoxin should bind preferentially to ribonucleotide reductase, while 
the oxidized form should bind to thioredoxin reductase. Thioredoxin 
has been crystallized from alcoholic solutions in the presence of cupric 
ions. The crystals belong to the monoclinic space group C2 with unit 
cell dimensions a = 89.5 Å, b =50.8 A, c=60.2 Å and 8 =113.5° and have 
two molecules of thioredoxin per asymmetric unit (723). 

Infection of E. colt with bacteriophage T4 am 122 induces the syn- 
thesis of a new ribonucleotide reductase and a new thioredoxin; these 
two proteins are able in the presence of E. col thioredoxin reductase to 
catalyze the reduction of ribonucleotides. Phage-induced ribonucleotide 
reductase has a specific requirement for phage-induced thioredoxin 
and is unable to use reduced E. coli thioredoxin as a hydrogen donor; 
conversely, E. colt ribonucleotide reductase is unable to use phage 
thioredoxin (47, 42). Like thioredoxin from E. coli, phage-induced 
thioredoxin is a small heat stable protein containing a single disulfide 
bridge made up of two half cystine residues. However, the two thiore- 
doxins differ in size, amino acid composition and redox potential. T4 
thioredoxin has a molecular weight of 10,400 and consists of a single 
polypeptide chain with 88—89 amino acid residues. It contains three 
methionine residues but is completely lacking in tryptophan. Indeed 
the amino acid sequence of the active center of T,4-thioredoxin is quite 
different from that of E. coli thioredoxin (724). The absence of trypto- 
phan is of interest because the two tryptophan residues in the bacterial 
thioredoxin are situated close to the disulfide bridge and have been 
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implicated in the reaction mechanism. The redox potential of T4 thio- 
redoxin is —0.23 V at pH 7.0 and 20°, significantly more positive 
than that of E. coli thioredoxin, indicating that oxidized T4 thioredoxin 
is able to oxidize reduced E. colt thioredoxin (724). 

Both E. coli and T4 thioredoxin have recently been purified by immu- 
noabsorbent affinity chromatography. The specific immunoadsorbents 
were prepared by coupling the «-globular fractions of rabbit anti-thio- 
redoxin antisera to Sepharose. This technique makes it possible to isolate 
both thioredoxins from phage-infected cells by two consecutive chroma- 
tography steps on the specific adsorbents (725, 726). 

Thioredoxin from yeast has been obtained in two forms (I and II) of 
which thioredoxin II has been purified to homogeneity (43). Both yeast 
thioredoxins are able to serve as hydrogen donors for the ribonucleotide 
reductase from E. coli. The molecular weight (12,600) of thioredoxin II 
is similar to that of thioredoxin from E. coli. Although both yeast thio- 
redoxins contain only one tryptophan residue and although their amino 
acid compositions differ markedly from that of E. col thioredoxin, the 
amino acid sequences around the disulfide bridge of these three thio- 
redoxins are identical: 


-Trp-Cys-Gly-Pro-Cys-Lys- (722): 


The oxidation reduction potential of the yeast thioredoxin-II (SH) 2/ 
thioredoxin-Sg couple is —0.24 V at pH 7. 

A thioredoxin has also been partially purified from L. leichmannii 
(39). Although this thioredoxin is similar in size (MW approximately 
12,000) to the E. col thioredoxin, it is not able to function as a substrate 
for the E. colt thioredoxin reductase. More recently thioredoxins have 
been purified to homogeneity from rat Novikoff ascites hepatoma (727) 
and from calf liver (45). The properties of these two thioredoxins are 
quite similar (MW 11,400 and 12,000 respectively). The amino acid 
compositions however, are different. For instance the Novikoff tumor 
thioredoxin contains six half-cysteine residues, whereas calf liver 
thioredoxin has only four. Both these thioredoxins have a tendency 
to aggregate in the oxidized form. This aggregation 1s probably due to 
mixed sulfide formation between the additional sulfhydryl groups. 
Like yeast thioredoxin II, calf liver thioredoxin contains only one 
tryptophan residue; unfortunately the tryptophan content of Novikoff 
tumor thioredoxin was not determined. 

Of the thioredoxins studied thus far only the T4 phage thioredoxin 
is unable to serve as the reducing substrate for E. colt ribonucleotide 
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reductase. T4 thioredoxin does not contain tryptophan and the amino 
acid sequence around the disulphide bridge is completely different from 
that of the thioredoxins from E. coli, yeast, Novikoff hepatoma, calf 
liver and rat liver (44). Presumably these thioredoxins are able to interact 
with E. coli ribonucleotide reductase because their active sites are very 
similar or identical, as is the case for the thioredoxins of E. coli and 
yeast. 


B. Thioredoxin Reductase 


Thioredoxin reductase catalyzes the reduction of the oxidized form of 
thioredoxin by NADPH (36). 


Thioredoxin-S3 + NADPH + H+ —— Thioredoxin-(SH)g + NADP+ 


The enzyme is specific for the reduction of the disulfide bond of 
thioredoxin and NADPH is the preferred hydrogen donor. The hydrogen 
transfer reaction is stereospecific for the B-position (Hs) of the nico- 
tinamide ring of NADPH (728). Although the purified enzyme shows 
no activity with glutathione and lipoate, thioredoxin reductase will 
catalyze the reduction of a variety of disulfides by NADPH in the pres- 
ence of catalytic amounts of thioredoxin. Thioredoxin reductase from 
E. coli has been purified to homogeneity in two laboratories (38, 729) 
and has been studied extensively (130—133). The molecular weight of 
the enzyme determined by ultracentrifugal methods was found to be 
66,000 (728). It contains two molecules of FAD per molecule of enzyme 
and its visible spectrum is very similar to the spectra of glutathione 
reductase and lipoyldehydrogenase. Neutron activation analysis showed 
that the enzyme does not contain a stoichiometric amount of a metal. 
The enzyme consists of two identical or very similar polypeptide chains 
held together by non-covalent bonds. Each of the subunits contains one 
molecule of FAD. The molecular weight of a subunit based on amino 
acid analysis and FAD content was 36,500. Each subunit contains 
three sulfhydryl groups and one disulfide bridge (8). Earlier determina- 
tions had shown the presence of only 4 half-cystines per mole of FAD 
(730, 131). The sulfhydryl groups can be reacted with alkylating agents 
only after complete denaturation of the enzyme with guanidine hydro- 
chloride suggesting that they are buried in the enzyme structure. 
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Thelander (131) has proposed the following structure for thioredoxin 
reductase: 


S S 
= E j. 


SH SH SH 
SH SH SH 


S S 


Each subunit of the enzyme functions in a series of oxidation-reduc- 
tion reactions as outlined below (730): 
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The normal catalytic process involves the sequence I--II>III+>IV-I, 
while the steps [V+V—>VI-III may also be part of the overall reaction. 
Thelander (137) has pointed out that the reaction mechanism of thio- 
redoxin reductase differs from those of glutathione reductase and 
lipoyldehydrogenase because during the anaerobic titration of thioredoxin 
reductase no flavin semiquinone intermediate could be detected. Although 
thioredoxin reductase contains two subunits only one octapeptide 
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containing the active center of the enzyme could be isolated, indicating 
that the amino acid sequence around both disulfides is identical. The 
amino acid sequence of the cystine peptide was determined to be (732, 
133): 

S — ə,—ooAs 


| | 
-Ala-Cys-Ala-Thr-Cys-Asp-Gly-Phe- 


Thus thioredoxin reductase, like its substrate thioredoxin, contains 
only two amino acid residues between the two half cystine residues and 
in both proteins the disulfide defines a small loop of only 14 atoms. 

Thioredoxin reductases have also been purified from yeast (43, 734) 
and from rat liver (735). Thioredoxin reductase from yeast is a flavo- 
protein with a molecular weight of approximately 75,000 and consists 
of two subunits, each containing one molecule of FAD. Although the 
amino acid composition of this thioredoxin reductase is quite different 
from that of E. coli, both enzymes contain 5 half-cystine residues and 
have almost identical absorption spectra. Like the E. coli enzyme, 
thioredoxin reductase from yeast is completely inhibited by p-chloro- 
mercuriphenylsulfonate (PCMS) only in the presence of NADPH suggest- 
ing that the yeast enzyme also contains a disulfide bridge at the catalytic 
site. 

Thioredoxin reductase from calf liver was found to reduce 5,5’- 
dithiobis-(2-nitrobenzoic acid) (DTNB) at the expense of NADPH even 
in the absence of thioredoxin as an intermediate electron carrier, indi- 
cating that this mammalian reductase may have a wider substrate 
specificity than the corresponding enzymes from E. coli and yeast. 

The thioredoxin reductases are highly specific for their thioredoxin 
substrate. For instance, yeast thioredoxin reductase is specific for the 
reduction of the homologous thioredoxins and does not interact with 
thioredoxins from E. coli or from T 4-infected E. coli. Furthermore E. colt 
thioredoxin reductase does not use the thioredoxins from Novikoff 
hepatoma, yeast, and L. leichmannii. 


VI. The Interaction of Ribonucleotide Reductases with 
Nucleotide- and Coenzyme Analogs 
A. Interaction with Nucleotide Analogs 


The ribonucleotide reductases appear to be highly specific for the 
B-D-ribofuranoside moiety of their nucleotide substrates. Neither the 
reductase from E. col nor the enzyme from L. leichmannii is able to 
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reduce the appropriate nucleotide of 3’-deoxyadenosine (736, 737). 
Whether the 3’-hydroxyl group is necessary for binding to the enzyme 
or whether it participates in the catalytic process has not yet been establi- 
shed. Analogs of ATP and CTP which contain 2-0-methylribose, 3-0- 
methylribose, arabinose, xylose (79) 2-deoxy-2-aminoribose, 2-deoxy-2- 
azidoribose and 2-deoxy-2-chlororibose (47) are not substrates even in 
the presence of a specific effector. Furthermore most of these nucleotides 
are unable to compete with ATP for the catalytic site. 

Requirements for the interaction of the heterocyclic base of the ribo- 
nucleotide with the catalytic site are much less stringent. Kinetic studies 
have indicated that all ribonucleotides bind to the same substrate site, 
which therefore is able to accommodate a variety of bases with different 
functional groups; specificity is mediated through allosteric modification 
of the substrate site. In addition to the common ribonucleotides ATP, 
GTP, ITP, CTP and UTP, the reductase from L. leichmannii is able 
to reduce the triphosphates of 1-ß-D-ribofuranosylbenzimidazole and 
of the antibiotics tubercidin, toyocamycin and sangivamycin (Fig. 3). 
Furthermore the rate of reduction of these four nucleotides is greatly 
increased by the allosteric effector dGTP (736, 737). In the ribonucleotide 
reductase system from E. coli the diphosphates of tubercidin and toyo- 
camycin, but not sangivamycin 5’-diphosphate, are substrates. The 
reduction of these two substrates is also specifically stimulated by dGTP, 
the allosteric effector of ADP reduction. These results suggest that the 
pyrimidine moiety and N-7 of the imidazole ring of the adenine nucleo- 
tides are not essential for interaction with the reductases. In contrast 
formycinyl triphosphate is not a substrate for the L. leichmannii enzyme 
(772). The inability of formycinyl triphosphate to act as a substrate 1s 
not surprising because the modification of a purine to a pyrazolopyrimi- 
dine base brings about quite drastic changes in the conformation of the 
entire nucleotide. This base modification affects not only the conforma- 
tion about the glycosyl bond but also the sugar phosphate backbone 
(738). 

Reduction of CTP by L. leichmannii: ribonucleotide reductase is 
specifically stimulated by its allosteric effector dATP. The structural 
specificity of the effector sites was studied with the same series of nucleo- 
tide analogs. Like 2’-dATP, 3’-dATP stimulated the reduction of CTP 
indicating that the 3’-hydroxyl group is not necessary for binding at an 
effector site. The pyrrolopyrimidine analogs, tubercidin- and toyocamy- 
cin triphosphate were also able to activate CTP reduction. Since these two 
nucleotides are also substrates for the L. leichmannii enzyme, they may 
interact with the regulatory site as the corresponding 2’-deoxyribonu- 
cleotides. Similarily in the E. cols ribonucleotide reductase system, 
3’-dATP, as well as the two pyrrolopyrimidine triphosphates mimicked 
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Fig. 3. Structures of nucleoside antibiotics 


the effect of 2’-dATP in inhibiting the reduction of 


do not act as allosteric modifiers, nor do they inhibit 
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all four common 
nucleoside diphosphates. Both 2” and 3’-O-methyl ATP stimulate CTP 
reduction, but do not promote hydrogen exchange between adenosylco- 
balamin and the solvent. In contrast ara- and xylo-ATP and ara-CTP 
the action of the 
physiological effectors, dATP and dCTP. However, ava-ATP and ara- 
CTP, but not xylo-ATP are able to promote the hydrogen exchange 
reaction. The results suggested that the reductase from L. leichmannti 
has two effector sites for the interaction with these nucleotides, one site 
which modulates enzyme-substrate interaction and another which 
affects the interaction between the enzyme and coenzyme (79). 
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The effect of ara-ATP and ara-CTP on ribonucleotide reduction has 
also been studied in crude extracts from Novikoff (739) and Ehrlich (740) 
ascites tumor cells. Inhibition by these nucleotide analogs appeared to be 
competitive with effector nucleotides, but not with nucleoside diphos- 
phate substrates. 

Irreversible inhibition of ADP and CDP reduction in Ehrlich tumor 
cells was caused by ATP-, AMP- or adenosine oxidized with sodium 
periodate. The inhibition probably involves Schiff base formation 
between the 2’, 3’-dialdehyde product and a lysine residue(s), because 
pyridoxal phosphate also inhibits ADP and CDP reduction (747, 742). 


B. Analogs of Adenosylcobalamin 


Numerous analogs of adenosylcobalamin have been tested for their 
ability to replace or to inhibit the action of the coenzyme in the adenosyl]- 
cobalamin-dependent ribonucleotide reductase reaction; the enzyme 
from L. leichmannii has been used in most of these studies. Kinetic 
studies have been used in most investigations of analog-enzyme inter- 
actions and thus the interpretation of data regarding the affinity of 
analogs for the reductase is subject to the limitations imposed on kinetic 
studies of a complex reaction. 

Three classes of analogs have been studied: a) analogs in which the 
5,6-dimethylbenzimidazole moiety of the nucleotide has been replaced, 
b) analogs with modifications on the periphery of the corrin ring and c) 
analogs in which the 5’-deoxyadenosyl moiety has been altered or 
replaced. Coenzymes in which the 5,6-dimethylbenzimidazole moiety 
has been replaced by benzimidazole, 5-methoxybenzimidazole, adenine 
and 2-methyladenine have been tested with ribonucleotide reductase. 
Adenosylcobamides with benzimidazole and adenine function as coen- 
zymes in the L. leichmannii system (52, 80), while those with 5-methoxy- 
benzimidazole and 2-methyladenine are active with the R. meliloti 
reductase (82). In contrast adenosylcobinamide, which lacks the entire 
nucleotide moiety does probably not function as a coenzyme because it 
did not yield the rapidly appearing paramagnetic intermediate (77). 

Studies with analogs modified in the corrin ring have shown that the 
generation of a negative charge on the periphery of the ring is detri- 
mental to coenzymatic activity. Adenosylcobalamin b-, d- and e-carbox- 
ylic acids show coenzyme activity only at very high concentration and 
are poor inhibitors of ribonucleotide reduction (72, 743). Blocking of the 
negative charge of adenosylcobalamin e-carboxylic acid by amidation, 
even with the bulky 2,4-dinitroaniline moiety, restores coenzyme 
activity. 
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Adenosylepicobalamin, in which the e-propionamide side chain is 
projected up rather than down relative to the plane of the corrin ring is 
completely inactive as a coenzyme, but does inhibit the action of adeno- 
sylcobalamin (744). X-ray analyses of cyanoepicobalamin have indicated 
that the e-propionamide side chain prevents the 5’-deoxyadenosyl 
moiety from occupying its normal position above the corrin ring (745). 

The cobalt-free corrinoids discovered by Toohey (746) do not show 
any affinity for the reductase of L. leichmannii (75). 

Cob(II)alamin was found to be a competitive inhibitor with respect 
to adenosylcobalamin of ATP reduction by the L. leichmannii enzyme. 
In the presence of 5’-deoxyadenosine the apparent Kr was decreased 
more than 10-fold; adenosine was able to enhance inhibition by cob(ID- 
alamin to a lesser extent, but 2’-deoxyadenosine, 3’-deoxyadenosine, 
and 4’, 5’-didehydro-5’-deoxyadenosine were much less effective. These 
studies indicate that cobalamin and nucleoside bind at sites that are 
occupied by adenosylcobalamin during catalysis and that the bind- 
ing of the cobalamin moiety is greatly enhanced by the presence of a 
specific nucleoside (706). 

Cyano-, hydroxo-, methylcobalamin, or their e-amidated derivatives 
are ineffective as coenzymes, but all exhibit some inhibition at high 
concentration. Alterations of the adenosyl moiety in the upper coordina- 
tion position of the coenzyme either diminishes or abolishes coenzyme 
activity (75, 747). Compared to adenosylcobalamin (relative activity, 
100) analogs containing the following nucleosides have the indicated 
activities: 3-ß-D-ribofuranosyladenine, 60; 9-f-p-ribofuranosylpurine, 
10; 2’-deoxyadenosine, 5; tubercidin, 3; L-adenosine, 1; and ara-adeno- 
sine, 0.6. In contrast analogs containing 1-f-D-ribofuranosylbenzimi- 
dazole, inosine, aristeromycin (Fig. 3), 9-(5’,6’-dideoxy-f-D-ribo-hepto- 
furanosyl)adenine, thymidine, uridine, cytidine and 9-(w-alkyl) adenines 
with 2 to 6 methylene groups are totally inactive (748). However all 
these analogs are inhibitors of L. leichmannii ribonucleotide reductase. 
The results suggest that although neither the 6-amino group nor N-7 of 
the purine ring are essential for activity, they nevertheless participate 
in the enzyme-coenzyme interaction. The results also indicate that even 
small changes in the ribofuranose ring have a profound effect on coen- 
zyme function. The most potent inhibitors of adenosylcobalamin action 
(Kx < Km) are the analogs derived from 2’-deoxyadenosine, 9-{5’,6’- 
dideoxy- B-D-ribo-heptofuranosyl) adenine and 9-(5’-hydroxypentyl) 
adenine. 

Two recently synthesized analogs of adenosylcobalamin, formycinyl- 
cobalamin and N&-ethenoadenosylcobalamin may be useful in studies 
of the mechanism of action of the coenzyme (749). The fluorescence of 
both nucleosides is completely quenched in the coenzyme form, but after 
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homolytic cleavage of the carbon-cobalt bond of these analogs by photo- 
lysis fluorescent nucleosides are regenerated. 


C. Inhibitors of Non-Heme Iron Ribonucleotide Reductases 


Several compounds capable of metal chelation inhibit the non-heme 
iron ribonucleotide reductases; with many of these chelators the mechan- 
ism of inhibition appears to be more complex than simple competition 
for the iron atom required by the enzyme. The most potent members 
of one such group of inhibitors are the thiosemicarbazones of 1-formy]l- 
isoquinoline (D and 2-formylpyridine (II). 


= 


II 
CH-NNHC-NH 


These compounds cause 50% inhibition of a partially purified 
reductase from Novikoff rat tumor at 10~8 to 10-? M (750, 151). Approx- 
imately the same degree of inhibition was observed with other mamma- 
Han reductases (752, 153), but the non-heme iron containing reductase 
from E. coli was not affected. The inhibition of the mammalian reduc- 
tases is only partially reversible (754). Since these compounds are 
strong metal chelators complexation of iron is probably involved in the 
mechanism of inhibition ; however excess Fe2+ does not reverse the inhibi- 
tion, and other evidence indicates that these compounds do not act 
solely by chelating free iron from solution thus depriving the enzyme of 
a cofactor (750, 151). Kinetic studies indicate no competition with 
respect to nucleoside diphosphate substrate, nucleotide effector, or 
magnesium ions, but partial competition for the dithiol substrate was 
observed. 

A high degree of structural specificity for inhibitor-enzyme interaction 
was demonstrated in studies where several isoquinoline and pyridine 
derivatives were tested with the reductase (750—755). The partially 
competitive relationship between inhibitor and dithiol substrate and the 
specificity of inhibitor binding are consistent with a model in which an 
inhibitor-iron complex binds, or free inhibitor coordinates with holo- 
enzyme, at or adjacent to a site normally occupied by dithiol substrate. 
Studies of this class of inhibitors have not only served as a probe into the 
structure of the mammalian reductases but also as a measure of differ- 
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ence with the E. coli reductase. The lack of inhibition on the E. colt 
enzyme suggests that the iron is bound in such a way that it is inaccessible 
to the inhibitor, which might indicate a substantial structural difference 
from the mammalian enzyme near the catalytic site. 

All non-heme iron containing ribonucleotide reductases are also 
inhibited by hydroxyurea and related hydroxamates, while the adenosyl- 
cobalamin-dependent reductases are not affected (27, 756). The mhibi- 
tion by these reagents can be partially reversed by excess Fet? or 
dithiols. Reaction of ribonucleotide reductase of E. colt with [14C]hydro- 
xyurea inactivated only the B2 subunit and this inactivation was not 
reversed by removal of the radioactivity (757). Inactivation by hydroxy- 
urea does not affect the iron content of protein B2, but involves the 
destruction of the stable free radical (66, 67). Reactivation can be accom- 
plished by removal of the iron and reconstitution of apoprotein B2 
with Fet+?. Hydroxyurea has been demonstrated to be a powerful 
radical scavenger in another system (758). 
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I. Introduction 


a) Iron in Biology 


Iron is especially important in the life processes as it is directly involved 
in the transport and storage of oxygen (hemoglobin and myoglobin), 
electron transport (cytochromes and ferredoxins), the metabolism of Ha 
by dehydrogenase, the utilization of peroxide (through peroxidases) and 
several other vital metabolic functions. 

It is safe to say that with a few possible exceptions in the bacterial 
world, there would be no life without iron. Generally, the importance of 
iron in various life processes comes about because of the ability of iron 
to exist in two stable oxidation states (Fe! and Fell!) and occasionally 
as Fe! and FelV. The oxidation-reduction potentials of Fe which occur 
in various metabolic functions can vary over a wide range depending 
upon the ligands attached to the iron. Nature has been very clever in 
utilizing molecules to fulfill various needs and it is amazing how fre- 
quently iron is of central importance. 

Iron is the most abundant element in the earth (29% by mass) — 
although most of it is in the core (about 25%). The remaining 4% is in 
the earth’s crust, mostly in the form of FeO and FesOs. 

The relevant biological questions are how the iron is converted from 
mineral to protein and once in the protein what is the specialized role 
played by the iron in the function of that protein? In this chapter we 
concern ourselves with the means by which iron is obtained and secured 
for the living organism through the iron storage and transport proteins, 
and in particular how Mössbauer effect and electron spin resonance 
spectroscopy can help to understand the electronic and molecular 
structure of these proteins. 

Since its discovery, Méssbauer spectroscopy has grown into a highly 
sophisticated technology in which precise measurements (10-3 ev) are 
made and detailed information about the environment of the Méssbauer 
nucleus can be found. In the case of 57Fe, the Méssbauer technique has 
seen its greatest development so that an accurate description of the 
chemical state of an iron ion can be determined which includes the details 
of the electronic states and their energies, the symmetries of the cry- 
stalline or molecular environment and the covalent sharing of the 
valence electrons. 

In Mossbauer spectroscopy of 57Fe, the radiation is specific for the 
57Fe which eliminates the complications caused by the presence of other 
atoms. Most elements do not have a suitable isotope available for Möss- 
bauer spectroscopy. Because iron is of such paramount importance in 
the molecules of life it is indeed fortunate that 5’Fe exhibits such a good 
Mossbauer resonance. If a biological molecule contains iron, it is almost 
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certain that the biological purpose for that molecule depends on the 
iron and its chemical state. This is to say that the iron is at the site of 
the biological action. 


b} Iron Transport Compounds 


Iron transport compounds are found low on the evolutionary scale, 
in the world of microbes and bacteria where there is fierce competition 
for the available supply of nutrients. These organisms require iron in 
small, but nevertheless, vital quantities primarily for growth and res- 
piration. The compounds which serve the function of iron accumulation 
for these organisms are called siderochromes (7). The siderochromes are 
divided into two categories according to their function: sideramines which 
are growth factors and sideromycins which are antibiotics. Both sider- 
amines and sideromycins have a high affinity for ferric iron which 
ensures efficient chelation of the metal even at low iron concentrations. 
The unusual properties of these compounds are discussed in the review 
by Neilands (1). 

There are also iron transport proteins transferrin, conalbumin, and 
lactoferrin which occur in birds and mammals and have many of the 
same properties as the siderochromes described in the preceeding para- 
graph, in particular a strong affinity for ferric iron and a much weaker 
affinity for ferrous iron. There are other similarities which will be revealed 
later in this discussion. The transferrins and their characteristics are 
discussed in the review by Feeney and Kumatsu (2). 


c) Spectroscopy in Biological Materials — What can be Seen? 


Many different spectroscopic techniques have been employed in the 
study of biological materials. These methods include optical absorption, 
infrared and Raman spectroscopy, electron spin resonance (ESR), nuclear 
magnetic resonance (NMR) and more recently Méssbauer effect spectro- 
scopy. 

Each of these spectroscopic methods involves a source of radiation, 
a sample to absorb or scatter the radiation and a detector arranged as 
in Fig. 1. The counting rate (the output of the detector) is measured as 
a function of the energy of the source radiation in absorption experi- 
ments with fy variable; or there is a fixed source frequency (A» is con- 
stant) and the scattered radiation Av’ is analyzed (as in the Raman 
Effect); or the detector output is measured as some external field is 
varied (as in ESR) while the source radiation is held fixed. 

In each technique there are two levels of information to be obtained. 
The upper level information is the observation of characteristic spectral 
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Fig. 1. Schematic diagram for nearly any spectroscopic experiment in which radi- 
ation is incident on the sample to be studied. The sample either transmits or scat- 
ters the radiation which is detected as a function of energy 


features such as the Soret band in the optical spectrum of heme proteins 
or the g=4.3 resonance in the ESR spectrum of the iron transport 
compounds. Another example would be the characteristic g=1.94 
resonance observed in the iron-sulfer proteins. Each member of a class 
of materials has some spectral feature in common and it is the signature 
which identifies the classification. This qualitative information is, of 
course, extremely valuable in the classification of a species and the 
strength of this signal can be used to measure the quantity of material 
present in a sample. 

On the other hand, there is a deeper level of information to be gained 
by the understanding of each spectral profile in detail, which leads to the 
fundamental knowledge of the electronic states involved and how they 
are affected by the atomic arrangements in a molecule. This problem 
is more easily visualized perhaps by the assumption that if one knows 
the electronic wave functions in a molecule and the interaction of radia- 
tion with these electronic states, then it is certainly possible to compute 
the spectral features that should be observed in an experiment. It is the 
work of a spectroscopist to do the experiment and follow this procedure 
in reverse as far as possible to obtain the molecular arrangement or 
crystal structure which in turn determines the details of the electronic 
states. 

It is interesting that no one of the spectroscopies mentioned above 
can give complete information about the electronic or molecular structure, 
but that they provide complementary information to each other. X-ray 
structures can stand by themselves, of course, but by themselves give no 
direct information about the electronic energy levels involved in the bio- 
chemical processes as in the absorption of light by chlorophyll for 
example. 
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In this article we will discuss the electronic states that are respon- 
sible for the features observed in ESR and Mossbauer effect spectroscopy 
of the iron transport materials. These electronic states are parameterized 
by a set of coupling constants which also characterize a class of materials, 
but we then are making comparisons between the electronic states of 
molecules rather than between the spectral features. If one knows the 
electronic state, then predictions of other phenomena can be made. If 
one merely knows the spectral features, then at best one can only classify 
the compound. 


IL. Mössbauer Spectroscopy 


a) Magnetic and Electric Hyperfine Interactions 


In recent years there have appeared several reviews in which Mössbauer 
spectroscopy of various biological iron compounds are discussed. For 
most of the experimental details, the reader is referred to the papers by 
DeBrunner (3), Weissbluth (4), Bearden and Dunham (5) and Lang (6). 
The review by Lang (6) is especially recommended since it deals with the 
paramagnetic hyperfine structure (PHS) observed in so many of these 
complexes and treats several different spin configurations. There are 
several excellent references concerning instrumentation; one of the most 
recent is by Kalvius and Kankeleıt (7). 

The electronic configuration and the local environment can be 
studied through the magnetic and electric hyperfine interactions which 
determine the features of the Méssbauer spectrum of the 5’Fe nucleus. 

The magnetic hyperfine interaction between the atomic nucleus with 
spin Í and its » electrons is 


n 


P Bau > [+ TS — = afi) 1- si} 


753 rg? 
i=l 
(1) 


The sum is over the z electrons of the iron ion. The vectors L, and s; 
are quantum mechanical operators which act on the electronic states 
while I acts on the nuclear wave function. 

The first term is the interaction of the nuclear magnetic moment 
(spin I) with the current loops made by the orbiting electrons with angu- 
lar momentum l;. The second term is the dipole-dipole interaction between 
the nuclear (I) and electronic (s;) spin magnetic moments. The last term 
is the Fermi contact term which describes the interaction of the nuclear 
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moment with the s-electrons which are polarized by the unpaired elec- 
trons of the ion with spin s;. For the high spin ferric ions in the iron 
transport compounds with electronic state 6S, only the Fermi contact 
term is important and Zx can be represented by AF- S. 

The electric quadrupole interaction for °’Fe is given by 


Hq = PUT — I +1)/3 + 2 IJ), n = — 2) 


where P = oer ; Q =quadrupole moment of 57Fe (Q =0.18 x 10-24 


cm?) and Vy; = Ak is the electric field gradient (EFG) in the z direction 
which is a measure of the distortion of the electric charge cloud about 
the 57Fe nucleus. In the absence of a magnetic interaction the EFG splits 
the nuclear excited state (J =3/2) into two (my; = +3/2 and my = + 1/2) 
while the nuclear ground state (J = 1/2) is unsplit. The resulting Möss- 
bauer spectrum is then a two-line pattern (Fig. 2). 
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Fig. 2. Resonance absorption between the nuclear ground state and the first excited 
state in 57Fe. The upper diagram shows the effect of an electric quadrupole splitting 
while the lower diagram shows the magnetic splitting of the nuclear levels 
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The magnetic and electric hyperfine interactions acting together give 
much more information than the sum of the two acting independently. 
It is our goal, therefore, to create experimental conditions such that 
the magnetic interaction becomes visible in the Méssbauer spectrum; 
the quadrupole interaction will always be seen if it is present. Even in 
powder or frozen solution samples, the relative orientation of the EFG 
and magnetic axes can be determined (8). The profile of the Mossbauer 
spectrum can be related to the electronic state of the iron and its ligand 
field for many of the different electronic configurations of iron (9—14). 
As an example of the simultaneous action of the quadrupole and magnetic 
hyperfine interactions, we show in Fig. 3 the spectra of a powder sample 





Fig. 3. Experimental data for a powder sample of 5*Fe diluted in K,Co(CN), taken 
at 4,2 K in applied fields of zero and 125 Gauss parallel to the gamma beam. This 
figure shows the broadening of the absorption lines due to the angular dependence 
of the hyperfine interactions 


of K3Fe(CN)¢ [diluted in KgCo(CN)¢]. This is a low spin ferric complex 
(S =1/2) with one unpaired electron. The solid curve in the figure is 
calculated from an electronic model (75) which takes into account the 
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spin-orbit coupling and the small departures from cubic symmetry 
which are manifest in the EFG, the magnetic hyperfine interaction, and 
the g-tensor. The analysis for strict cubic symmetry of a low spin ferric 
complex is given in Ref. (76). The behavior of the Méssbauer spectrum 
under variable local symmetry has also been examined (75, 76). 

One remaining electron-nuclear interaction which is measureable 
is the isomer shift which is the centroid of the absorption spectrum. 
This is a measure of the electronic charge density at the 5?Fe nucleus, 
and can be used to group iron compounds according to their oxidation 
state. The isomer shift varies because of the shielding effects of the 
valence electrons on the s-electrons at the nucleus. 


b) Paramagnetic Hyperfine Structure (PHS) 


In paramagnetic iron molecules, the unpaired electrons bound to the 
iron can give rise to a magnetic interaction with the nucleus which splits 
the nuclear states. The resulting magnetic structure in a Mössbauer 
spectrum together with the electric quadrupole interaction gives infor- 
mation about the symmetry of the iron site and also some indication 
of the strength of the interaction with the local environment. This 
information can be obtained even from randomly oriented samples such 
as powders or frozen solutions. 

For a non-Kramers paramagnetic ion (an even number of electrons) 
in a site of non spherical symmetry, there is no electronic magnetic 
moment in the absence of an external field or a cooperative magnetic 
interaction (77). This is to say that the expectation values <s> and 
<I> are zero and therefore the expectation value of Eq. (1) is zero be- 
cause the nuclear magnetic moment has nothing to which it can couple. 
This is why there is no PHS observed in the zero field spectra of high spin 
ferrous ions (S =2) which are certainly paramagnetic. However, the 
application of an external magnetic field will create a magnetic moment 
in each electronic state, which in turn will produce a hyperfine field 
which acts on the nucleus. 

For Kramers ions (odd number of electrons) there is a non-zero mag- 
netic moment (the expectation value <S> is not zero) and magnetic 
hyperfine structure can be observed even in zero magnetic field. We will 
be concerned here with the spectra in iron transport proteins where the 
iron is in the high spin ferric (6S) configuration with five unpaired 3d 
electrons, one in each of the five 3d orbital states leaving no orbital 
angular momentum (L =0) and spin S =5/2. This is a spherically sym- 
metric electronic state which in first order gives no contribution to the 
hyperfine field from the orbital and dipole terms so that only the Fermi 
contact term is unequal to zero. The EFG from this electronic state is 
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zero as well. If the electronic spin changes direction, then so does the 
hyperfine field sensed by the 57Fe nuclear moment. When the spin 
direction changes rapidly compared to the Larmor precession of the 
57Fe nucleus, then the field sensed by the nucleus is time averaged to 
zero. The electronic spin lattice relaxation can be slowed by going to low 
temperatures to reduce the vibrations of the lattice which can flip the 
electronic spin direction. The paramagnetic iron centers can be separated 
by large distances in these macromolecules so that the spin-spin relax- 





-10 O +10 
SOURCE VELOCITY (MM/SEC) 


Fig. 4. Experimental data for ®’Fe bound to ferrichrome-A taken by Wickman 
et al. (18). These data show temperature effects due to changing spin-lattice re- 
laxation. The spins are slowed enough at the lowest temperature for there to be 
complete coupling between the electronic and nuclear magnetic moments 
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ation effects are small. In the limit where the electronic relaxation time 
t is very long compared to the Larmor precession time ty of the 57Fe 
nucleus, the nuclear magnetic moment can couple strongly to the elec- 
tronic moment giving rise to sharp magnetically split Mossbauer pat- 
terns. However, as the electronic relaxation time + is shortened, the 
Mossbauer resonances become less sharp and a very broad background 
appears. This is because the magnetic hyperfine field changes its orien- 
tation during the precession of the nucleus. At faster relaxation rates 
(t<tzy) the magnetic hyperfine field is averaged to zero so that only 
the electric quadrupole interaction is left to split the Mossbauer spectrum. 
(If there is a non-zero field gradient) One of the first examples of this 
behavior was observed in ferrichrome A (Fig. 4) an inactive sideramine, 
and the essentials of the time dependent behavior were put forth by 
Wickman et al. (78). 

If an external field is applied at low temperatures, PHS can still be 
observed even in the fast relaxation regime if the populations of the elec- 
tronic spin states are unequal (79). However, we are concerned with the 
long relaxation time limit since the iron transport compounds are often 
magnetically dilute and the spins are slowly relaxing at low temperature. 

The zero field data may be only partly resolved at low temperature in 
some cases such as ferrichrome A and mycobactin P, since the size of 
the molecule may not be large enough to separate the paramagnetic 
centers and prevent spin-spin relaxation. In such cases it is sometimes 
possible to dilute these molecules in a diamagnetic solvent such as alco- 
hol, glycerol, or DMSO to achieve the necessary separation. The heme 
proteins are large enough to achieve the necessary separation naturally. 
The application of a magnetic field can help to give better resolution in 
the spectrum (20, 9), if it can overcome the coupling between different 
paramagnetic spins. 


III. Electron Spin Resonance (ESR) 


The ESR technique has been employed for many years in the study of 
biological materials and several books dealing with ESR and its appli- 
cations are in existance (27—25). A comprehensive treatment of the 
entire field of ESR by Abragam and Bleaney (26) has recently appeared. 
The reader is encouraged to consult these references for the experimental 
details and the general discussion of spin resonance. 

The goal in this section is to discuss the information that can be ob- 
tained from ESR experiments, to compare this with the information ob- 
tained from Méssbauer experiments and present a more complete picture 
of the electronic states that exist in the iron transport compounds. 
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a) Electronic Structure 


In a conventional ESR experiment, radiation of fixed frequency > is 
sent into a paramagnetic sample where it is resonantly absorbed whenever 
the radiation encounters a pair of energy levels which are separated by 
energy Av, causing a transition in the absorbing ion from the lower 
electronic state to the higher electronic state. The conservation of angular 
momentum and of parity determine the selection rules which govern the 
transition probability and thus the intensity of the absorption in the 
spectrum. 

In order to examine the electronic energy levels, the resonance can 
be tuned in by changing the energy level separation to match the micro- 
wave quantum (hy) which is easily accomplished in paramagnetic materi- 
als by the use of an applied magnetic field which splits those energy 
levels with magnetic moments (Fig. 5). A typical ESR spectrum will 
contain one or more resonance peaks as a function of the applied field 
which can be interpreted in terms of the energy level structure of the 
paramagnetic ion. This is the information that is obtained and is para- 
meterized in terms of several constants which appear in the spin Hamil- 
tonian to be discussed in detail in a later section. One is dealing with 
microwaves (hy~ 0.3 to 2.5 cm-1) which have very small quantum ener- 
gies so that in order to observe a resonance, one usually has to start out 
with a pair of energy levels that are degenerate or at least come close 
enough under the action of a magnetic field so that the radiation Ay can 
be absorbed (see Fig. 5). 


increasing magnetic field —— 


microwave P 
quantum 


absorption 
derivative 


v field 


Fig. 5. As the applied magnetic field increases in strength, the electronic energy 
levels split at a rate given by the “‘g-value”. At resonance, the energy of the micro- 
wave quantum matches the energy separation between the states and can be ab- 
sorbed 
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For this reason, ESR is much more easily done in Kramers ions which 
have an odd number of electrons and at least a two fold degeneracy in 
the absence of an applied magnetic field (Kramers theorem). This re- 
stricts one to those paramagnetic iron materials which have half integer 
spin (1/2, 3/2, 5/2) and eliminates those iron materials with spins of 
0, 1, 2 unless special efforts are made which we won't discuss here. 


b) Comparison with Méssbauer Effect Experiments 


Méssbauer spectroscopy in 5?Fe on the other hand, has radiation whose 
absorption depends only on the presence of a nucleus (°’Fe) with the right 
set of energy levels in the absorber. Thus Mossbauer spectroscopy will 
detect the presence of iron with any electronic configuration including 
diamagnetic iron and can also determine the chemical state of the iron 
being observed. Méssbauer spectroscopy or nuclear gamma resonance is 
done by varying the frequency of the source radiation and keeping the 
external field fixed in contrast to the fixed frequency variable field ESR 
experiment. Thus one could observe a set of electronic levels at a con- 
stant field value (chosen by the experimenter) or even in zero field. This 
simplifies the analysis as we shall see later. 

ESR does have one clear cut advantage which should be mentioned 
and that is the precision with which measurements can be made. The 
ESR resonance line is much narrower and more sensitive to small changes 
in the electronic state. This has been shown in the case of transferrin (27) 
where small differences were observed in the states of the two iron sites 
which were not detectable with Méssbauer spectroscopy (28). Further- 
more, ESR can be done on any isotope of iron so that samples need not 
be enriched with 5’Fe as in the case of Mossbauer experiments. 


c) Complementary Information by ESR and Mössbauer Spectroscopy 


Since one is dealing with the same set of energy levels in each experiment 
one should look for a model consistant with each experiment. In the 
case of the iron transport compounds the experimental results are 
parameterized by the coupling coefficients D, A, u, A, P in the spin Hamil- 
tonian to be discussed in the next section. The parameter P in 5? Fe can 
be measured only with Méssbauer spectroscopy since it comes about 
through the excited nuclear state which is not available in an ESR ex- 
periment. 

The magnetic hyperfine coupling constant, A, is in principle measur- 
able in both ESR and Mossbauer effect, but will be observed only in the 
case of 5?Fe which has a nuclear magnetic moment. Since 57Fe is only 
2% naturally abundant the ESR spectrum with unenriched samples 
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will not have much of a contribution from the 5’Fe ions. It would be 
possible to use isotopically enriched iron to enhance the magnetic hyper- 
fine contribution to the ESR spectrum. The parameters A and z reflect 
the symmetry of the electronic environment and cause admixtures in 
the electronic wave functions giving rise to the features in an ESR or 
Méssbauer spectrum. The actual energy level separation is scaled by 
D which can be determined by thermal population effects in the spectrum 
but because of experimental difficulties (such as a changing cavity 
resonance with temperature in an ESR experiment) D can be found more 
easily from thermal population effects in the Mössbauer spectra. How- 
ever, if the microwave quantum is large enough so that transitions be- 
tween different Kramers doublets are possible, then it should be possible 
to measure D quite accurately with ESR. The mixing parameters A 
and u can also be determined from Mossbauer experiments but not to 
such a high accuracy as in ESR. 


IV. Spin Hamiltonian with § = 5/2 


a) Description of the Interactions 


The detailed structure of the Mössbauer and ESR spectra for the iron 
transport compounds can be described in terms of a spin Hamiltonian 
H with effective spin S =5/2 for the high spin ferric ion. These para- 
meters can give information about the site symmetry of the iron. Since 
there is no orbital angular momentum in the 65 state, the effective spin 
is the same as the real spin of the iron ion. There is spherical symmetry 
in the 65 state to first order, but spin-orbit coupling to excited (non- 
spherical) orbital states gives rise to asymmetries about the iron site 
which are reflected in the spin Hamiltonian. The general form of the spin 
Hamiltonian which we will use here is a quantum mechanical operator 
X which acts on the electronic states |ms> and nuclear states |mr>. 


5 
æ =D\si—> Male “(se it Ge sl — 3) a 


4+ AI-S+28,H-S + PLL — t AB] —pybyH-I (3) 


The terms in D, A and u give the zero field splitting of the six spin 
states of the 65 state. The quadratic interaction described by D and A 
usually dominates except when the applied field makes the Zeeman term 
larger (Fig. 6). In several of the iron transport compounds, the zero 
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field splitting is comparable to the Zeeman energy in a field of a few 
kilo Gauss so that there may be considerable mixing of the states and 
perturbation theory is not applicable. Usually the quartic term has been 
neglected (u =0) but there have recently appeared some cases which 
are strongly dependent on u (29, 30). (Note: A=E/D and u =a/D where 
E and a are the coupling constants used by other authors.) 

This splitting in zero field reflects the symmetry of the local crys- 
talline environment and it arises out of the spin-orbit interaction which 


EZD = 0.33 
f 0.62 9.69 
— 
3.53 D 


4.286 4.286 4.286 


3.53 D 





Fig. 6a and b. A diagram showing the behavior of the three Kramers doublets tn an 
applied magnetic field for (a) axial (E/D =0) and (b) rhombic (£/D = 1/3) sym- 
metry. The ESR spectra would look quite different for the two cases 


couples the ground (®S) state with higher orbital states which in turn 
are more sensitive to distortions in the local environment. The origin 
of the zero field splitting in S-state ions has been studied (37) rather 
extensively but there is still no agreement about the specific origin of 
D, A and u. 

The term with u as a coefficient represents the cubic crystal field 
energy and arises from coupling with higher excited states (32). This 
contribution can usually be neglected except for those ions with nearly 
cubic symmetry. If there is rhombic symmetry then the (7, &, ¢) axes 
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coincide with the (x, y, z) axes of the zero field splitting. However, in 
trigonal symmetry the z-axis is in the (777) direction with respect to 
to the (n, £, Ë) axes. 

It might be worthwhile to point out that in the wide variety of 
crystal fields possible, the symmetry at the metal site can be orthorhom- 
bic which has fewer symmetry operations than an axially symmetric 
complex for example. However, the strength of the axial distortion might 
disturb the local environment more severely (as in the heme complexes) 
in contrast to the weak rhombic departures from cubic symmetry in the 
sideramines discussed later. The electronic wave functions which deter- 
mine the hyperfine interactions depend on the parameters A and u 
which reflect the degree of asymmetry in the local environment. 

The Zeeman interaction 2 BH-S is isotropic for the ®S state and 
describes how the energy levels behave in a small external field. Small 
corrections due to admixtures of excited states can destroy the isotropy. 

similarly, the magnetic hyperfine interaction [Eq. (1)] can be written 
as AI-S if the A tensor is isotropic although anisotropy can arise from 
non-zero orbital and dipole contributions due to admixtures with higher 
electronic states. Only small departures from isotropy have been found 
in A for one of iron transport compounds (28) with the spherically sym- 
metric §S state. 

The electric field gradient (EFG) is zero in the êS state with spherical 
symmetry, but of course, there are contributions from the charge asym- 
metries about the iron site and from the anisotropic admixtures mention- 
ed above. Since the EFG has its origin in the same asymmetry as D 
and A, the asymmetry parameter n can be constrained by n=34 (33). 
The last term in Eq. (3) is the Zeeman interaction of the nuclear magnetic 
moment with the externally applied field and is negligible for fields less 
than a few kilo Oersted. Mössbauer data in powder samples or frozen 
solutions can give valuable information because the simultaneous action 
of EFG and magnetic hyperfine coupling yield a spectral profile that 
can be interpreted in terms of the electronic states of the iron complex. 


b) Tetragonal Symmetry (Ax 0, u~ 0) 


In the heme proteins, there is a large axial crystal field which gives a 
value of D (~5 to 10 cm!) large enough to dominate all the other terms 
in the spin Hamiltonian. These give Méssbauer spectra such as those of 
Lang (34) for cytochrome-c peroxidase fluoride, in which only the ground 
state is populated at helium temperature. ESR experiments have been 
done on several heme complexes (35) and the g-values are interpreted 
in terms of the energy level diagram of Fig. 6a. The strong anisotropy 
in the g-values is indicative of the strong axial ligand field experienced 
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at the iron site. The quartic term in the spin Hamiltonian is quite 
negligible in these cases (u =0). 


c) Rhombic Symmetry (A >0, u#0, z = <100 >) 


In each of the iron transport proteins there is a rhombic crystal field 
which is strong enough to make E comparable to D, (Actually, D is rather 
weak in these cases.) In fact E = D/3 or 4 = 1/3 is the condition for which 
the extreme rhombic asymmetry exists (33, 36). A value for A>1/3 
merely implies that the principal axis of the zero field splitting is no 
longer the z axis. 

Of the Kramers doublets in extreme rhombic symmetry (4 = 1/3) 
Fig. 6b, the middle doublet has an isotropic g-value (g = 4.3) which means 
that the magnetic moment in those states will follow the field in any 
direction and so the contribution of that doublet to the Mössbauer 
spectrum indicates the hyperfine field to be antiparallel to the applied 
field. The g = 4.3 resonance in ESR is found in each of the iron complexes 
considered here. However, the magnetic moment of the lowest doublet 
will continue to point in the y direction unless the field is almost directly 
along the x or z axes and the moment of the highest doublet will likewise 
stay in the z direction; unless the applied field strength becomes compar- 
able to the zero field splitting energy and causes mixing of the states. 
The “g” values of these states which describe the rate at which a ficti- 
tious spin 1/2 would split in an applied field are (with A = 1/3), 


Es gz= gy ~ 06, gz = 9.6 (highly anisotropic) 
E o gz = gv = gz = 4.286 (isotropic) 
Ey gz = gz = 0.6, gy = 9.6 (highly anisotropic) 


A value of ¿> 1/3 implies that the principal distortion is along the 
y-axis as explained by Blumberg (36). The inclusion of the quartic term 
in the spin Hamiltonian doesn’t change his argument since the quartic 
term has cubic symmetry. We will transform all reported values for A 
and D so as to compare like quantities. 

Although it has been the usual practice to neglect the quartic term in 
the spin Hamiltonian it has been found recently that it cannot be neg- 
lected, especially in these iron transport complexes in which the mag- 
nitude of D is small and more nearly comparable to a. If one includes 
the quartic term, rather striking effects can be seen in the Méssbauer 
spectrum and although another parameter should, of course, allow a 
better fit to the data, this is not to say that the interaction does not 
exist. In order to carry out a calculation one must decide whether the 
z axis (of the EFG and of DS) is in the <100> or <111> direction 
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relative to the £, 7, ¢ axes of the cubic term. The axes (š, n, ¢) are usually 
taken as the metal-ligand directions, but there is no reason a priori not 
to choose another set as long as it has cubic symmetry. 

If we first consider z to be in the <100> direction, we then have a 
situation with orthorhombic symmetry where x, y, z are the symmetry 
axes. It has been found that the inclusion of the quartic term with 
coefficient u in the spin Hamiltonian allows a range of values for A= E D, 
and still maintains an isotropic g=4.3. I one demands that the g’s in 
one Kramers doublet be isotropic, then it follows that g = 4.3 (actually 
g = 4.286) for that doublet and a linear relationship [Eq. (4)] between 
D, E and a=uD follows from this demand (29, 30). 


{4a = 3D — 9E} or {4u =3— 94} (4) 


Thus, we can have an isotropic g = 4.3 even for near axial symmetry 
(A--0) as long as u =3/4. Over the years, g =4.3 has come to be known 
as the signature for extreme “Rhombic’”’ iron complexes, but now it is 
clear that this is not necessarily so. Several experimental examples such 
as enterobactin and mycobactin-P have shown this to be true. However, 
the pure axial case with E=0 will still give the electronic states as 
described in Fig. 6a. 

If (A=1/3, u=0) then the energy levels are separated by a constant 
amount (3.53D) as shown in Fig. 6b. However, if A# 1/3 and u is con- 
strained by Eq. (4), the isotropic g =4.3 is maintained (in the middle 
Kramers doublet) but the energy level separation is varied as A changes. 


d) Trigonal Symmetry (A=0, #40, z= <111>) 


If we have z in the <111> direction (trigonal symmetry) and if we 
again demand that g be isotropic in one doublet, we find that g=3.3 
and a = (—9)D. In this case we begin with trigonal symmetry (A = 0). 
This clearly is no help to us since experimentally we find g = 4.8 in ESR 
and the Mössbauer results also indicate rhombic symmetry. Since several 
of the iron complexes discussed below are thought to have nominal 
trigonal symmetry (based on other experiments) we have a bit of a 
dilemma. 


V. Calculation of Spectra from the Spin Hamiltonian Model 


a) Calculated ESR Spectra for S =5/2 


We begin with the spin Hamiltonian in Eq. (3) with the omission of the 
terms involving the nuclear spin (T) leaving only those interactions which 
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determine the electronic behavior. Since an ESR experiment is done in 
an applied magnetic field of at least a few hundred Gauss, the nuclear 
moment with an effective field strength of about ten Gauss at the site 
of the electron is negligible compared to the external field. 

To calculate an ESR spectrum we first of all assume the applied 
field to be in some direction (6, ó) with respect to the crystal or molecular 
axes. The matrix elements of the spin Hamiltonian are loaded into a 
(6 x 6) Hermitian matrix (Tables 1 or 2) which is diagonalized numeri- 


Table 1. The spin Hamiltonian in orthorhombic symmetry (z = <001>) 


| 5/2> | 3/2> | 1/2> las Jas jae 


10 I —_ y5 
a — 0 10 A 0 a 0 
3 7 9h | 2 
2 3 x" y5 
ie En oe 0 31/22 0 — 
g, V z” 
8 z 
E ph 0 3/24 Q 
8 = 
2 3 
— E 0 
3 z” 
Waa 
s 2" 
Table 2. The spin Hamiltonian in trigonal symmetry (z = <111>) 
| 5/2> | 3/2> | 1/2> | — 1/2> | — 3/2> | — 5/2> 
10 1 aS 2/5 
See 0 10A — 0 0 
sig y z” 
2 = 
Ta +4 0 3/24 0 0 
8 2 u 215 
yo E 0 312 À —— 
3 i 3 
8 2 
= _. — _— 0 10A 
3 3” Y 
2 
—_ 0 
3 + p 
10 1 
3 3” 


All matrix elements are in units of D. (A= E/D and u=a/D). 
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cally on a computer to find the energies and wave functions. Then the 
transition probability for the absorption of radiation between any 
pair of these levels is computed, a lineshape function with width T' is 
folded in and the contribution to the absorption of the radiation hy 
at that field strength and at that field direction is determined. This calcula- 
tion is repeated for many different directions over the unit sphere 
(although rhombic symmetry allows us to do only one octant of the unit 
sphere since all octants are equivalent) and all contributions are added 
with equal weights to simulate the spectrum for a randomly oriented 
absorber. This angular average is not appropriate, of course, if one has 
single crystal absorbers. At this point the absorption of hy at one field 
value has been calculated. The next step is to increase the applied field 
strength and repeat everything from the beginning. It is somewhat 
unfortunate that the ESR is very sensitive to small changes in the 
direction of the applied field because this makes it necessary to do the 
powder averages over one octant of the unit sphere on a much finer 
grid (i.e. to do separate computations for small changes in the field 
direction). In some of the examples studied, a linewidth (full width at 
half max) of 100 G was used and a grid of 400 field directions distributed 
uniformly over one octant of the unit sphere was necessary to get rid 
of computational artifacts. This makes the calculation very costly in 
terms of computer time. Sucha calculation with 25 different field strengths 
takes 47 min of CPU time on the CMU 1108 computer. 

An X-band (hy = 0.3 cm-1) spectrum computed from 500 G to 4300 G 
is shown in Fig. 7. Because we are dealing with iron complexes in rhombic 
environments, there are few shortcuts that one can take to decrease the 
computation time. 

The profile of the calculated ESR spectrum depends strongly on the 
spin Hamiltonian parameters D, A, and u which are adjusted to maximize 
the agreement with experiment. In the dashed curve of Fig. 7 we see 
the influence of the parameter u=a/D in the g =4.3 region. 

It is found that the graininess in the calculated ESR spectra is a 
problem if the derivatives of the absorption spectra are compared. For 
this reason, it would be more economical (given the cost of computing 
time) to compare the absorption spectra, even though this is in contrast 
to the usual experimental arrangement in ESR experiments. 


b) Calculated Mössbauer Spectra 


We assume that the electronic spin relaxation is very slow so that we 
don’t consider time dependent effects in the spectrum. This can be 
arranged experimentally by magnetic dilution and low temperatures. 
The electronic basis states are the six states |n, > =|—5/2>, |--3/2>, 


77 


W. T. Oosterhuis 


D= O.34 cm, E=0.090cm"', a=0.050 cm’! 


D=034 cm’, E=0090 cm'!, a = 0000 


absorption spectrum 


derivative spectrum 





APPLIED MAGNETIC FIELD (kOe) 


Fig. 7. Calculated ESR spectrum for a randomly oriented absorber at X-band 
(hy =0.3 cm!) with D=0.34 cm}, A= EjD = 0.265, u= aJD = 0 and with a 
linewidth of 100 Gauss folded in using a Lorentzian line shape. The dashed curve 
shows the effect of the fourth order term in the g = 4.3 region. The intensity of the 
dashed curve is reduced on the figure by a factor of two 


---, |-+5/2>. These together with the nuclear basis states [mr> form 
product states |ms>>|m;> and the spin Hamiltonian operates on these 
product states. The spin Hamiltonian matrix which results from Eq. (3) 
is given in Table 3 for the nuclear ground state. The corresponding 
24 x 24 matrix for the nuclear excited state is similar but not explicitly 
given here. The various terms mix the basis states and the (zm Xn) 
matrix is diagonalized to find the energies and eigenstates. [n = (2S + 1) 
(2I +1). 

This procedure is followed for the nuclear excited state (7 =3/2) and 
the ground state (I = 1/2) and then the 288 transition probabilities and 
line positions are calculated with a Lorentzian line shape of width F 
folded in to give the Méssbauer spectrum. This has been done on the 
computer assuming the applied magnetic field and the gamma beam in 
some fixed direction relative to the crystal axes. Then this calculation 
is repeated for all other directions of the field and gamma beam and the 
spectra are added with equal weights to simulate the spectrum for a 
randomly oriented absorber. It generally takes about ten minutes 
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Fig. 8. Mössbauer spectra for S?Fe bound to mycobactin-P in a frozen methanol 
solution at T = 1.5 K. The data in zero field and 1.3 kOe are shown. The solid curve 
is calculated using D = 0.22 cm-1, A=0.45, u =—0.26, Aex = 1.51 mm/sec and 
P = —0.20 mmjsec 


for the computer to generate one spectrum with 100 different orienta- 
tions of the magnetic field over one octant of the unit sphere. 

The profile of the resulting spectrum depends on the spin Hamil- 
tonian parameters which are adjusted to maximize the agreement 
between the calculated spectrum and the experiment. 

Mossbauer spectra for mycobactin P, a growth factor (sideramine) for 
mycobacterium Phlei, which is slightly soluble in ethanol were observed 
at 4.2 K (Fig. 8), and the spin Hamiltonian parameters D, A, u, A and P 
were found to give a reasonably good simulation of the data (37). Notice 
the effect of the applied field in stablizing the spin states as evidenced 
by the decrease in the broad background of the zero field spectrum. 

In the slow relaxation situation, there is one Méssbauer spectrum 
associated with each electronic state and the spectra due to the various 
electronic states are superimposed on each other with weights given by 
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Table 3. Explicit matrix elements for the S = 5/2 electronic state and the T = 1/2 
nuclear ground state <m,sm;|4# |msmy> 


| 5/2, 1/2> | 5/2,—1/2> | 3/2, 1/2> |32%,—1/2> | 1/2, 1/2> | 1/2, —1/2> 





5A/4+10D/3 Hp- V5$H_- 0 V10E 0 
+ af2 
+5He + Hnz 
—54/4+10D/3+a/2 //5A/2 /5BH- 0 /10E 
+ SH ez—H nz 
3A /4—2D/3—3a/2 Hnr- VS$H- 0 
+3Hez + Hnz 
—3A /4—2D/3—3a/2 /2A VsS$H- 
+ 3Hez + yz 
Where: 
B H, = He A/4—8D/3 +a + Hez + Hnz Hn- 
H_=H,—iH, —4/4—8D/3+ 
EnßnH- = Hn- 











their Boltzmann factors. Thus one can get significant changes in the 
Mössbauer spectrum by depopulating some of the electronic states at 
low temperatures. An example of this is in enterobactin shown in Fig. 9. 
If this applied field is small enough so that there is no significant mixing 
of the electronic states by the applied field, then each member of a 
Kramers doublet will yield the same Méssbauer spectrum although 
different Kramers doublets give different spectra. 

The Mossbauer spectra calculated for each Kramers doublet assum- 
ing A=? =0.46 and u= —0.27 is shown in Fig. 9 where there is a 
small field of 1.3 kOe assumed to be perpendicular to the gamma beam. 
Good agreement with experiment is obtained. The reader can observe 
the different contributions from each Kramers doublet in the composite 
spectrum. Each contribution is weighted by the appropriate Boltzmann 
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32.122] 1 os] 32, 15. =S 12 | 52 1/2 S52 12 


0 0 V5a/2 0 0 0 
0 0 0 V5a/2 0 0 
3VZE 0 0 0 V5a/2 0 
0 3/2E 0 0 0 V5.«/2 
38H_ 0 3/2E 0 0 0 
3412 3ßH_ 0 3/2E 0 0 
a—H a + Hre Hn- 88H- 0 V10E 0 

— 3A /4—2D/3—3a/2—3H2+Hnz Hn- V5ßH_ 0 

3A /4—2D/3—3a/2—3He2,—Hnz \5A/2 V5ßH_ 


—54/4+10D/3 +a]2—5He + Hne Hn- 


54/4+10D/3 + a/2-—-5Hee—Hnz 


factor. Similar spectra are obtained by Bock and Lang (38) for defer- 
oxamine. 

In Fig. 10 are presented some computer generated Mössbauer spectra 
with different degrees of rhombic asymmetry (å < 1/3 and u =0) so that 
the effect of the variation of A can be seen. Both Mössbauer effect and 
ESR show easily measurable differences between spectra calculated 
from models with z = <100> and z= <111>. The Mössbauer spectra 
are really quite different for the two models as seen in Fig. 11 and the 
ESR spectra give an isotropic g =3.3 for the trigonal case (if u =—'9) 
and isotropic g =4.3 for the rhombic case [if u =0.75 (1—34)]. Thus we 
have an experimental method to distinguish between these two possib- 
ilities. It is striking that the Mössbauer data in the cases ferrichrome A, 
enterobactin and benzohydroxamic acid are best fit with z in <100> 
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Fig. 9. Méssbauer data for a methanol solution of 5?Fe bound to enterobactin in an 
applied field of 1.3 kOe at 4.2, 10.2, 18.9 K. The solid curves are calculated from the 
electronic model with rhombic symmetry and the parameters given in Table 4. 
The contributions from the three Kramers doublets are given at the top in des- 
cending order of energy 
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Fig. 10. Calculated Méssbauer spectra using the rhombic model with D = 0.50 cm-1, 
Ais varied, u = 0, A = 1.5 mm/sec, P = 0.25 mm/sec, [H| = 1:3 kOe and T = 4.2 K. 
The various degrees of distortion from axial symmetry (4 = 0) to extreme rhombic 
distortion {A = 1/3) are shown 


direction and with u# 0. This would imply rhombic symmetry in the 
spin Hamiltonian and also in the site symmetry. However, other experi- 
mental evidence for enterobactin (39), ferrichrome A (40), and benzo- 
hydroxamic acid (47) indicates the molecular structure to be trigonal! 
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Fig. 11. Calculated Mössbauer spectra with the same spin Hamiltonian parameters 
D = 0.50 cm-1, A=0.0, u = 0.75, A = 1.51 mm/sec, P = — 0.15 mm/sec in rhombic 
(left) and trigonal (right) symmetry. The bottom spectrum is a Boltzmann average 
of the contributions from each of the three Kramers doublets shown above. We as- 
sume for this calculation that the sample is at 4.2 K in an applied field of 1.3 kOe 
perpendicular to the gamma beam. One can observe the marked difference in the 
spectra for the trigonal and rhombic model calculations 


Thus there is a dilemma at this point: How can one get g =4.3 in ESR 
and Mössbauer spectra indicative of a rhombic environment when the 
complexes are known to be trigonal (axial)? It is impossible to get g =4.3 
from a spin Hamiltonian with trigonal symmetry; and the Mossbauer 
data don’t reflect trigonal symmetry! It may be that the relationship 
between the spin Hamiltonian symmetry and the crystal symmetry 
is not as strong as we have thought. Or perhaps the ESR and Mossbauer 
data are very sensitive to small departures from axial symmetry which are 
not picked up by the X-ray techniques. 

Many of the iron transport compounds have hydroxamic acid 
ligands which bind the ferric ion and so synthetic iron (III) hydroxamic 
acid compounds have been prepared. 
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FERRIC BENZOHYDROXAMATE 
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Fig. 12. Mössbauer data for ferric benzohydroxamate showing agreement with the 
rhombic model. The zero field calculation assumes a randomly oriented field of 
about 30 Gauss to be acting of the unpaired electrons of the iron due to neigh- 
boring nuclear moments 


Experiments on these synthetic tri-hydroxamic acids have shown the 
rhombic symmetry effects in the Mössbauer data (Fig. 12). Other 
experiments using Mössbauer spectroscopy on low symmetry sites 
have been published by Lang et al. (33) in EDTA and by Zabinski et al. 
(42) in protocatechuic acid 4, 5-oxygenase, and ESR by Peisach et al. 
(43) in rubredoxin. 

If the iron is complexed by three identical bidentate ligands, there 
are two possibilities for the ligand field as shown in Fig. 13. There is the 
trigonal (facial) arrangement which has a true threefold axis (either left 
or right handed); and there is also the “rhombic” (meridional) conforma- 
tion whose ligand field can be represented by a Hamiltonian with rhombic 
symmetry. If we assume that the electronic state has the same symmetry 
as the conformation of the iron complex then we can distinguish between 
the trigonal and the rhombic conformations using either Mössbauer or 
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TRIGONAL " RHOMBIC " 


Fig. 13. Diagrams showing how asymmetric bidentate ligands can be arranged in 
either trigonal (facial) or rhombic (meridional) symmetry about the iron 


ESR spectroscopy. The difíerences in the Màssbauer spectra for the 
rhombic case with z in the <100> direction and the trigonal case with 
z in the <111> direction are shown in Fig. 11. 


VI. Method for Analyzing Mossbauer Spectra with Paramag- 
netic Hyperfine Structure 


The usual procedure to follow in the analysis of the Mössbauer data 
for the spin Hamiltonian parameters is to go first to the data taken in a 
small applied field of several hundred Gauss. This is done in order to 
break the coupling between neighboring paramagnetic ions which give 
rise to spin-spin relaxation and also to break the (Back-Goudsmit) 
coupling between the electronic spin and the ligand nuclear magnetic 
moments (44). In the absence of an applied field, this ligand nuclear 
moment effect can have drastic effects on the Méssbauer spectrum. We 
see this in the case of transferrin (28) (Fig. 14), where the spin Hamil- 
tonian parameters (D, A, P, A) used to fit the data in an applied field 
can also be used to calculate the spectrum expected in zero applied 
field. The result of the calculation is in sharp disagreement with the 
experiment in zero field even though the spectrum from the same sample 
in an applied field gives good agreement for the applied field. 


a) Neighbor Nucleus Effects in Zero Field 


The ligand nuclear moments (J;) interact with the electronic spin which 
is covalently shared to some extent between 57Fe and the ligands. The 
proper term to add to the Hamiltonian [Eq. (3)] is DI; B; S where B; 
is a coupling tensor for the 7th ligand nucleus. 
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Fig. 14. Méssbauer data for transferrin in variable conditions and the calculation 
from the electronic model. Each spectrum on the right was calculated with 
the assumption of a small randomly oriented field acting on the electronic moment. 
The data on the right were collected in exactly zero applied field 


This could be done on a computer and has been fairly successful 
(45). However, if there are several ligand nuclear moments then the 
computation rapidly gets out of hand because the dimension of the 


matrix to be diagonalized goes as (2I ve + 1) (2S + 1) TI(2I, + 1) so that 
j 
it is useful to approximate the sum of the ligand interactions acting on 


the electron spin S, by a random magnetic field Hr with strength of 
several Gauss. Then one simply adds the term g$Hr :S to Eq. (3) and 
calculates the spectrum. This is done for varying field strengths Hg and 
indeed this does improve the agreement with experiment (44) (Fig. 14). 
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VII. Experimental Studies 


a) Ferrichrome A 


Ferrichrome A is one of the most studied iron transport compounds. 
It is found in the smut fungus Ustilago sphaerogena and has a structure 
similar to several molecules which are growth factors for several micro- 
organisms. The molecule consists of a ferric ion bound by three biden- 
tate hydroxamic acid residues which form a left-handed propeller as 
established by the X-ray crystallographic studies by Zalkin, Forrester and 
Templeton (40). There is a three-fold axis through the iron corresponding 
to the axis of the propeller. The ferric complex is a red-brown color like 
most of the natural iron complexes discussed below. The iron is bound 
very tightly as indicated by the dissociation constant K where pK = log 
K= 32. Nevertheless iron is found to exchange rapidly with free iron in 
a solution (7). 


Electron spin resonance (ESR) experiments were done on ferrichrome 
A by Wickman, Klein and Shirley (46) and a strong resonance was 
observed at g =4.3. Their results were interpreted in terms of a “rhombic” 
spin Hamiltonian in which D >0 and A=E/D =0.25 + 0.04. 


Méssbauer experiments were then done by the same authors (78) 
who were not only among the first to observe paramagnetic hyperfine 
structure and spin relaxation effects in Mössbauer spectra but also 
formulated a theory to describe the effects of time dependent hyperfine 
fields in the Méssbauer spectrum. The data in Fig. 4 show the effect of 
increasing the electronic spin relaxation time by lowering the tem- 
perature. The spin Hamiltonian parameters are found in the long relaxa- 
tion time limit to be 220 kOe/unit spin for the hyperfine field (A ~ 1.50 
mm/sec) and D~ 1 cm! estimated from thermal depopulation effects 
in the Méssbauer spectrum. The zero field data still show a broad para- 
bolic background which is indicative of the existance of relaxation 
effects even at the lowest temperature and concentration. The applica- 
tion of a small magnetic field can help to sharpen the resonance lines 
(20) and facilitates the determination of the hyperfine coupling constant 
A which is found to be 1.47 mm/sec for the nuclear excited state of °’Fe 
or 216 kOe/unit spin. A refined analysis (62) of Méssbauer data gives 
D =0.53 cm-1, 4=0.51, u =— 0.133, and P =— 0.50 mm/sec. 


Dowsing and Gibson (47) have deduced values of D =0.5 +.03 cm-1 
and 4 =0.28 + .02 from their ESR experiments. 
Recently Brackett, Richards and Caughey (48) have done far infrared 


measurements on ferrichrome A and report that D = —0.27 cm! and 
A =0.285. 
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b) Deferoxamine 


Deferoxamine (desferrioximine B) is a growth factor (sideramine) for 
several micro-organisms in which three hydroxamic acid residues 
bind ferric iron very tightly with pK ~ 30.5 (49). The iron complexing 
ability of this substance is used in cases of severe iron poisoning to chelate 
excess iron and is available for clinical use under the name “‘desferal’’. 
The ferric complex is red-brown and magnetic susceptibility measure- 
ments yield uer =5.86 Bohr magnetons at room temperature. The 
molecular complex is probably such that the hydroxamic acid ligands 
form a propeller arrangement about the iron. 

A rather complete Mössbauer study of the 5?Fe complexed by de- 
feroxamine has been done by Bock and Lang (35) and their results are 
in excellent agreement with a model spin Hamiltonian in which the 
coupling constants have the values D=0.5 cm-l, A=y7/3=0.46, 
Aex =1.474 mm/sec, P =—0.30 mm/sec (or AEg=0.77 mm/sec in 
absence of magnetic interactions). The asymmetry parameter À = 7/3 > 
1/3 merely means that the principal distortion axis is in the y-direction 
instead of the z-direction for the axis system assumed for the spin 
Hamiltonian. We will use an axis system in which 0<A<1 and D>0 
so that we can compare with other molecules. Some of the data and the 
fitted curves show the typical features associated with these iron trans- 
port compounds. 


c) Various Hydroxamic Acid Complexes of 57Fe 


The natural hydroxamic acid complexes of iron have been the central 
topic of discussion here, but since the naturally occuring compounds are 
somewhat complicated by the organic residues, it is probable that some 
valuable insight into the iron-binding abilities of the natural compounds 
can be gained through the study of some model hydroxamic acid com- 
plexes. 

Again, these ferric complexes are very stable (6K ~28) as are the na- 
tural iron chelating compounds and they bind ferrous iron relatively 
weakly (50). In an acid medium (pH ~2-6) the hydroxamic acid forms 
a deep purple 1:1 complex with ferric iron while in neutral or basic 
solution a brown-red 3:1 complex is formed in which the hydroxamates 
act as three bidentate ligands which occupy the six octahedral posi- 
tions about the iron. The 1:1 purple complex has an absorption maxi- 
mum at 5000—5200 A depending on the nature of the hydroxamic acid 
ligand. The 3:1 red-brown complex (it becomes yellow-orange on dilu- 
tion) has an absorption maximum usually between 4250—4400 A with 
many of the same spectral features as the naturally occurring complexes 
of ferric iron. 
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Mossbauer experiments on a few of the concentrated powder hydrox- 
amic acid complexes with Fe3+ were reported by Epstein and Straub (57) 
who observed only a small variation in the spectral features by compar- 
ing the different hydroxamic acid complexes. They studied complexes 
with benzo, acetyl, oley, napthal, nicotinyl, and salicyl hydroxamates. 
Epstein and Straub concluded that there is a high degree of ionic binding 
to the iron (thus a small amount of covalency) based on the isomer 
shifts which are among the largest observed for ferric complexes so far. 
This conclusion is in agreement with the observations by Spartalıan 
and Oosterhuis (29, 30) in which paramagnetic hyperfine structure was 
observed and the magnetic splitting was among the largest observed 
for any ferric ion so far, which is interpreted as being due to very little 
delocalization of the magnetic electrons. | 

Epstein and Straub note that there are two geometrical isomers 
possible for these MA3B3 complexes as shown in Fig. 13. It is possible 
to determine the differences between the spectral profiles of the facial 
or trigonal conformer and the meridional or rhombic conformer as 
discussed above. 

A crystal structure has been determined for the Fe(III) (benzohydro- 
xamate)3 complex by Lindner and Gottlicher (41) and shows the 3-fold 
propeller coordination of the ferric ion. This may be the case for the 
solid but one wonders if this coordination persists in solution or frozen 
solution since the spin Hamiltonian parameters from Mossbauer experi- 
ments seem to indicate rhombic symmetry. These parameters are listed 
in Table 4. If in fact there is a 3-fold symmetry in the molecular arrange- 
ment, then one must question the relationship between the symmetry 
of the spin Hamiltonian and the symmetry of the local crystal or ligand 
field (52). 


d) Mycobactin P 


Mycobactins are growth factors for most species of mycobacteria such 
as Mycobacterium johnei, M. phlei, M. Tuberculosis (53) and their 
production by the mycobacterium is stimulated by iron deficiency in 
the growth medium. They were originally thought to act as scavengers 
for ferric iron, but recent work by Ratledge (54) has shown that they 
are located in the lipid rich regions of the cell wall and participate in 
the passage of iron through the cell wall by the capture of Fe3+ ions 
at the outer surface and releasing the iron to the cytoplasm inside as 
Fe?t, 

The mycobactins consist of two hydroxamic acid residues which bind 
the ferric ion as well as a third bidentate ligand with a nitrogen and an 
oxygen at the iron site. Snow (53) has deduced a probable molecular 
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structure in which there is a charge associated with three of the oxygen 
ligands. The crystal structure of mycobactin-P has been determined 
recently by Hough and Rogers (55) and has verified the molecular 
arrangement proposed by Snow in which the charged oxygen ligands are 
in a meridional configuration as contrasted with the facial arrangement. 
The details of the crystal structure of mycobactin-P show the ferric 
ion to be exposed because of the strained octahedral coordination of 
five oxygens and one nitrogen. The dimensions offer some insight into 
the high stability of the ferric complex and also for the ease of release 
according to Hough and Rogers (55). 

Mycobactin-P is a molecule of molecular weight 925.9 which readily 
forms a red-brown complex ion with iron (absorption maximum at 
4450 A) which is quite soluble in chloroform and somewhat less soluble 
in alcohol. It is insoluble in water. The ferrimycobactin is extremely 
stable with pK~30. Mössbauer experiments have been done (37) on 
the frozen alcohol solution and the data were fit using the spin Hamil- 
tonian [Eq. (3)]. The best agreement with u=0 was obtained with 
D =0.34 cm-1, 4=E/D=0.27, AE =— 0.76 mm/sec, and A =221 
kOe unit spin. 

However, subsequent ESR experiments (29) have shown a rather 
narrow resonance at g = 4.3 which is incompatible with 4 =0.27 and u =0 
(Fig. 15). A value of A=0.33 would be better suited for the narrow 
(isotropic) g =4.3 resonance, but isin disagreement with the Mössbauer 
experiment. These contradictory experimental results were reconciled 
by the inclusion of the quartic term (wu 0) in the spin Hamiltonian. 
A new set of parameters was found to give agreement with both ESR 
and Mossbauer experiments. (They are listed in Table 4). The low 
symmetry of the spin Hamiltonian is consistant with the low symmetry 
molecular environment determined from the X-rays. 


APPLIED FIELD (kOe) 
080 105 130 155 180 205 





Fig. 15. ESR data at X-band for mycobactin-P showing the large unsplit resonance 
at g = 4.2 
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e) Enterobactin 


Enterobactin (MW 669) is derived from certain enteric bacteria such as 
Escherichia coli and is classified as a siderochrome by Neslands (7) in 
which the iron is bound by three bidentate catechol groups instead 
of the hydroxamic acid ligands described above. There are three “charg- 
ed” oxygens which bind the iron either in a facial or meridional con- 
figuration. 


C=O Ice-o, 
Fe | Fe 
—o c-0 
Pi I 
hydroxamic acid catechol 


Fig. 16. Catechol and hydroxamate ligand binding to ferric iron 


Enterobactin is described (7) as the prototype of the catechol sidero- 
chromes and is produced by enteric bacteria when placed in an iron- 
deficient medium. Enterobactin has a high affinity for iron pK >25 and 
is soluble in ethanol (or water). ESR experiments indicate (56) the 
characteristic g = 4.3 resonance found in the other high spin ferric com- 
pounds mentioned in this chapter which have been described by a 
“rhombic field” spin Hamiltonian and would indicate a meridional 
configuration. Proton and Carbon 13 NMR studies have been done by 
Llinas, Wilson and Neilands (39) which indicate the ferric ion to be in 
a propeller type of configuration leaving the iron on a 3-fold axis of the 
molecule. This would of course indicate the facial configuration in contrast 
to the ‘‘rhombic field” discription used to interpret the g=4.3 ESR- 
signals. 

Mossbauer effect experiments have been done on 57Fe-enterobactin 
and show unusually well resolved paramagnetic hyperfine structure and 
the spectral features have been reproduced extremely well with a spin 
Hamiltonian with rhombic symmetry (57). An attempt to fit the Möss- 
bauer data with a trigonal spin Hamiltonian corresponding to the 3-fold 
propeller arrangement fails completely. Some of the spectra and the 
computer fits are shown in Fig. 9 and the parameters are listed in Table 4. 
In this figure, the temperature dependence is due to the population of the 
excited electronic states which then make a larger contribution to the 
composite spectrum at higher temperatures. 
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f) Transterrins 


The transferrins are iron chelating agents found in the various fluids 
of vertebrates including humans (2). They are powerful iron chelators 
with binding constants for ferric iron of PK ~30-—31. They have mole- 
cular weights up to about 80,000 and two iron binding sites per molecule. 
Other transition metals can occupy the two iron sites although iron will 
displace them. 

The transferrins include serum transferrin found in blood serum, ovo 
transferrin or conalbumin found in egg whites and lacto transferrin 
found in the milk of mammals. Serum transferrin has the primary role 
of transporting iron throughout the system from the iron storage areas 
to the bone marrow for incorporation into hemoglobin in the immature 
red blood cells. It also appears to have the function of sweeping free 
iron ions out of the serum to prevent iron poisoning. 

The functions of ovotransferrin and lacto transferrin are hypothesized 
to be the inhibition of bacterial growth in the egg white or milk respect- 
ively, perhaps by limiting the available iron necessary for bacterial 
growth. 

A question of recent interest is whether the two iron sites in each 
molecule are equivalent (58), Mössbauer effect studies (28) and ESR 
studies (59) of human serum transferrin could detect no difference in the 
spectra of the two iron sites, but recently a very careful study by Aasa 
(27) has shown small differences in the ESR spectrum due to a small 
change in the ESR parameter 2. The Mossbauer data were shown to be 
in good agreement with a calculation by Spartalian and Oosterhuis (28) 
and more recently by Tsang, Boyle and Morgan (60). The spin Hamil- 
tonian parameters for each experimental study are presented in Table 4. 

The ESR experiment by Aasa (27) yields two sets of spin Hamilto- 
nian parameters — one for each iron site. A Méssbauer spectrum cal- 
culated from each set of Aasa’s parameters has been found to be no 
different from the calculation assuming identical sites. This indicates 
that the ESR experiment is much more sensitive to small effects than 
the Mössbauer technique. 

The zero field data show the effects of neighboring nuclear magnetic 
moments (Fig. 14) which are represented by a random magnetic field 
of about 4 Oersted acting on the electronic moment of the iron. Aasa’s 
experiments (27) on the transferrin of Hagfish show it to be quite differ- 
ent from that of human serum transferrin and distinctly shows the 
differences in the two iron sites. Conalbumin or ovotransferrin is similar 
to serum transferrin in that they both form pink complexes with ferric 
iron with absorption maximum at 465 mm and both bind iron tightly 
in two different sites. A recent Mössbauer study by Aisen, Lang and 
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Woodworth (61) has shown that differences in the two iron sites are 
detectable in this molecule but the data are not fit as yet. Other ESR 
studies of conalbumin (45) show the g=4.3 resonance. No ESR or 
Méssbauer studies of lacto transferrin have been reported yet. 

Each of the other iron binding agents discussed above have either 
catechol or hydroxamic acid residues. No hydroxamic acid residues are 
found in transferrin although there are most certainly some nitrogens 
as well as oxygens gripping the iron which also binds a bicarbonate 
ligand. 


g) Discussion 


Each of the natural iron binding compounds above has a very strong 
affinity for ferric iron, but relatively weak attraction for ferrous iron. 
Since these molecules most certainly give up the iron at some point in 
the biological cycle, the importance of binding is of crucial interest here. 

The mechanism by which these compounds give up the ferric ion is 
not yet understood nor is it understood whether or how the iron gains 
entrance through a cell wall or in what chemical state the iron is as it 
passes through the cell wall. One interesting application would be to 
study the mechanism by which iron enters the immature red blood cell 
for incorporation into hemoglobin. If the electronic states of the initial 
(transferrin) and final products (hemoglobin) are known then it may be 
possible to observe other chemical states or to deduce the mechanisms 
by which the iron enters the cell. 

There are many other iron complexes of a similar nature yet to be 
investigated including many model complexes such as cupferron and 
ferric acetyl acetonate as well as the hydroxamates and the catechols. 

We see in Table 4 the similarities of the spin Hamiltonian components 
among the classes of the iron proteins and transport compounds and the 
differences between the classes. Each of the natural iron binding com- 
pounds (1—15 in Table 4) has A with values ranging from 0.25 to 0.51 
and in those cases where a non-zero u was considered in the analysis 
we find u in the range from —0.13 to —0.35. Although each of these 
iron complexes has a strong, nearly isotropic ESR resonance at g =4.3 
it would appear to be too much of a coincidence for each complex to 
have the same ligand field symmetry as given by 4. The consideration 
of the fourth order terms in the spin Hamiltonian allows a wider variety 
of ligand fields and still allows the isotropic g =4.3 signal. 

It is also noteworthy that the spin Hamiltonian (and therefore the 
electronic wave functions) have rhombic symmetry although in several 
cases there have been other experiments which indicate trigonal symme- 
try in the local environment for the iron atom. 
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The heme complexes show axial symmetry in the spin Hamiltonian 
as expected and D is an order of magnitude larger than in the iron 
transport compounds. The hyperfine fields are considerably larger in 
the iron transport compounds compared to nearly all other high spin 
ferric complexes. This indicates a high degree of localization on the iron 
of the magnetic valence electrons. The quadrupole splittings are of the 
same magnitude as other high spin ferric complexes but have the extreme 
rhombic symmetry in the electric field gradient. 

The study of the natural iron-binding compounds has brought 
forth several new and interesting aspects concerning the electronic 
states of the ferric ion and has stimulated some new thinking with regard 
to the ligand fields acting on the metal ion. 

It is also interesting to speculate about the feed back system in which 
the low iron stress causes an organısm to produce much larger than 
normal quantities of a siderochrome. This question and the question 
of how a siderochrome binds to or passes through a cellular membrane 
are undoubtedly related to the molecular conformation which can be 
studied via the electronic states using Méssbauer and ESR Spectroscopy. 
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I. Introduction 


The general versatility of the Mössbauer technique in biology was de- 
scribed in detail by Debrunner (1) and Lang (2). Of the elements which are 
important in biology only a limited number is applicable as Mössbauer 
isotopes. These are K49, Fe57, Zn6”, [127, and [129. The spectroscopic 
difficulties with K4° and specifically with Zn®’, however, restrict the 
Mossbauer technique to the investigation of biological compounds con- 
taining only either iron or iodine. Despite this limitation, the first appli- 
cations of Méssbauer spectroscopy to be used in biophysics by Gonser 
et al. (3) (1963, 1965) Karger (4) (1963), and Lang et al. (5) (1966) released 
an avalanche of work dealing with “Mossbauer Spectroscopy of Biological 
Materials”. All relevant publications have been collected in an index 
under this title by May (6). 

The purpose of the present work is to describe the applicability of 
the Mossbauer effect to heme compounds, particularly with regard to 
more recent work in this field. Two earlier contributions to this series 
(Structure and Bonding) partly dealt with the same subject: M. Weiss- 
bluth (7) in Vol. 2 (1967) with “The Physics of Hemoglobin’, and A. J. 
Bearden and W. R. Dunham (8) in Vol. 8 (1970) with “Iron Electronic 
Configurations in Proteins: Studies by Mössbauer Spectroscopy”. After 
a short description of the effect and the experimental set-up (Section II), 
the study of iron heme compounds is considered (Sections III-VIII). 
Sections IV-VII represent an attempt to investigate the heme geometry 
of deoxy (Mb, Hb), oxy- (MbO2, HbOg), and carbomonoxy-(MbCO, 
HbCO) myoglobin or hemoglobin by comparing experimental and theo- 
retical Mossbauer data. From the work of Parak et al., which is described 
in Section VIII, we may learn about the difficulties and expectations 
concerning y-ray scattering experiments with myoglobin single crystals. 

For completeness the reader should refer to two review articles, by 
Lang (9) and Münck (10), which are in preparation. Both articles deal 
mainly with the interpretation of magnetically hyperfine split Méssbauer 
spectra of ferric heme proteins and other biological molecules. 

Also in the present volume of “Structure and Bonding’, W. T. 
Oosterhuis contributes to an understanding with his article on “The 
Electronic State of Iron in Some Natural Iron Transport Compounds; 
Determination of Méssbauer and ESR Spectroscopy”. 


II. Mössbauer Effect 


1. Recoilless Emission and Absorption 


Since Mössbauer’s first publication (77) on “recoilless emission or ab- 
sorption of y-rays’’ several extensive articles and books (72-24) have 
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been published on the Mössbauer effect. The application of the effect 
has been discussed (25-37) at several international conferences. A very 
illustrative description of the method was given by Frauenfelder (32) 
during the symposium “Mössbauer Spectroscopy in Biological Systems” 
at Urbana, Illinois (USA), in March 1969, which is essentially given here: 

“Consider a cannon fixed rigidly to the earth, and a target at a dis- 
tance Ey. If we trace a large number of shots and measure their distances 
E, we find that not all have flown the same distance; there exists a strag- 
gling as sketched in Fig. la. We characterize the distribution by the full 


JIE 





Fig. 1. Model for explaining the Méssbauer effect. [Taken from Ref. (32)] 
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width at half-height and we call this quantity, I’, the “natural line width”. 
J(E) is the number of projectiles observed at the distance E. Next, we 
mount our cannon on a small boat and try to hit the same target again. 
We expect the same natural line width and indeed it appears. However, 
the small boat recoils in order to conserve momentum and the shots fall 
short (Fig. 1b). The distance they fall short, R, can be calculated from 
energy and momentum conservation. Actually, the situation is usually 
not nearly as clean. Unless the lake is quiet, the distance distribution will 
be much wider than shown above. In a storm, the boat will pitch and roll, 
and the distribution will experience an additional broadening, which we 
call the “Doppler broadening” (Fig. 1c). Shots fired from a cannon 
mounted on a small boat in a rough sea experience both a recoil (R) and 
a Doppler broadening (D). Is there a way to avoid the recoil R and the 
Doppler broadening D? Rudolf Méssbauer found the answer in a some- 
what different context: If we wait till the lake is frozen, Doppler broad- 
ening and recoil will disappear!” 

With the simple picture described so far, part of the Mossbauer effect 
can be understood in a straightforward way: 

A nucleus in an excited state corresponds to the cannon and the 
nuclear y-decay to its firing. Excited nuclei of atoms can emit during 
their decay from their excited state to their ground state. The y-radia- 
tion has a specific energy, E,, and a certain energy spread (line width), 
T, which is determined by the mean lifetime 7 of the excited nuclear 
state through the Heisenberg relation I’: + ~}. For the most commonly 
used Mössbauer isotope, Fe5’, E, and T are 14.4 keV and 4.6: 10-% eV,. 
respectively. During the emission process a recoil energy Er = E2/2 Mc? 
(M = mass of the emitting nucleus, c = velocity of light, ER =2 + 10-3 eV 
for Fe57) is transferred to the nucleus, thus decreasing the energy of the 
emitted y-quantum by the same amount En because of energy conserva- 
tion. During the absorption process the y-quantum “knocks” at the 
absorbing nucleus and again loses the same amount of energy. The 
resulting total energy loss of the y-quantum during the emission—ab- 
sorption process is therefore 2 Er, and since 2 Eg =4- 10-3 eV is con- 
siderably larger than T’ = 4.6 - 10-9 eV, resonance absorption is impossible, 
provided that source and absorber have a velocity of v ~ 0 with respect 
to each other. 

The basis for an observation of the nuclear y-resonance emission and 
absorption (Méssbauer effect) is the fact that for nuclei of atoms bound 
tightly in a solid, the recoil energy can be taken up by the crystal as a 
whole under favorable conditions (‘cannon frozen in the lake” analogy). 
For this purpose the atomic force constants of the emitting or absorbing 
atom bound in the solid must be strong enough (t.e. a high Debye 
temperature 0p is required); and, furthermore, the temperature of the 


105 


A. Trautwein 


solid must be low enough to prevent excitation of lattice vibrations. 
Then a fraction f of y-quanta does not experience recoil energy losses. 
The quantity f is called the Debye-Waller factor (33) or the recoilless 
free fraction (see Table 1, first column). 

Because of the narrow energy line width of the Fe? y-radiation small 
relative displacements between source and absorber nuclear levels 
can be measured by the Mossbauer effect. The position of the nuclear 
energy levels is, influenced by electrostatic and magnetic interactions 
between the atomic electron shell and the Mössbauer nucleus. Three 
different types of interaction are of interest and will be described briefly. 


2. Isomer Shift 


The Coulomb interaction between the nuclear charge distribution and 
electronic charges, which penetrate into the nucleus, is proportional 
to the amount of electron density, e(o), within the nuclear volume. 
This electrostatic interaction alters the energy separation between 
nuclear ground state and excited state. Since o(o) is influenced by 
the chemical surrounding of the atom under study, one expects different 
nuclear level shifts for the same Mossbauer atom in different compounds. 
For example, for a metallic source of radioactive Co®? in Rhodium and 
an absorber of ferrous high-spin Fe5? in deoxygenated myoglobin, the 
different transition energies Es and E, (Table 1) give rise to a net isomer 
shift of ó =E,4~—Es, of the order +0.8 mm/sec. 


3. Electric Quadrupole Splitting 


The interaction of the nuclear quadrupole moment Q with the local 
electric field gradient (EFG) at the nuclear site causes a splitting of the 
nuclear energy levels. Table Ic illustrates the situation for an Fe®? 
absorber atom embedded in a crystal with lower than cubic point 
symmetry for the iron lattice site. Since the nuclear ground state of Fe®? 
has zero quadrupole moment it remains unsplit. The excited state of 
Fe5’, characterized by Q > 0, however, splits into two sublevels. 


4. Magnetic Hyperfine Splitting 


If a nucleus possesses magnetic dipole moments in its ground-excited 
states it can interact in both states with magnetic fields. The effective 
magnetic fields at the nuclear site, Her, may be produced either by 
applied external fields, Ho, or by intrinsic fields, Hint, due to contribu- 
tions of the electron shell. Under the action of Hett the ground state of 
the Fe5? nucleus splits into two, and the excited state into four sublevels 
(Zeeman splitting). According to selection rules for magnetic dipole 
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radiation six transitions are usually possible then from the ground state 
to the excited state in an Fe57 absorber nucleus, as indicated in Table 
id. (A more general situation, where we may expect eight transitions, 
is described in IV. 2.1 and Fig. 11). 


5. Apparatus 


A Mössbauer experiment is performed in practice (Fig. 2a) by mounting 
the radioactive source S on a velocity drive L. By applying well-known 
velocities v to the source one can modulate the energy of the emitted 





Jr (v) 1a 


ze a) 


Vr y 


Mb single crystal 


g., 
m 


N 


b) 





Js iv) fy 


Vr V 


Fig. 2. Apparatus (a) Méssbauer transmission experiment, (b) Mössbauer scattering 
experiment 
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radiation through the Doppler effect. The velocities are defined as positive 
if the source approaches the absorber. In a Mössbauer transmission ex- 
periment (Fig. 2a) the y-beam with intensity Jo passes through an 
absorber. The transmitted intensity /r(v) is counted as a function of the 
source velocity v by a detector D and is stored in the memory of a multi- 
channel analyser MA. The velocity of S and the channel number of MA 
are adjusted to each other by electronic synchronization ES. Single 
y-line sources are usually used to simplify the absorption spectra. Re- 
sonance absorption of y-rays occurs for specific velocities vr, namely 


under the resonance condition Es (1 +=) = Ea. The decrease of the 
C 


transmitted y-intensity Jr(v) by approaching resonance velocities vr 
is described by Jr(v) = Jo — Jo u É — vr)? + Cu and shown in 
wt 
Fig. 2a. 
Mossbauer scattering experiments (Fig. 2b) are characterized by 


scattered y-intensities /s(%) = / o= [e — %r)2 + (5) ] ~" In this case 
wt 


the resonantly absorbed y-quanta are reemitted isotropically. These 
reemitted (‘‘scattered’’) ys are then counted by the detector D which 
is shielded from direct irradiation of the source S by tantal (Ta) or lead 
(Pb) screening. In Section VIII, however, where a short review is given 
on the Mössbauer scattering technique for the solution of the phase 
problem for protein-structure determination, we shall be concerned 
with y-waves rather than with y-quanta, since interference phenomena 
of radiation are adequately described in the wave picture. 

More detailed descriptions of Méssbauer equipment may be found in 
the literature (78, 34, 35). 


III. Experimental Considerations 


From Section II we know that the key point of the Mossbauer effect 
is the recoilless emission and absorption of y-rays. This has the important 
consequence that the Mössbauer nuclei must be bound in a solid system, 
or equivalently, that the myoglobin or hemoglobin samples cannot be 
studied in aqueous solution. They have to be frozen, crystallized, lyophil- 
ized or, at least, they have to be embedded in a material with high 
viscosity. Because of these unphysiological conditions under which 
biological compounds are investigated one may doubt the reliability 
of experimental Mössbauer results with respect to the real physiological 
state of a molecule. This same problem arises in other experiments: 
susceptibility, ESR, X-ray diffraction, X-ray, and neutron small-angle 
scattering. There, however, one can argue: (a) that the protein structure 
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of hemoglobin in aqueous solution as inferred by X-ray and neutron 
small-angle scattering (36—38) is in agreement with the atomic coordi- 
nates determined by X-ray diffraction analysis from hemoglobin single 
crystals (39, 40), and (b) that the electronic structure of the heme iron 
is not detectably changed by going from liquid to solid state as exem- 
plified by susceptibility measurements (47). Thus, we may assume that 
the molecular or electronic structure under certain circumstances under- 
goes only undetectably small changes upon freezing or crystallizing the 
sample, As further confirmation for the validity of this assumption 
we found the quadrupole splittings of an Mb single crystal at 77 °K and 
of a frozen Mb solution at 77°K (from which the single crystal was 
originally grown at around 300°K) to be equal within the error limit 
(42). In order to conserve the “native” state of the sample it is very 
important, however, to freeze the Mb or Hb solution as fast as possible. 
Slow freezing will allow the growth of relatively large crystals of ice, 
changing the concentration and pH-value of the solution. Fast freezing 
probably prevents these changes which have been found to affect the 
Mossbauer spectra significantly (43—43). 

After cooling the sample to the temperature of liquid Na we some- 
times found it very important not to raise the temperature of the sample 
above about 200°K for a long time. After keeping an in vitro Fe5? 
enriched Mb sample for about 20 h at 220°K we found a significant 
amount of heme iron irreversibly transformed to a different electronic 
structure, indicated by different Méssbauer spectra (Fig. 3). Similar 
changes have been observed by G. Lang with Hb (46). This strange be- 
havior of a frozen myoglobin or hemoglobin sample at relatively low 
temperatures may correspond to the anomalies which have been observed 


We 





v 


Fig. 3. Móssbauer spectra of frozen Mb solution: (a) at 77 °K, (b) at 220°K after 
keeping the sample for about 20 h at this temperature 
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with high-spin ferrous and ferric ions in ice (47): at 183°K the Möss- 
bauer lines vanish, that means the f-factor which is a measure for the 
recoilless absorption of y-rays becomes zero, indicating rapid molecular 
motion in the sample; increase in the temperature causes the ice to 
recrystallize and the Mössbauer absorption lines reappear. Similar 
effects have been found in frozen solutions of iron compounds at about 
190 °K (45) (Fig. 4). Further systematic studies are necessary to elucidate 





~200 -150 -100 -50 0 °C 


Fig. 4. Temperature dependence of Mössbauer parametrs ofe frozen solution of 
(a) Fe Cle, (b) Fe (NHa)2 (SO4)2, (c) FeSO4, and (d) Fe (ClO4)a J’: line width. 
J: intensity of the absorption lines. A Eq: quadrupole splitting. [Taken from Ref. 
(54)] 
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whether the phase transition of ice or the melting of glass phases can 
affect the bio-molecule and generate denatured protein chains which 
can then be detected by Méssbauer spectroscopy through the change of 
electronic state of the heme iron. Due to this problem it might be wise 
to limit the study of temperature-dependent Méssbauer parameters of 
biological samples to a temperature of T < 200 °K. 

Upon lyophilization the conformation of a protein can change too. 
From Mössbauer investigations of Hb (human) (47, 48) and Mb (horse 
heart) (48) in frozen solution and in dehydrated from it is known that the 
high-spin ferrous state in Hb is converted to an equal mixture of ferrous 
high-spin and low-spin iron in anhydrous Hb, whereas Mb and anhy- 
drous Mb both contain ferrous high-spin iron only (Fig. 5). 


Effect [%] 0 


Hb 





AHb 








Fig. 5. Quadrupole splittings of Hb (human) at 80°K, of AHb (human) at 80 °K, 
of Mb (horse) at 83°K, and of AMb (horse) at 83°K. [Taken from Ref. (48)] 
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Further experimental problems upon Mössbauer inspection of heme 
proteins can arise because of the relatively low abundance of Fe57 
within natural iron (2.2%). We therefore give an estimate of the amount 
of biological material which is necessary to produce a feasible effect 
in a reasonable amount of time. The resonance cross-section oo for 
Fe5? is 2-10-18 cm?. We divide this value by a factor of 8 (a factor 
of 2 comes from the finite line-width of the spectrum which is at least 
double the natural line-width; a second factor of 2 because we assume 
the spectrum to consist of two absorption lines; and another factor of 
2 because we assume the product of recoilless emission (fs)- and recoilless 
absorption (fa)-probability to be 0.5). Due to this factor of 8 and the 
limited abundance of only 2.2% of Fe5? we derive an effective resonance 
cross-section of Sett ~0.5- 10-26 cm?. In order to measure an absorption 
effect of 2.5%, we therefore need about 1017 Fe5’ atoms per cm?, which 
is equivalent to 5- 1018 iron atoms per cm? (now including all isotopes 
or iron). Using a solution with a 10% concentration of hemoglobin 
(0.1 g Hb/cm$) from which a sample, 0.5 cm thick, was prepared (con- 
taining H20 and Hb) the total amount of 0.05 g Hb/cm? in the absorber, 
corresponded to about 5: 1017 Hb molecules per cm? or 2 : 1018 iron 
atoms per cm?. Comparing this value with the foregoing values we find 
that a hemoglobin absorber with an effective thickness of 0.5 g/cm?, 
with a concentration of 0.1 g Hb/cm?®, and with a natural abundance of 
2.2%, Fe5? produces an effect of roughly 1%. Since we are limited by 
the absolute thickness of an absorber due to strong mass absorption 
of the 14.4 keV Fe5? Mossbauer line (about 60% of the 14.4 keV y-beam 
passing through 0.5 g/cm? of hemoglobin or myoglobin solution is 
absorbed), and because it is desirable sometimes to work with less 
concentrated material and within a reasonable time limit for the measure- 
ment, it is necessary to enrich the sample with Fe®’. With in vitro 
enrichment artificial abundances up to 90% of Fe5? may be obtained, 
(49) whereas with im vivo enrichment artificial abundances may be up 
to about 30% (44). 

Assuming that our absorber contains enough Fe5? to produce an 
effect of 2.594 we can estimate the time of measurement if we know 
the counting rate N’ per sec of the 14.4 keV y-radiation (which depends 
on the activity of the source and on the geometry of the experimental 
set-up), and if we limit the accuracy of the measured effect to 5%. 
Under this condition we need N counts in each of the 400 channels of 
the multichannel analyzer and from Fig. 6 we see that the total error 
of the measured effect is given by 2 VN/0.025 N = 0.05, leading to 
N =2.5-106 counts per channel. With a source activity of 100 mCi, 
a source detector distance of about 12 cm, a detector window of 2.54 cm 
diameter, and an absorber with effective thickness of about 0.5 g/cm? 
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N statistical error: +VN 


25 % effect corresponds to 0.025 N 


N{1 -0.025) BE statist. error about: + VN” 





Vr V 


Fig. 6. Typical single line absorption spectrum 


of Hb solution a counting rate N’ = 5 : 105 of 14.4 keV y-particles per 
sec was measured. In order to collect 109 counts per 400 channels we 
therefore needed approximately 35 min. In cases where we are concerned 
with magnetic hyperfine split spectra, the time required for measure- 
ment easily increases by a factor of 20 because the two absorption 
lines for which this estimate was carried out may spread out to a variety 
of absorption lines with much smaller intensity. 


IV. Experimental Results of Frozen Solutions and Their 
Computer Analysis 


1, Experimental Results 


From the application of the recoilless y-absorption technique to MbCO, 
HbCO, MbOs, HbO2g, Mb, and Hb Mössbauer results have been derived 
which are presented in this section. Fig. 7 shows Mossbauer spectra 
of a frozen solution of MbCO (pH 7.0) at 4.2 °K (50), a typical candidate 
for ferrous low-spin state. Curve (a) coresponds to a measurement with 
zero applied magnetic field. Assuming a Lorentzian line shape (see 
Table 1), a least-squares fit to this spectrum leads to the following 
values for the Méssbauer parameters: /' (line width) =0.328+0.011 
mm/sec, 6 (isomer shift, relative to iron metal) =0.266 +0.010 mm/sec, 
and AEg (quadrupole splitting) = 0.363 + 0.006 mm/sec. Curve (b) 
shows a spectrum taken with a magnetic field of Ho = 47 kOe ap- 
plied perpendicularly to the y-beam. Both spectra have been found 
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Fig. 7. Experimental Mössbauer spectra of frozen solution of MbCO, (a) T =4.2°K 
Ho =0. (b) T=4.2°K, Ho = 47 kOe, Hol y. The additional resonances at —0.9 
mm/sec and +1.33 mm/sec are due to oxygenated myoglobin (MbCOg) which 
remained in the sample through the preparation of MbCO. The broad shoulders 
in (b) are probably due to the impurity effect of MbOg. [Taken from Ref. (50)] 


to be nearly temperature independent within the wide range of 4.2 °K < 
T< 200 °K. 

Since MbO% is a relatively unstable compound, which easily shows 
autoxidation to metmyoglobin, most of the investigations concerning 
oxygenated heme proteins have been made using HbOg of different 
species. Fig. 8 shows Mössbauer spectra of oxyhemoglobin in red cells 
at 4.2 °K with: (a) Ho =0, and with (b) Ho =30 kOe applied perpendic- 
ularly to the y-beam (5). From analysis of curve (a) the following param- 
eters are derived: F= 0.44 mm/sec, ô (relative to iron metal) =0.2 
mm/sec, and AEg=2.2 mm/sec. A4Eq of HbQz is found to be tempera- 
ture dependent, falling to about 1.9 mm/sec at 195 °K. 

Spectra of deoxygenated hemoglobin Hb, a typical candidate for a 
ferrous high-spin ground state, are presented (5, 57) in Fig. 9. Curve 
(a) corresponds to Ho =0 and 4.2°K, (b) to Hp =60 kOe, Ho 1 7, and 
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Fig. 8. Experimental Méssbauer spectra of HbOg in red cells. (a) T = 4.2 °K, Ho = 0, 
(b) T=4.2°K, Ho =30 kOe, Hol 2. [Taken from Ref. (5)]. Absorption lines in 
(b) approximately correspond to transitions a-f for Vzz<0 in Fig. 10 


T =4.2°K, (c) to Ho =60 kOe, Ro 12, and T=16°K, and (d) to 
Ho =60 kOe, Ho L y, and T=35°K. From the analysis of curve (a) 
the following parameters are derived: T'=0.46 mm/sec, ô (relative to 
iron metal) =0.9 mm/sec, and AEg =2.4 mm/sec. To demonstrate the 
strong temperature dependence of electronic properties of Mb or Hb of 
different species, Fig. 9e shows experimental AEg(T) data of Hb (rat, 
pH 10) (52) and Mb (sperm whale) (#2, 53). 


2. Analysis of Magnetic Hyperfine Spectra 


2.1 Diamagnetic Heme Proteins 


Starting with the simplest case, when no electric quadrupole interaction 
is present and when either an externally controllable magnetic field, 
Ho, is applied or, as in the case of metallic iron, when there is an internal 
magnetic field, Hint, we shall be concerned with Zeemann splittings of 
the nuclear levels as indicated in Table 1. The energy level positions Em 
are then dependent upon the effective magnetic field Her =H + Aint, 
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Fig. 9. Experimental Mössbauer spectra of Mb and Hb. (a) Hb at 4.2°K, Ho = 0, 
(b) Hb at 4.2°K, Hp =60 kOe, olL ¥, (c) Hb at 16°K, Hp =60 kOe, Hoi >, 
(d) Hb at 35 °K, Hp = 60 kOe, Ho L >. [Taken from Refs. (5, 57)]. (e) Temperature 
dependent quadrupole splittings for Mb and Hb 


(as detected at the nuclear site of Fe5’) through the relation Ey = — gn 
un Her I. Here gy is the nuclear g factor, which is different in the 
ground state and the excited state: gy =0.1805 and gy = — 0.103 for 
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Fe5?, The quantity un =ehl2 Mc=3.15-10-12 eV/Oe is the nuclear 
magneton, and Jz=+1/2 and I; = +3|2, +1/2 stand for the mag- 
netic components of the nuclear spin I in the ground and excited 
state of Fe5’, respectively. In diamagnetic compounds we do not expect 
electronic contributions to the effective magnetic field at the Fe? 
nucleus, i.e. Hett is equal to Ho. In the two diamagnetic compounds 
discussed here, namely MbCO and HbOz, electric quadrupole interaction 
must be considered in addition to nuclear Zeemann splitting. From 
Fig. 8b for HbO% it is obvious that the applied field Ho causes only a 
small magnetic hyperfine splitting of the nuclear levels, AEm (+3/2) 
and AEm (+1/2), compared to the quadrupole splitting AEg of Fig. 8a. 
Under the approximations AEg> AEm (+3/2), with AEm (+3/2) = 
— gyn un Ho ( —3|2 — 3/2) =3 gn un Ho, Ho =30 kOe, and y (asymmetry 
parameter, see Table I) = 0, the nuclear levels [I =1/2, fz7=-+1/2> 
for the groundstate, and |3/2, +1/2> and |3/2, +3/2> for the excited 
state will be split into |1/2, +1/2> and |1/2, —1/2> by 0.566 mm/sec, 
into |3/2, +3/2> and |3/2, —3/2> by 0.969 mm/sec, and into |3/2, 
+1/2> and |3/2, —1/2> by 0.323 mm/sec, respectively. Acocrding to 
the sign convention for Vzz from Table 1c we then get a different energy 
sequence of the eigenfunctions |3/2, Iz> for Vzz>0 than for Vzz< 0. 
Further, taking into account that we are concerned in the Fe5? Möss- 
bauer transitions with magnetic dipole transitions, for which we have to 
consider the selection rules A/z=0, +1, transition or absorption lines 
are obtained as indicated in Fig. 10. Comparing this result with the 
experimental magnetic hyperfine pattern for HbOs in Fig. 8b we derive 
that the EFG component Vzz of Fe5? in HbO% is negative. This informa- 
tion is important for the further interpretation of the results. The 
reason that the experimental spectrum of Fig. 8b shows unequal inten- 
sities for the transitions e and f, or for d and a, compared to the ideal 
case of Fig. 10, is due to the fact that 7 in HbO; is slightly larger than 
zero, and hence the real nuclear levels will be linear combinations of 
basis functions |I, Iz>, 

For the analysis of the magnetic hyperfine spectrum of HbCO (Fig. 
7b) the approximation AE g > AE (+3/2) is no longer valid. In this 
more general case the Fe®? nuclear eigenstates were derived with the aid 
of the hamiltonian. 


A = pO BI — 10 +1) + mE) — gn uw Aol (D 


All parameters of Eq. (1) have been defined in Table 1 as well as 
the nuclear magnetic moment operators Í, and /, which may be sub- 


stituted by f,=1/2 (ËI +Í and f,=1/2 (I; —L.), respectively. 
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Fig. 10. Determination of sign Vzz under conditions 7 =0 and AEg>AEy 


Together with the definition that the quantization axis z is parallel to 
the direction of Ao Eq. (1) becomes 


eO V > 
f = op (SHEE +) +U(R+E)] Ten. 0) 


Only for the special case that Vzz || A, and n=0, the resulting 
matrix <I, /z|A| I, I;> is diagonal and the states |I, Iz> are eigen- 
functions of A. Assuming the EFG to be axial in its main axis system, 
sv ge [7 = K =0), but having nonparallel Vzz and Apo, 

z z 
the transformation 7’+2Z=Hp/|Ho|will cause a non-zero asymmetry 
| Vex — Vyy| 

A I 
“low symmetry’’-case the operators Í H and {2 produce off-diagonal 
matrix elements <I, Iz+2|H|J,Iz>, and the eigenfunctions to A 
contain admixtures with AI; = +2. Thus, new eigenfunctions |J, «> 
can be expanded in terms of the basis functions |Z, Iz>:|J, «> =>Ca,1, 
[I ‚Iz>. Iz 

In order to find the nuclear eigenstates Z7 and the expansion coef- 
ficients, C; r, the two matrixes < Ig, Iz |H|Ig,Iz> for the ground 
state and </e, Iz |i | Ze, [z> for the excited state of the Fe5? nucleus 
have to be diagonalized. To calculate transition (absorption) probabilities 


parameter 7 = #0 in the new axis system x, y, z. In this 


Ves = | < Íe, a| P | Ie, x > |2 = | 2 Ce’, Iz Cx, Iz < Tle, Iz | 
|? [Ie Iz > |? L (3) 
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during the absorption of an unpolarized single y-line with energy Ey 
and y-propagation direction Ë we have to bear in mind that |< |P? | >|? 
depends on the angle 0 between Fo and &. For Fe5’, as already mentioned, 
the radiation is of a pure magnetic dipole type with allowed transitions 
specified by Iz — Iz =0' +1. The 6-dependence of |< |P | > |? is then 
given by 




















3 3 1 1 2 3 3 1 1 2 
<+ +719 12. +3 > =\<2.-71 912. a> ~ 
(1 + cos? 6}/8 
3 1 1 1 2 3 1 1 1 2 ; 
la ale aer u 25 ae ~ sin? 
<p tg lp +> =|<2 ZIP 15-3 >] ~ sincere 
3 1 1 1 2 3 1 1 1 2 
eaa 4 -|< 4. +3 1912. -3> — (1+ 
cos? 6)/24 . (4) 


So far we have considered only a line absorption spectrum with in 
general eight lines (Fig. 11) having energies EY, and E?, («=1,2 and 


Je: 5 | a 
J: 7 beg 


Fig. 11, Generally possible nuclear transitions E2+ EN, for a = 1,2 and e = 1,2,3,4 
in the case of Fe?” 


a’ = 1,2,3,4), and intensities Vg’, as given by Eqs. (3) and (4). In order 
to get a Mossbauer spectrum /(£) a Lorentzian line shape characteristic 
of the energy widths of emitting (Mem) and absorbing (Dabs) nuclei has 
to be folded into the line absorption spectrum: 


| > Ca’, Iz Ca, Iz < Te, IZ |P] Ig, I;> |° 


iger > Iz > (5) 


Ter? E, Erë + (Bš 


+ 
a sa 





with T = Fem + Taps, and Era=Er Pr: 
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A corresponding relationship for /(v) results from replacing the 
energy unit “eV” by the Doppler velocity units “mm/sec” : 1 mm/sec «> 
4.8 - 10-8 eV. Eq. (5) describes an absorption spectrum of a single crystal- 
line absorber, since only a single orientation of the EFG tensor Vy, 
relative to Hy has been regarded so far. For the simulation of spectra as 
derived from frozen solutions or polycrystalline samples, all orientations 
of Vya relative to Ho must be taken into account by an appropriate 
average over possible single-crystal spectra of Eq. (5). 

The spectrum shown in Fig. 7b was analyzed by calculating a Möss- 
bauer spectrum for polycrystalline MbCO along the lines decribed above. 
Fig. 12 shows calculated spectra obtained by using the input parameters 


Velocity in [mm/sec] 


Fig. 12. Theoretical Méssbauer spectra for frozen solution of MbCO, calculated 
along the lines described in IV.2.1. (a) Vzz> 0, n=0. (b) Vzz> 0, 2 =0.4. (c) 
Vzz<0, n = 0. (d) Vzz<0, 7 =0.4. All spectra have in common: J'= 0.328 mm/sec, 
ô = 0.266 mm/sec, A Eg = 0.363 mm/sec, Ho = 47 kOe, and Ho | z. [Taken from 
Ref. (50)] 
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T =0.328 mm/sec, ó =0.266 mm/sec, Vzz = with Á Eg = 


eQ ( +2) a 

0.363 mm/sec, Ho =37 kOe, and Ho | y. Curves (a) and (b) are for 
Vzz>0 and correspond to n=0 and n=0.4, respectively; curves (c) 
and (d) are for the same n-values but with Vzz<0. A comparison of the 
experimental spectrum of Fig. 7b with the calculated spectra leads to the 
conclusion that Vzz at the Fe5? nucleus in HbCO is positive. Again this 
information will be important for the further interpretation of the results 
as shown in Section VII. 


2.2. Paramagnetic Heme Proteins 


In computing Méssbauer spectra of paramagnetic compounds (54) it 
should be taken into account that due to spin-lattice or spin-spin relax- 
ation the paramagnetic centers may fluctuate between their eigenstates. 
In the cases of (i) fast, and (ii) slow relaxation the calculation of a Möss- 
bauer spectrum remains relatively simple [Kündig (54)]. In the first case 
(i) the spectrum is determined by the thermal average of the hyperfine 
fields, in the second case (ii) the spectrum results from a superposition 
of spectra, which correspond to the hyperfine fields of the electronic 
states, weighted according to thermal population. The intermediate case 
becomes much more complicated [Blume (54)]. As an example of mag- 
netic hyperfine interaction in paramagnetic heme compounds we describe 
case (i) and refer to the literature for other cases. 

The magnetic hyperfine interaction in paramagnetic compounds con- 
sists of nuclear and electronic components. The strength of this inter- 
action in the case of Fe57 corresponds to a nuclear field of order 10 Oe 
acting on the electrons, or to an effective electron field of a few 100 kOe, 
acting on the nucleus. Thus, in general, the magnetic hyperfine inter- 
action hamiltonian Hy: operates on basis states which have both elec- 
tronic and nuclear components. Because the contribution of the hyper- 
fine interaction to the electronic energy in ferrous iron is usually small 
compared to typical electronic zero field splittings, the electronic and 
nuclear parts of the problem can be solved separately. Only in the case 
of ferric complexes, which have half integral spin, does the Kramers 
degeneracy prevent the separation of the two problems. A magnetic 
field of Ho> 10 Oe (Hy =10 kOe corresponds to some cm-1), exter- 
nally applied to the spin degenerated system, however, will lift the Kra- 
mers degeneracy and again makes the problem separable. Under this con- 
dition, in Ap, the spin operators $,, etc. are replaced by spin expectation 
values <S,>, etc.; then we are only concerned with a nuclear spin 
hamiltonian. 
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The Fe5’? nucleus now experiences a net expectation value of the elec- 
tronic spin which is a thermal average taken over all spin states being 
populated by temperature: 


_A a) 
exp ( 


(6) 





<Si>av = 2 <S>N exp ( 





In Eq. (6) the index i designates the components x, y, z of the electron 
spin, and N is the sum of all the populated states; AE y is the energy of 
the Nth electronic spin state relative to the lowest one. The nuclear 
Hamiltonian, containing the magnetic hyperfine interaction between 
nucleus and paramagnetic valence shell, the nuclear Zeemann term 
(—gnunfAol), and the quadrupole interaction Ha, is given by 


H = < Š See AT — gnun Hol + Ho 
SSA 7) 
= — NUN [Ho - =") r +A, \ 
ENUN 
with A representing the so-called magnetic hyperfine tensor. The 
quantity — <> is commonly called the internal magnetic field Hin 
N YN 


at the nucleus, produced by the electronic system. At this point we should 
bear in mind that this quantity entering (7) usually applies well for ferric 
low-spin (S = 1/2) and ferrous high-spin (S =2) compounds, since their 
paramagnetic centers are characterized by fast spin-lattice relaxation 
(case i). In some specific ferrous high-spin compounds (Zimmermann (54) ) 
and in most of the ferric high-spin (S =5/2) compounds, however, we 
are concerned with slow spin-lattice relaxation (case (H)). One ferric 
high-spin compound for which Eq. (7) applies in order to calculate its 
magnetic hyperfine split Mössbauer spectrum is hemin, because of fast 
spin-spin relaxation. Thus, the nucleus experiences an effective magnetic 
field Hess, which is the vector sum Her =Ho+Hint. Observe that 
Hint in general is quite a complicated quantity and depends on the 
A-tensor, through <Š >ay on the temperature because of the thermal 
averaging procedure, on the construction of the energetically low lying 
molecular orbitals and their mutual spin-orbit coupling, which itself may 
be dependent on o. To specify a theoretical Méssbauer spectrum ac- 
cording to Eq. (5) completely, the parameters <Sz>ay, <Sy>av, 
<S,>ay, Az, Ay, Az, Vzz, n, Ho, ô, T, the nine Euler angles defining 
the orientation of the A-tensor, of the EFG tensor Vpp, and of the 
vector <Š >ay with respect to H oare needed. The task is then to measure 
Mössbauer spectra at various temperatures in applied magnetic fields 
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of different strengths and directions and to find a consistent set of param- 
eters describing all spectra. Fortunately some of these parameters are 
known from theoretical estimates or from other measurements. The 
expectation values <S;>ay may be approximated in terms of zero- 
field splitting D and rhombicity E, which sometimes are available from 
far-infrared and susceptibility measurements, and in terms of the g- 
tensor, which for compounds containing ferric iron often is known from 
ESR-measurements. Endor measurements can provide the absolute 
values for the components of the A-tensor of ferric iron and its orien- 
tation relative to the g-tensor. High temperature Mossbauer data give 
the line width J’, the isomer shift ó, and the quadrupole splitting A Eg, 
and if single crystals are available, Vzz and », and the orientation of 
V pp With regard to a fixed framework can be derived. Computer programs 
available to handle the complicated problem of simulating Mössbauer 
spectra of paramagnetic heme proteins have been written separately by 
G. Lang (55) and E. Münck (56). 

Generally the spin lattice relaxation between spin states of high- 
spin ferrous iron is known to be fast. If the relaxation time is small 


: : : h 
compared with the precession time t of the Fe®? nucleus (z = -E where 


AE is a typical energy splitting of about 10-7 eV of the nuclear levels), 
the nucleus then sees a net expectation value of the electronic spin, 
<Š >av, which is zero. This is the case for deoxymyoglobin, deoxy- 
hemoglobin, reduced cytochromes, and reduced peroxidases etc., which 
show with zero applied field (Hp =0) only a quadrupole doublet (see 
Fig. 9a). With applied field (Ho # 0) the spins align more or less along 
nm. contribution to the effective 
EN HN 

magnetic field Herr at the Fe5? nucleus. 

Fig. 9b shows the experimental Mössbauer spectrum of Hb at 4.2 
K with Ho =30 kOe and # | 7. So far, it has not been possible to com- 
pute a theoretical Mossbauer spectrum, which fits the experimental pat- 
tern of Fig. 9b, mainly due to the fact that most of the parameters re- 
quired for such a fit have not been available. 

Magnetic measurements (Ro Z 0) of high-spin ferrous systems at 
elevated temperatures sometimes have the advantage that the thermal 
average <S>ay is zero. Figures 9c and d show spectra of Hb (57) meas- 
ured at 16°K and 35 °K, respectively, with Ho =60 kOe perpendicular 
to the transmitting y-rays. The spectra are similar to that found for 
diamagnetic HbO% (Fig. 8b), indicating that <Š >ay is approximately 
zero. Only the pronounced isomer shift of ó = + 0.9 mm/sec (relative to 
iron metal) indicates that we are concerned with a high-spin ferrous com- 
pound, for which an isomer shift (relative to metallic iron) of the order of 


Ho, thus producing a non-zero 
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+1 mm/sec is characteristic. From the line-triplet on the low-energy 
side (Fig. 9c) it appears that the sign of Vz, for Hb is positive (cf. Fig. 10). 
This conclusion agrees with the findings from single crystal measure- 
ments with Mb (42), which are described in the following section (V). 

From inspection of the spectrum of reduced cytochrome P450 (Fig. 
13), measured by Münck et al. (57) at 200 °K with Ho =45 kOe parallel 
to the transmitting y-rays, it was concluded by a similar argument (line 
triplet on the high-energy side) that Vz, is negative (57). A most recent 
computational analysis of 4.2 °K- and 220 °K-spectra of reduced cyto- 
chrome P450 by Münck (private communication, May 1974), however 
indicates that the “triplet”’ argumentation concerning sign Vz, is very 
deceptive, at least for this compound; the overall theoretically derived 
intensities now seem to favor Vz, > 0. 
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Fig. 13. Experimental Mössbauer spectra of frozen solution of reduced cytochrome 
P 450 at 200°K, Ho = 45 kOe, ol| ?. [Taken from Ref. (57)] 


In order to complete the discussion of magnetic hyperfine interactions 
in paramagnetic heme proteins as detected by Mössbauer spectroscopy, 
reference should be made to the work of Champion et al. (58). There the 
influence of the halides F, Cl, and I on ferric chloroperoxidase was 
studied. The results presented in this work indicate that halide anions, 
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which are substrates for chloroperoxidase halogenation reactions, can 


form complexes with the enzyme and cause significant changes in the 
heme environment and the spin state of the iron. Fig. 14a-c show signif- 
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Fig. 14. Experimental and theoretical Méssbauer spectra of ferric chloroperoxidase 
and its halide complexes. (a) Native chloroperoxidase, pH 3, T=4.2 °K, Ho = 
1.3 kOe, Ho | ¥. An identical spectrum was found for the chlorine-complex measur- 
ed under the same conditions. (b) Iodine complex, pH 3, T=1.5°K, Ho = 1,3 kOe, 
Ho 1 2. (c) Fluoride complex, pH 3, T =4.2°K, Ho =1.3 kOe, Ho | 2. [Taken 
from Ref. (58)] 
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icant changes in the Méssbauer spectra of chloroperoxidase upon complex 
formation with different halides; this feature is remarkable especially 
under the aspect that the halides are probably not directly coordinated 
to the iron. The change in the magnetic hyperfine patterns in Fig. 14a-c, 
which have been simulated by Champion et al. (58) along the lines de- 
scribed above (solid lines in Fig. 14a-c), might suggest the “presence of 
transferred hyperfine interactions between the halides and the iron” 
(58). 


V. Experimental Results of Myoglobin Single Crystals and 
Their Analysis 


1. Experimental Methods and Results 


Myoglobin (instead of hemoglobin) was chosen for single crystal Möss- 
bauer measurements (42) because myoglobin contains only two hemes 
within one unit cell, one being transformed into the other by a subsequent 
rotational (180° around the crystallographic b axis) and translational 
transformation. This condition together with the geometric arrangement 
shown in Fig. 15, with the y-ray direction always perpendicular to the 
b-axis, has the advantage that identical Mossbauer spectra are measured 
for the two heme groups (59). 

Single crystals were grown from a solution of sperm whale myoglobin. 
The iron of the myoglobin was enriched to 80% with Fe57. A met myo- 
globin (met Mb) crystal with a diameter of about 0.5 mm was fixed in a 
Debye capillary with the aid of two smaller Debye capillaries. The 
reduction of a met Mb crystal was accomplished by washing the crystal 
in the capillary with a solution of 4 M (NH4)2504-+ NasSsO4 in water. 
The capillary was then closed by wax and mounted on a goniometer 
head. The orientation of the deoxy-myoglobin (Mb) single crystal was 
determined on a precession camera (50), with the sample kept at 77 °K 
by washing the capillary with a fine spout of liquid nitrogen. After the 
orientation was determined the whole goniometer head was placed in 
liquid nitrogen until it was used for measurements. During the transfer 
of the goniometer head to the cryostat the capillary was carefully washed 
again with liquid nitrogen. The crystal remained at 77 °K throughout the 
Mossbauer experiments. After that the goniometer head was back — again 
the sample was maintained at 77 °K — to the precession camera to re- 
check the orientation of the crystal. 

During the freezing procedure the crystal water expands, therefore 
the crystal is somewhat damaged. Since X-ray pictures on the precession 
camera were also obtained after freezing, the crystal remained satis- 
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Fig. 15. (a) Schematic of the orientation relations in the Mb single crystal experi- 
ments. (b) Projection of the heme on the a’b plane. Heme nitrogens are denoted 
by nos. 1 to 4. @ is the angle of rotation around the crystallographic b axis 


factory for the purpose of the experiment; the misorientation of the 
whole single crystal remained smaller than 5°. 

The single crystal orientation in the experimental arrangement is 
shown in Fig. 15. The crystallographic axes are represented by a, b, and 
c. To simplify the description of the experiment a’ is introduced, so that 
the axes a’, b, and c are mutually orthogonal. The propagation direction 
of the y-ray, 7, remained perpendicular to b throughout all experiments. 
The crystal was rotated about the b axis by the angle gç. The angle 
g =45° corresponds to the situation 7 || a’, and g =135° to Jljc. The 
heme plane is also shown in Fig. 15. When projected on the a'b plane, it 
makes an angle of 22° to b (Fig. 15b). 
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A typical spectrum of an Mb single crystal is shown in Fig. 16. The 
small contamination at velocities —0.913 mm/sec and + 1.357 mm/sec 
is due to MbOg. Table 2 summarizes line widths /', quadrupole splittings 
AEg, and isomer shifts ó (relative to metallic iron) for Mb and MbO> 
as derived from the single crystal spectra, for Mb in frozen solution and 
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Fig. 16. Experimental Mössbauer spectrum of a Mb single crystal with orientation 
b17 and g = 156° (g = 45° corresponds to a’ || 7 and g =135° to c || 7). Side peaks 
originate from MbOg. [Taken from Ref. (42)] 


Table 2. Experimental®) Mössbauer line widths J’, quadrupole splittings AEg 
and isomer shifts 6 for oxygenated and deoxygenated myoglobin 


Material I’ (mm sec-1) AEg (mm sec-1) ôb) (mm sec-1) 
mbO2°) 0.402 + 0.017 2.269 + 0.017 +0.222 + 0.017 

Deoxy mb$) 0.402 + 0.017 2.286 + 0.017 +0.900 + 0.017 
mbOe4) 0.36 + 0.03 2.27 + 0.03 +0.22 + 0.03 

Deoxy mb®) 0.37 + 0.03 2.24 + 0.03 +0.86 + 0.03 


a) All experiments at 77 °K. 

b) Isomer shift relative to «-Fe at 300 °K. 
e) Single crystal A, see text and Fig. 16. 
d) Polycrystalline powder. 

e) Frozen solution. 
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for MbOg in polycrystalline powder. Table 3 contains the intensity ratios 
R (see Eq. 8) for the eight different orientations of Mb single crystal A 
and for two different orientations of Mb single crystal B. In Fig. 17 the 
values from Table 3 are plotted in a R-p framework (closed circles). 


Table 3. Experimental ratios of quadrupole split Mössbauer 
lines of deoxy mb single crystals under various orientations 


A B 
g) RP) ga.) RP) 
45 0.823 + 0.140 5047 1.065 + 0,192 
66 1.017 + 0.163 
88.5 1.121 + 0.179 
111 1.100 + 0.176 
135 0.875 + 0.140 135 +7 0.783 +- 0.140 
156 0.705 + 0.119 
178.5 0.658 + 0.118 
200.5 0.700 + 0.126 


a) Angles given in degrees. 
b) R is the ratio of intensity at position 42.043 mm/sec. 
and of intensity at position —0.243 mm/sec: R = 
I (+2.043) | r (-0.243) - 
Deoxy mb single crystal (A) was calibrated on precession 
camera (see text). The R (g) — values are plotted in Fig. 17. 
c) Deoxy mb single crystal (B) was calibrated by visual 
adjustment only; therefore we give an error for g. 


3.0 
2.8 


2.6 





o 


0.6 so re 
45° 90° 135° 180° 225° 
} | $=- 


{a c) 


Fig. 17. Relative line intensity R (see Table III and text) versus angle g@. [Taken 
from Ref, (42)] 
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2. Analysis of Myoglobin Single Crystal Mössbauer Spectra 


The relative line intensity R of a quadrupole split spectrum is determined 
by the angle 0 between 7 and the principal axis of the electric field gra- 
dient z=Vzz. From Eqs. (3) and (4) and assuming the assymmetry 
parameter 7 to be zero, the intensity ratio R of the two absorption lines 
of the quadrupole split Mb single crystal spectrum is derived to be: 


2 
R = I3/2 1 + COS 0 (8) 
11/2 


2 * 
— + sin? 0 
3 


Because 7 has little influence (67) on R, it is assumed to be zero for 
a first estimate. The angle between z and the 5 axis is designated by «. 
The three angles under consideration g, 6, and « are related to each 
other as shown in Eq. (9): 


cos 0 = sin < : sing. (9) 


From the R-values in Table 3 and Fig. 17 it is obvious that Vzz is 
positive. By rotating the z-axis around 0 there will be 0 =90° for one 
specific g and therefore R will be either 0.6 or 1.667, depending on wheth- 
er Vzz is positive or negative. As the R-values in Fig. 17 never exceed a 
value of about 1.2 it is concluded that Vzz is positive. 

If it is assumed that z is perpendicular to the heme plane, then, under 
the condition 7|]a’ (p =45°), R is calculated to be about 2.3 (indicated 
by a square in Fig. 17). However, if it is assumed that z is in the heme 
plane parallel to a'b, then R is about 0.73 (as indicated by a circle in 
Fig. 17). In both cases R = 0.6 should be found for + ||c (g =135°, also 
indicated by square and circle in Fig. 17), because the z's become per- 
pendicular to 7. Assuming z is in the heme plane but perpendicular to a'b, 
R is calculated to be R =0.6 for 7||a’ and R =3 for #||c (indicated by 
triangles in Fig. 17). The experiment (closed circles in Fig. 17) obviously 
shows quite a different behavior. 

The minimum and the maximum in R occur at g = 180° and g = 90°, 
respectively. At these positions 7 forms an angle of 45° to both the a’ 
and c axes. The minimum R ~0.6 at ç = 180° indicates [Eq. (8)] that 
0 =90° or 7 1z. At p=90°, R is about 1.13 and the smallest angle, 
6 =50°, between 7 and z is obtained. Under this condition, ¥, z, and 0 
are all in one plane and therefore 0 -+ |< =90°. It follows that the angle 
between z and b is « =40° (+8°), which is considerably larger than the 
angle of 22° indicated in Fig. 15. Calculation of R using Eqs. (8) and (9) 
and taking «=40°, yields the solid line drawn in Fig. 17, which is in 
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good agreement with the experimental points. In Ref. (42) the angle 
a =40° +8° was derived assuming 7 to be zero. A most recent com- 
putational analysis of the experimental data of Table 3 by Y. Maeda 
et al. (62) by taking use of Eq. (5) with Hg setting zero, however, indicates 
that x =40° +8° is consistent also with an n-parameter in the range 
0< 7< 0.3. 

Now it is possible to locate z in the single crystal. z forms an angle of 
40° +8° with the b axis, and when projected on the a'c plane forms an 
angle of about 45° to both axes. Then, in regard to the heme plane two 
alternatives of the z or Vzz orientation are consistent with the measure- 
ments: Vzz forms an angle of about 45° with the heme plane or it is 
oriented rather close to the heme plane and is directed toward one of the 
nitrogen atoms of the heme ring. The latter seems more likely considering 
the structure around the Fe nucleus. In Section VII we shall discuss a 
further confirmation of this assumption as derived from molecular orbital 
calculations. 


VI. Interpretation of Experimental Quadrupole Splittings 
in Crystal Field Approximation 


The interpretation of experimental quadrupole splittings, 4 Eq, of 
Hb was first based on electric field gradient contributions from Fe 34 
single-electron orbitals weighted by an appropriate Boltzmann factor 
(3, 63). However, it turned out that a single electron in a Fe 34, orbital 
would produce a 4 Eq value which is nearly twice the experimental 
value measured for Hb at 4.2°K. Therefore, it was investigated which 
of the various electron terms within the 348 configuration may lie close 
enough in energy to the high-spin ground state in Hb, in order to mix 
negative EFG contributions into the positive ground state EFG via 
spin-orbit coupling between ground and excited electron states (52). 
Otsuka (64) and Eicher et al. (52) have shown independently that in a 
ligand field of tetragonal symmetry (taking the heme normal as symmetry 
axis) the energy levels °Bo, 5E, 3E, and 1A ı may fulfill this requirement. 
This situation is illustrated by the results in Fig. 18, which have been 
computed by taking into account the proper tetragonal point symmetry 
(Cay) of the heme iron and the Coulomb repulsion of Fe 3d electrons. 
From Fig. 18 it is obvious how the many-electron term energies depend 
on the single electron energies e; (t =1, 2, 3), which are defined by e1 = E 
(dzz, dyz) —E (dzy), 89 = E (d3z2--2) — E (dzy), and e3=E (dz2_y2) — 
E (dgy). For energies e3 > 20000 cm —! the spin triplet 3E is found to be 
ground state, which is the case for Fe in iron phthalocyanine (65); for 
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Fig. 18. Crystal field results for energy positions of low lying many-electron states, 
relative to the lowest 3E term within the 3d configuration, plotted as a function 
of single-electron energy ea. The solid and dashed lines represent two different 
calculations using ei = 600 cm-1, e3 = 22000 cm-1, and &ı = 600 cm-1, e3 = 16000 
cm-1, respectively. [Taken from Ref. (52)] 


strong ligand fields characterized by £2 > 10000 cm-1, e3 315000 cm-t 
we expect the spin singlet as ground state, which is the case for the 
diamagnetic compounds MbCO and HbCO. In the weak field approxima- 
tion (small &) obviously a spin quintet will be ground state. Using the 
estimate (66—69) for s; (Fig. 19) from various experiments with heme 


; ybC- Y) 
H = estimated from magnetic 
r ee u ae susceptibility on me 
Fo 1 (374 ferric Myoglobi 
T 2y3 a=] j ee 
Pa 1 Absorption spectra 
3d u 1,000 cm”'=- and CD spectra on 
KNN ferrous Cytochrome C 
xU S yz 
N 1 magnetic susceptibility 
` — 
x 2000 a on ferrous Myoglobin 69 
Iron ion in 


free iron ion ferrous Myoglobin 


Fig. 19. Differences of single-electron 3d energies in reduced heme proteins. The 
arrows indicate the procedures of the estimates of these energies. [Taken from Ref. 


(66) ] 
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compounds it seems reasonable to take into account spin-orbit inter- 
action between the terms 5Bo, 5E, 3E, and 14, in order to calculate 
A Ea(T) for Mb or Hb. The solution of the corresponding secular problem 
for the low lying levels 5By, 5E, 3E, and 1A, is expressed by the energies 
ej; and the 22 spin eigenvectors |ex> are certain linear combinations 
of the base vectors each contributing to 4 Eg by typical high-spin, 
intermediate-spin, and low-spin values. In order to evaluate the tempera- 
ture dependence of A Eq it is necessary to average the EFG components 
Vzz and 7° of the individual substates |e,> according to Boltzmann 
statistics. A detailed description of the evaluation of A Eg(T) will be 
given in the following section. 

The basic results from Eicher’s and Trautwein’s calculations (52) of 
A Ea(T) for Mb and Hb follow from inspection of Figs. 20 and 21. 
Comparing experimental and theoretical A Eg(T) data the following 
energy sequence of many electron terms relative to the ground state °B3, 
is derived: E (5B2)=0, E (3E) ~250 cm-1, E (1Aı) ~2000 cm-1, 
and E (5E) ~2500 cm-1. The energies may be altered due to rhombic 
perturbations, which split the doublets 3E and 5E by about 120 cm-t 
and 200 cm}, respectively: E (5Bo)=0, E (3E) ~250+60 cm-t, 
E (141) ~2000 cm-1, and E (5E) ~2000+ 100 cm-t. This term scheme 
for Mb or Hb with the relatively low lying spin triplet 3E, however, 
must be questioned, because: (i) the recent single crystal Mössbauer 
investigation of Mb (42) indicates that Vzz is more likely to be parallel 
to the heme plane than perpendicular, in contrast to the former 
assumption (52), (ii) temperature-dependent susceptibilities y (T) 
derived with this energy term scheme for Hb are inconsistent with 
experiment (70), and (iii) taking account of direct ligand influence 
on the heme iron only through an isotropic covalency factor (an assump- 
tion which is typical for crystal field approaches) seems to be a too 
rigorous assumption. 

Crystal field calculations by Huynh et al. (70) on the basis of Eicher’s 
et al. (52) approach, but taking Vzz > 0 and Vy; as being parallel to the 
heme plane, led to a new energy term scheme for Hb which is consistent 
with both the temperature dependence of the quadrupole splitting in the 
range of 77°K < T <200°K and the magnetic susceptibility, x (T). 
The energies corresponding to Huynh’s et al. A Eg(T) and y (T) fit curves 
are: 


E (5E) =300 + 346 cm}, E (5B) =0 cm-1, E (143) =287 cm-1, 
E (3£)=1026+207 cm-t, and A=76 cm-1 (A: spin-orbit coupling 
constant). 


To study the reliability of this result further it remains to discuss 
the direct ligand influence on the heme iron. This task leads us to the 
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Fig. 20. Temperature dependent experimental quadrupole splittings of Hb (rat). 
The theoretical Á Eg(T) curves result from crystal field calculations (see Section 
VI) taking (a) £1 =228 cm! and A (spin-orbit coupling constant) =82 cm]; 
(b) e1 =268 + 139 cm-1, A=82 cm-1; (c) e1 = 1000 +-980 cm-1, A=82 cm-l; 
(d) si = 234 cm-1, A=77 cm-1. [Taken from Ref. (52)] 


estimate of Méssbauer parameters and susceptibility data on the basis 
of molecular orbital (MO) calculations, which include the three-dimen- 
sional structure of heme compounds beyond the heme iron and its imme- 
diate neighbors and which moreover, include all relevant valence elec- 
trons. 
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Fig. 21. Crystal field results for the temperature dependent quadrupole splitting 
in Hb, involving spin-orbit interaction between energetically low lying levels 
5Bo, BE, 141, and 5E. Their energies relative to the groundstate ®Bg are (for E 
(3E), E (141), and E (5E) in cm-1): 


(a) 269 2085 695 (c) 256 1883 1737 
(a) 338 1390 695 (d 247 1883 2432 
(b) 283 1730 1042 


VII. Interpretation of Experimental Mössbauer and 
Susceptibility Data by Molecular Orbital Theory 
1. General Considerations 


Since heme complexes involve to varying extends covalent binding 
between iron and ligands and charge rearrangements within the ligand 
structure, a theoretical description of the electronic behavior of such 
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complexes should not be limited to the valence electrons of iron only. 
Rigorous all-valence electron treatments of systems as large as hemin, 
however, are prohibitive in cost, even at the self-consistent field (SCF) 
level of approximation. It is even extremely costly to pursue semi- 
empirical antisymmetrized molecular orbital theory or approximations 
based on neglect of differential overlap, such as the CNDO/2 method (77). 
The key simplification appearing to permit calculation at reasonable 
cost is to avoid the construction of an explicit antisymmetrized many- 
electron wave function. Methods in this category include those of the 
Hiickel type (72, 73) and the multiple scattering x. method of Johnson 
et al. (74). 


2. Hiickel and Configuration Interaction Calculations 


The following theoretical studies are based on an iterative extended 
Hiickel method described in detail previously (72). The Hiickel one- 
electron Hamiltonian H was assumed to have diagonal elements Hag = 
— qa —gaÄ xa, where aq is the “a” parameter (ionization potential) for 
atomic orbital a (of real form) in a neutral atom, ga is the net charge 
of the atom at which orbital a is centered, and Aa, is an additional 
parameter controlling the variation of Haa with net atomic charge. 
The off-diagonal elements Hg, were determined from Haa and Hp, by 
the Cusachs approximation (75) Hay=Say (1-3 |Sav|) (Haa + Ho). 
Here Sa» is the overlap integral of atomic orbitals a and b, calculated 
assuming Slater-type orbitals (STO’s). In evaluating Aa», all aperations 
are performed in a coordinate system passing through the orbital 
centers so as to retain full invariance with respect to rotation of the 
molecule in the coordinate system. 

Iterations were carried out until the occupied-orbital solutions of 
the matrix equation (H —E,S)e, =0 were consistent with the net atomic 
charges used in calculating H. Here E; is the orbital energy of the 
molecular orbital whose atomic-orbital coefficients are represented 
by the vector e;. The overlap matrix S is built from the Sap as already 
defined. The iterative process must be sophisticated enough to over- 
come the strong nonlinearity of the theory; a combination of damping 
and second-order extrapolations has been used and thereby convergence 
to +0.01 electronic charge units in atomic charge in five to ten iterations 
was achieved. The net atomic charges are obtained by apportioning 
exactly-computed overlap charges among the atoms involved. 

The calculations include all valence orbitals, z.e. 1s for H, 2s and 26 
for C, N, and O; and 3d, 4s, and 44 for Fe. All orbital parameters have 
been taken from literature and were adjusted to a consistent parameter 
set in a variety of MO studies (72, 76—80). 


138 


Méssbauer Spectroscopy on Heme Proteins 


The Hiickel model is inadequate for distinguishing spin eigenstates 
and for the calculation of energy differences between them, but Mb and 
Hb exist in “high-spin” states necessitating the use of singly-occupied 
molecular orbitals (MO’s). For such problems the Hiickel calculations 
are supplemented by proceeding to a limited configuration interaction 
(CI), simplifying the computations by using the Htickel SCF MO’s to 
construct configurations. As even a limited CI study is expensive for 
large systems, and since the present needs are to characterize the MO 
structure and to compare different spin states, attention is confined 
to calculation of energies relative to that of a “reference configuration’. 
If the reference configuration has energy Eo, energies relative thereto 
can be obtained as eigenvalues of H — EoS, where H and S now refer 
to many electron space-spin states of appropriate symmetries. As Eo 
can be written in terms of the integrals occuring in the elements of H, 
the elements of H — ES will depend only upon the orbitals affected by 
passage from the reference configuration to the configurations actually 
used in H and S. A more detailed description of this CI procedure was 
given earlier by Trautwein and Harris (76), and was applied to study 
Mössbauer parameters of several iron-containing compounds. 


3. Heme Models for Molecular Orbital Calculations 


The model which represents the heme group in MbCO, HbCO, MbO», 
HbOg, Mb, and Hb for the MO calculations is shown in Fig. 22a. This 
planar representation of the heme (73) is based on Koenig’s (87) X-ray 
structure analysis of hemin. Fig. 22b indicates the steric arrangement 
in a plane perpendicular to the heme. The proximal histidine group 
was represented by a water molecule. To study the ligand influence 
upon the electronic structure of iron in MbCO, HbCO various assumptions 
were made as to the Fe—C bond length, the Fe—C—O bond angle £, the 
position of the Fe atom relative to the heme, and the occupancy of one 
of the two axial ligand binding sites of iron (see Table 7). Corresponding 
changes will be studied for MbOə2, HbOg. It is specific to the three 
sterostructural models used for Mb, Hb (Fig. 23) that the ferrous iron is 
pentacoordinated and out of the heme plane. This assumption seems 
to be reasonable from the interpretation of X-ray data of Muirhead 
et al. on deoxyhemoglobin (82). The MO calculations for Mb, Hb have 
been carried out with the axes convention as defined in Fig. 23; this is 
only as a reminder that the EFG component is probably oriented more 
parallel to the heme than perpendicular to it. For MbCO, HbCO the com- 
mon convention used is that the z-axis is parallel to the heme normal. 
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Fig. 22. Molecular structure of the heme group (a) and of MbCO or HbCO (b) for 
MO calculations. The proximal histidine group was replaced by a water molecule 


4. Isomer Shift 


The isomer shift in Méssbauer spectroscopy, 6, depends upon the electron 
density at the nucleus under study, e(0). Values of ó are proportional 
to changes in o(0) from a reference compound: 6=«-@(Q), where « 
is known as the calibration constant, having a value of x= — 0.249 + 
0.025 ao mm/sec (83, 84) for Fe5?. An estimate of ó or o (0) from approx- 
imate MO calculations is based on ideas developed by Flygare and 
Havemeister (85). Since only Fe 3d, As, 4p, and ligand 2s, 24 atomic 
orbitals (AO’s) are included in the iterative extended Hiickel (IEH) and 
CI-calculations, the Fe 1s, 2s, 3s contributions to o (0) have to be estim- 
ated by taking into account: (i) the lack of orthogonality of Fe Is, 2s, 
3s AO’s to IEH MO’s or CI MO’s &;, and (ii) the shielding effects on 
l, (0)|2 due to the specific electronic configuration 3474s" of the com- 
pound under study. The electron density at the Fe nucleus may then 
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a N 


Fig. 23. Stereostructural models for Mb or Hb: (1) a = 0.4 A, b = 2.09 A; (2) a = 0.4 
Å, b = 2.34 A; (3) a=0.8 Å, b = 2.34 Å. The insert defines the orientation of Vzz 
with respect to the heme coordinate system 


be written (83): 


3 
o(0) =r] ¥as(0)|* +2 - > |#a0)|2 , with (10) 


$s (0) = Nas(Was(0) — > < gil Pss > ,(0)) , 


t+3¢ 

$o (0) = N 2s( F 2s(0) = 2 < gil Fas > 8i(0)) , and (11) 
¿+2ç 

$13(0) = Nas(Pis(0) — Z < gil Pas > g (0). 
ilg 


In Eqs. (10) and (11) Yans refers to Fe AO’s and ¢; to IEH MO’s or 
CI MO’s, depending on the level of approximation. For the diamagnetic 
(low-spin ferrous) compounds MbCO and HbCO the use of closed-shell 
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IEH MO’s will be a reasonable approximation, whereas one must take 
open-shell CI MO’s for the description of paramagnetic (high-spin 
ferrous) compounds like Mb and Hb. The wave functions ó; and di 
in Eq. (11) are given by & =ġ3s +ó and di =des +-$ss +. Thus in 
o(O)all cross-terms Wys(0)Yms(0) which have been attributed to be 
significant (83, 86, 87) are included. The Fe 4s population factor z, the 
normalization factor Nas, and the MO’s ¢; are determined through 
IEH MO calculations for MbCO and HbCO and through CI MO calcula- 
tions for Mb and Hb. The evaluation of overlap integrals <¢;|Wns> 
is based on Hartree-Fock atomic orbitals (88). The values Y„,(0) in Eqs. 
(10) and (11) were taken from fully relativistic Dirac-Fock calculations, 
carried out by Freeman et al. (89). The actual values of F „,(0) used either 
for MbCO, HbCO or Mb, Hb were determined for the appropriate 
electronic configuration 3d94s" of iron by interpolation between the 
values given in Table 4. The direct ligand contributions Y1jg(0) at the 
iron nucleus entering Eq. (11) have been estimated by using STO’s 
and were found to be negligible. 


Table 4. Electron densities (in ag?) for iron as derived from fully relativistic mixed 
configuration Dirac-Fock calculations. (Taken from Ref. 84) 


Configuration 3474s1 = 3d? 3d64sl 346 345451 345 

|, (0)|2 6799.018 6799.054 6798.881 6798.917 6798.653 6798.688 
|Y2,(0)]2 644.299 644.259 644.231 644.169 644.225 644.153 
| ¥%35(0)|2 87.753 87.477 88.727 88.480 90.380 90.206 
4, (0)|° 3.092 5.575 8.213 


Comparing electron densities ọ(0) of similar systems, t.e. e(0) oÍ 
MbCO, HbCO with o(0) of Mb, Hb, it may be believed that systematic 
uncertainties in ọ(0) for example due to the approximate character of 
the IEH and CI calculations are cancelled in evaluating A 0(0). 

The calculated change in ọ(0) under complex-formation of Mb, Hb 
with CO is then compared with the difference in experimental isomer 
shifts, Aó =6mp, up — Ompco, moco. By applying the Eq. 4d = « - o(0) 


it is possible to compare the ratio £ with the known value of «. In 


Table 5 the electron densities gns(0) calculated as described above for 
MbCO, HbCO (model 9, see Fig. 22b and Table 7) and for Mb, Hb 
(model 1, see Fig. 23) together with the calculated orbital occupancies 
3d%4st and the experimental values of ó are summarized. Taking Ao 
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Table 5. Experimental isomer shifts ó, calculated electronic configurations 3da 4s" 
for iron, and caculated electron densities 025(0)°) at the Fe nucleus for MbCO, 
HbCO and Mb, Hb 


Compound exp. 68) Configuration?) 015(0) 025(0) 035(0) 04s(0) 
(mm/sec) 


MbCO, HbCO 0.27 +0.03 348.36 450,18 13600.9 1280.9 184.6 0.82 
Mb, Hb 0.93+0.03 346.25 450.22 13599.8 1283.2 180.46 1.16 


a) Values of 6 versus metallic iron are mean values taken from MbCO, HbCO and 
Mb, Hb of various species. 

b) Configurations 3d@ 4s7 as derived from IEH MO calculations for MbCO, HbCO 
and from CI MO calculations for Mb, Hb. 

e) Values of ọņs(0) as derived from Eq. (10) and (11) of the text. 


and Aó from Table 5 the isomer shift calibration constant is derived as 


æ = — 0.253 + 0.023 ae mm/sec, which is in reasonable agreement with 
the value «= —0.249 40.025 aj mm/sec from other investigations 
(83, 84). 


From examination of the heme geometries of Figs. 22, 23 and the 
results listed in Table 5 it is obvious that the drastic decrease of ó, 
t.e. the drastic increase of ọ(0) during the reaction of CO with Mb 
or Hb, is mainly due to the increased overlap distortion contributions 
— > <¢di|ss>¢:(0) >0 and Ns, >1 to dgs(0), since in <¢;|¥35> the 
contributions <¥coļ¥3s> for MbCO, HbCO are considerably larger 
than the corresponding contributions <Yy,0|W3s > for Mb, Hb. In other 
words: due to the stronger interaction between ligands and iron in MbCO, 
HbCO compared to Mb, Hb the Fe 3s contribution $35 (0) to e(0) is 
considerably higher in MbCO, HbCO than in Mb, Hb. 


5. Quadrupole Splitting and Electric Field Gradient Tensor 


From the calculated electronic MO structures for the various models 
of Fig. 23 one obtains quadrupole splittings 


1 211 
AEq = > e Q Vzz (1 +7) l2 (12) 


The nuclear quadrupole moment Q is taken to be +0.21 b (76, 77, 
79, 90). The main component of the EFG tensor, Vzz, is defined by 
IVzzI> |Vyy|>|V22| and results together with the asymmetry param- 
eter 
| Vex — Vyy| 

|Vzz| 





n = (12a) 
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from diagonalizing the EFG tensor V yg, which contains valence and 
lattice contributions 


Vog=Voq se Va i 2.9 =% pi ma (13) 


val 


The tensor V pg has been calculated with the equivalence 


Vy = <||EFG||>¢ (1 — R) <73> Lg. (13a) 


The reduced matrix element <||EFG||> takes the value 2/7 for d 
electrons, 6/5 for $ electrons, and 0 for s electrons. The quantity e re- 
presents the (positive!) elementary charge, and <778> is the radial 
factor resulting from taking the expectation value of (8z2 —7?)/r®. The 
<v-3>gq4 are taken from estimates (97) based on Clementi’s (88) atomic 
Hartree-Fock wave functions: 4.49 a.u. for Fe configuration 3d’, 5.09 
a.u. for 346, and 5.73 a.u. for 345. The actual value of <v-3> gg is 
determined for each model from its calculated Fe 34 orbital occupancy 
by interpolation between the foregoing values. The Sternheimer shielding 
correction (92) is taken as (1 — R)3q==0.68. Since we find for the three 
models of Mb, Hb the configuration 3¢6-25 from MO calculations we 
derive for AEQ (Mb, Hb) =(1/2) e2Q (1 —R)3a<r-3>34 4/7 the value 
4.1 mm/sec. For Fe 4% electrons the quantity (1 — R)4p <77? >4p Is 
taken to be 1/3 of the corresponding quantity for the 3d electrons (77, 
93). Finally the tensor Lpg results from an appropriate summation over 
expectation values <¥p |p| Pn >, namely 


Lya = 2 <Pr |in] n> , (13b) 
h’h 


where y, and Z, are single-electron atomic orbitals (AO) of iron, and 
the single-electron operators lpg acting an ¥p are: 


lng = Z (lala + 14) = 5 + 154. (14) 


For iron 3d, 4s, and 4 AO’s all possible matrix elements are 
summarized in Table 6. 


6. Spin-Orbit Coupling and Temperature-Dependent Electric Field 
Gradient Tensor 


In high spin ferrous compounds we are concerned with electronic states 
IT,y,S,ys> characterized by irreducible representations T' and by total 
spin S =2. All irreducible representations T of the regarded symmetry 
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group Cey are one-dimensional, and thus y can be neglected. For con- 
venience the spin states |S =2,ys> are linear combinations of the spin 
states |S =2,ms>, which transform according to irreducible representa- 
tions of Cay. This can be achieved similar to the construction of d- 
orbitals |yz>, |xz>... from |L =2,mr>. 

The five-fold spin degeneracy of each state is lifted by spin-orbit 
interaction: 


Aso, = — 41 - Š, (15) 


where the coupling constant A is commonly related to the free ion 
value, Ap =103 cm-1, by a reduction factor «?:4= a2 Ag (94). The 
eigenvectors jez> of this problem, having energies E, are certain 
linear combinations of the base vectors |I’,S, y |> 


le> = È Cry, S, y>. (16) 


Ps 


The coefficients C7 Taye and the corresponding energies Eg, result 
from diagonalizing the spin-orbit interaction matrix <I”,S,ys|As.o.|T, 
S,ys>. Since we are interested in quadrupole splittings, A Eg, we cal- 
culate for each state |ex> its relevant EFG tensor: 


Vig = <ealV palea> , (17) 


where Vyq is defined through Eqs. (13) and (14). In order to evaluate 
the temperature dependence of the EFG one must average the com- 
ponents Vpg of the individual substates |e,> according to Boltzmann 
statistics: 


<Vog>r = 2, Vya exp ( ar] 2, exp ( (18) 


Diagonalization of <V yg> m leads to Vzz(T), to n (T), to the orienta- 
tion of Vzz(T) and finally through Eq. (12) to A Ea(Z). 

The spin-orbit interaction matrix </”,S, Ys IH s.o.|T',S,ys> and the 
EFG tensor elements <J”, S, vsa LS, ys> may be simplified by 


<I",S,ys|AL S| PS, y> = à <P |L] T> <S, y6|$|S,ys> 
<I”,S, ys |V palT, S, ys> = <T |P p| T> Öyiys . 


The evaluation of matrix elements <I” |O|I'>, with Ô =V pq or Lp 
are discussed in the following section. 
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7. Spin-Orbit Matrix and Electric Field Gradient Tensor Derived from 
Many-Electron Molecular Orbital Wave Functions 


From the Hückel one-electron wave functions ¢,= J akn Wn, with 
h 


dr representing MO’s and t, AO’s, configurations (CFG) are built up 
by constructing appropriate Slater determinants |CFG> ={&k1, éz2,..., 
dim}. Here óx¿ is a molecular spin orbital (MSO). If the CFG’s describing 
high-spin ferrous states are different from each other by only a single 
MSO, which may be denoted by x, we write |CFG, z> ={dr1, $z, . 
Pkm—1, Pk}. 

Taking into account Coulomb interaction (configuration interaction 
CI) one gets new terms |CI,/> which are in general linear combinations 
of |CFG, >: 


..} 


ICI, 2> = > mg|CFG, k>. (20) 
k 


With the wave functions defined in Eq. (20) we now calculate matrix 
elements of the operator Ó, which stands here for the EFG tensor 
operator V pg or the orbital angular momentum. We get 


<CI,/’|Q|CLI> = > mr mg <CFG,R'|Q|CFG,k>. (21) 
k',k 
The matrix elements <CFG,X |O|CFG, k> in Eq. (21) are repre- 
sented by (95) 
m-1 
<CFG,&|Q|CFG,k> = <#|0|é> + I <¢uil@lbusc> (22a) 
i=1 
for È =k, and 
<CFG,& |O|CFG,4> = <ér |Oldg> (22b) 
for Å Æ £. 


Using the convention k =&, k; the quantities <ó: |Ö]ldr> are given 
in terms of AO’s by 


<br (Olp > = > arw arn <Pr |Q|Pr> . (23) 
h’,h 
Thus Eq. (21) becomes: 


<CL/’|Q|CL/> = 2 Cun wn <Pr Ô, Pr> , (24) 


where the coefficients Cn wh) are defined by 
Can awn = svi Cwn + Can (h'h) > 
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Ch’n = > Akk Ah, 
k 

Cw ih) = > Mrk Meh’ 2, MikMen . 
k k 


Now, concerning <¥|O|¥,> of Eq. (24) we distinguish three 
different matrix elements <P, On >, <Vn:|O|Pa >, and <P Ô| 
Wr > where the indices h, # stand for Fe AO’s and Z”, h” for ligand 
AO’s, respectively. With this convention Eq. (24) can be rewritten into: 


<CL/’|Q|CLI> = > Canwn <Pr 0) Pr > 
hh 
+ 2, (Conan + Canam) <Pr Q| > (25) 


+ > Cw (h’’h’’’) <P yp ||, > š 
hh” 


The “+” sign in Eq. (25) corresponds to the case where Ô represents 
a real operator (EFG), the “—” sign to the case where Ô represents 
a complex operator (orbital momentum). The matrix elements 
<t, |Ô|Pn >, representing cross-terms between iron AO’s and ligand 
AO’s, can be transformed by introducing a completeness relation of 
orthogonal terms |¥z>, which contain the iron AO’s [Fp >: 


<p: Ô| Pr” > = > <Pr ÔER > <Ph| Pr > : (26) 
h 


Since Q does not lead out of the subspace of |Pw >, the summation 
can be restricted to iron orbitals, and we get: 


<Por Ol Pr" > = > <Pr On > San , (27) 
h 


where Syn’ stands for the overlap integrals between iron AO’s t, and 
ligand AO’s Yp”. 

From Eqs. (25) and (27) matrix elements for the spin-orbit inter- 
action and for the EFG tensor can be calculated. For the evaluation of 
matrix elements <A,.o. > from Eq. (25) we use the “ —” sign and neglect 
the pure ligand contributions <Yp |Ô|Fn >. In the case of the EFG, 
Ving is derived using the ‘“-+” sign in Eq. (25). Two parts can be distin- 
guished: (i) the so-called valence contribution 


1 
Ya = > (Cwnwn +> (Carp0 + Car pü nn) San) 
h',h h” 


<Por |P pal Pr > 


(28) 
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with <p |P pln > being obtained from Eqs. (13a), (13b), and (14), 
and (ii) the so-called lattice contribution 


t 
Vor = > Cenarny <Pn |P y|, >. (29) 


h ty h’’’ 


Assuming the ligand charges qa to be localized at lattice points (Rza, 
Rya, Rza) we can approximatel L: of Eq. (29) to 


Vpg = (1—yoo) > ga (8 Rpa Ryn — Öpg Ra) | Ra, (pq =x, 9,2), (30) 
& 


where ga and (Rra, Rya, Rza) are the charge and the cartesion coordinates 
of the ath ligand, and (1 — y,,) represents the Sternheimer antishielding 
(92) of iron-core electrons due to ligand charges. 

Information is obtained from the MO calculations concerning the co- 
efficients Cart) (h’ h), Cw (hR) and C DR RD the overlap Shp’, 
and the atomic charges ga. Thus the problem of calculating 4 E Q(T), 
n(T) and the orientation of the EFG component Vzz is easily solved 
using the formalism described in the foregoing sections. 


8. Discussion of Molecular Orbital Quadrupole Splittings for MbCO, 
HbCO 


From IEH MO calculations for MbCO, HbCO total energies E, quadru- 
pole splittings 4 Eq and atomic charges gq of iron and its nearest neigh- 
bors have been derived for ten different molecular structures (96). The 
results in Table 7 show that the original molecular geometry of Zerner 
et al. (73) (no. 1 in Table 7) results in a 4Eq far from the experimental 
value (4 EQ, exp ~ 0.35 mm/sec). Adding a water molecule to the system 
lowers E but increases A Eg (no. 2). Rotating the water around the z-axis 
in steps of 20° has no significant effect on E or A Eg. Moving the iron 
into the plane of the four nitrogens reduces both E and A Eq (no. 3). By 
taking away the water molecule and reducing the Fe—C distance to 1.64 A 
(no. 4) E increases but A Eg decreases. The total energy is decreased by 
changing the heme geometry from planar (Z) to nonplanar (K) (no. 5), 
and A Eg was lowered further. Therefore it was used for the remaining 
calculations f molecular geometry K, without water, and iron in the 
plane of the four nitrogens; all three conditions have the effect of reducing 
A Eg in the direction of the experimental value. Further reduction of 
AEg and E was obtained by reducing the Fe-C distance to 1.54 A. 
We then rotated the CO molecule from 8 = 180° to 135° in steps of 
15° (nos. 6-9). The total energy increased slightly. However, our estimate 
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of A Eg indicates that a bending angle 8 ~ 135° is more favorable than 
a linear geometry. A further decrease of the Fe-C distance seems un- 
desirable, because both E and A Eg increase again (no. 10). For all geom- 
etries, the net charge calculated for the iron atom is in agreement with 
previous calculations for other low-spin compounds (73, 76). Calculated 
spin multiplicities (96) and the 4Eq data of Table 7 for the ground 
electronic state both favor a bent geometry (no. 9) for Fe-C-O in MbCO, 
HbCO, in opposition to the fact that the linear geometry has a slightly 
lower calculated total energy. The indicators for the bent geometry are 
probably more reliable than the total energy, as simple approximate 
closed shell MO studies are known to be far better at describing gross 
charge distributions than they are at predicting electron repulsion ener- 
gies. Incidentally, it should be noted that the present results suggest 
that it might be unwise to attach too much significance to the energy 
minimum found recently by Halton (97), and which we also reproduced 
(98), (at a closed shell level of approximation) for bent Fe-O-O in oxy- 
hemoglobin. 

Overall, the results for the bent geometry are remarkably close to 
those found for the linear geometry by Zerner et al. (73, 91). This near 
coincidence is due to the fact that their study used different atomic 
parameters than ours. Zerner’s parameters were chosen to fit transition 
energies of the porphyrin system itself, while ours (72, 76-80, 96) were 
based on a total of approximately 60 molecules of diverse structure, 
non-heme compounds, and not all containing iron. We therefore argue, 
that our parameter set may be the more reliable for discussing changes 
in predicted spectral properties with changes in geometry. A bent struc- 
ture is also consistent with X-ray studies on the carbon monoxide com- 
plex of erythrocruorine (99) which indicate an Fe-C-O angle of 145° + 15°. 
However, note that our calculations do not take into account any side 
chains which might stabilize the CO in a bent geometry. 

The final conclusion from a MO interpretation of Mossbauer results 
then is that in MbCO, HbCO the bonding between the heme iron and the 
proximal histidine is very weak, the iron is situated in or nearly in the 
heme plane, in agreement with recent findings (99), the Fe-C distance 
is somewhat smaller than the typical iron-carbonyl distance of 1.84 A; 
and a Fe-C-O bond angle of <180° is more likely than a linear geom- 
etry. 


9. Remarks Concerning MbO2, HbO2 


An estimate of AEg for MbO2, HbO> by Werssbluth et al. (97) on the 
basis of MO results of Zerner et al. (73) for a heme model with O> binding 
parallel to the porphyrin ring (Griffith model), led to AEg=2.1 mm/sec 
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with Vzz<0. This is in agreement with experimental findings (see 
Fig. 8b and 10); Vzz, here, is oriented parallel to the heme plane. 

Preliminary MO results for MbO2, HbO% on an approximate closed 
shell level have been presented recently (98). These results indicate that 
there are energetically close-lying, unoccupied MO terms above the high- 
est doubly occupied MO’s. Recent Heitler-London calculations on the 
bonding of oxygen to hemoglobin by Otsuka et al. (100) indicate, that 
paramagnetic electronic configurations indeed play an important role 
for the explanation of the diamagnetism of MbO2 or HbOs. For these 
reasons the use of configuration interaction (CI) seems necessary to 
derive a satisfactory picture of structure and bonding of Fe-Os in MbO3, 
HbO;. Such CI calculations for several stereostructural arrangements of 
Oz with respect to the heme are in progress. 


10. Discussion of Molecular Orbital Quadrupole Splittings and 
Susceptibilities for Mb, Hb 


For the following discussion spin-orbit interaction has been restricted 
to states which are defined by Fig. 24 and which have been found to be 
lowest in energy by CI calculations. According to Cay symmetry of the 
three models (Fig. 23) presently under study these states may be 
characterized by irreducible representations 5A 2, ®Bı, and 5B2. To study 
the effect of an energetically low-lying, low-spin term on 4E Q(T) and 
n(T) we added the state 14, (Fig. 24) to our investigation. Since the 


IA, 

—— e = 045864| 32-1?> - 0.7977 | #-y’> + iig. 
T ome +— — Ó, = 0.7267 | 327-1? > - 041781 - y> +lig. 
+— — — —- Pey = 0.9944 yz > +ligand contributions 
— — — — b. = 0.9223|xy > + 
++ — — +4 Ó, = 0.9174[xz > + 


—u eee 


doubly occupied MO states Q; i=1-61 


Fig. 24. Many-electron states which are used to calculate AEq(T), n (T), and x (T). 
The coefficients of d44n(n = 1, ..., 5) represent the situation for model (1) of Fig. 23 
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four states have different irreducible representations the coefficient 
matrix mi of Eq. (20) will remain diagonal, independent of the energy 
separations of the four states. Thus, one gets the simple relation 
ICI, > =| CFG, &> instead of Eq. (20). This allows variation of the 
mutual energy separations between 542, 5B, °Bg and 14,1 without 
changing the states |CI,/>. 

For each of the three models in Fig. 23 calculations along the lines 
described above have been carried out (707) to obtain AF Q(T), n(T), 
sign Vzz, and the orientation of Vaz relative to the heme coordinate 
system. Some from the variety of resulting AEg(T) curves for various 
energies E (5Ag), E (?Bı), E (5B2) and E (141) are shown in Figs. 25-27 
for spin-orbit coupling constant A=90 cm-t, and in Figs. 28 and 29 
for A=75 cm-t. The “+ or “—” sign at each curve indicate the 
calculated sign of the EFG component Vzz. The energies, asymmetry 
parameters, and orientations of Vzz corresponding to each AE Q(T) 
curve of Fig. 28a and 29a are tabulated in Table 8. Experimental 
AEg{T) data for frozen Mb solution (53) are indicated by closed circles, 
for a Mb single crystal (42) by “x” and error bars, and for frozen Hb 
(rat) solution (52) by open circles and error bars; the sign of Vzz was 
experimentally determined to be positive from Mb single crystal Möss- 
bauer measurements (42) and from magnetic hyperfine split Méssbauer 
spectra of Hb (57) (Figs. 9c and d). 

Comparing theoretical AE (T) curves with experimental 4 £ Q(T) 
data from Mb and Hb reasonable coincidence of computational and 
experimental results is found only for curve b in Fig. 28a. From inspection 
of Table 8 it is obvious that this AEgQ(T) curve corresponds to an EFG 
component Vzz being oriented along the heme axis y. This result agrees 
well with the recent findings from Méssbauer investigations of Mb single 
crystals at 77 °K (42). There an angle of x =40°-+8° was found between 
Vzz and the crystallographic b axis of the Mb single crystal, indicating 
that Vzz possibly is oriented along one of the two heme axes y or z. 

The asymmetry parameters 7(7) (from Table 8) which are related 
with the A Eg(T) and fit curve b of Fig. 28a also agree with the Mb 
single crystal Mössbauer results. In Ref. (42) the angle «=40°+8° 
was derived assuming 7 to be zero. A most recent computational analysis 
of our former Mb single crystal Mössbauer data by Maeda et al. (62), 
however, indicates that «=40°-+8° is consistent also with an n-para- 
meter in the range of 0< 7<0.3. 

The energy separation E (14,)-E(2B2) which was varied in some 
cases (Figs. 25-27) between 50 and 10,000 cm-1 is sensitive only to the 
slope of the theoretical AEg(T) curves for temperatures T > 50 °K, 
and, further, AEg(T) curves up to T=200°K are identical for 
E(1Aı) — E (2B) =300 cm! and £F(1A1) —E(5Bs)=10,000 cm}. 
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Table 8. Energies of spin quintet states (in cm~1), temperature dependent quadru- 
pole splittings (in mm/sec), and asymmetry parameters. The angles a and ß (in 
degrees) defining the orientation of Vzz with respect to the heme coordinate system 
of Fig. 23 are for all curves 28a and 29a and for temperatures 4.2 < T < 200 °K 
given by « = 90°, B =90° 


Figure Curve Energies of parameters for temperatures T 


5Ba 5B, A5 42 50 75 100 150 200 


a 0 150 150 AEq 2.157 2.118 2.139 2.101 1.934 1.748 

n 0.461 0.436 0.420 0.419 0.447 0.494 

b 0 150 175 AEg 2.228 2.211 2.235 2.196 2.019 1.820 

28a N 0.355 0.341 0.327 0.325 0.341 0.377 
c 0 150 200 AEg 2.299 2.303 2.330 2.290 2.103 1.891 

N 0.249 0.246 0.236 0.230 0.235 0.261 


a 0 120 160 AEq 2.242 2.181 2.204 2.175 2.025 1.850 

29a N 0.369 0.326 0.309 0.307 0.320 0.336 
b 0 120 180 AEg 2.313 2.269 2.294 2.265 2.112 1.931 

n 0.421 0.371 0.352 0.352 0.372 0.394 


Though we find AE. (T) curves corresponding to E(1Aı) — E(?Ba) = 
100 cm-1, which fit the experimental 4 Es (T) data, we present arguments 
why an energy separation of E (141) — E (5B2) 3300 cm-1 seems to be 
more realistic. In addition to the static molecular structure of Fig. 23 
one has to consider lattice vibrations which have the effect that the 
potential V of the iron contains higher terms than the static ligand field 
potential V( (702): 


V = VO + VM e + VO er 4 ++, (31) 


The quantity e in Eq. (31) represents an isotropic average of the 
strain tensor, and VD, VO) ... describe orbital-lattice interactions. 
The transition probability of a direct process between states |a > and 
|b > is given by 


Way = y |<a, N — 1|VO) e | 5, N > |2 o (E) (32) 
= |<a|V0O |b>[2|<N —1|#|N|2o (E), 


where N = occupation number of phonons, and o (E) =phonondensity 
of final states. Wap is zero if |a > stands for a spin quintet and |b > for a spin 
singlet, because VD is spin independent and spin-orbit coupling mixes 
high-spin and low-spin states via intermediate spin-triplets only. Thus 
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0 50 100 150 200 TIK) 





0 50 100 150 200 T(K) 


Fig. 25. Temperature dependent quadrupole splittings for model (3). 4 = 90 cm~1. 
Experimental AEg(T) data for frozen Mb solution [Ref. (53)] are indicated by closed 
circles, for a Mb single crystal [Ref. (42)] by “x” and error bars, and for frozen Hb 
(rat) solution [Ref. (52)] by open circles and error bars. Energies 5Bg, 5B, 542, and 
14; in cm-1l corresponding to AEg(T) curves are: 


(25a) a 0 150 100 50 (25b) a 0 50 300 50 
b 0 150 100 10000 b 0 50 300 10000 
c 0 200 100 50 c 0 100 300 50 
d 0 200 100 10000 d 0 100 300 10000 
e 0 300 100 50 
f 0 300 100 10000 
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AEe 


a) 





0 50 100 150 200 T(K) 


1.9 





0 50 100 150 200 T(K) 


Fig. 26. Temperature dependent quadrupole splittings for model (2). ¿ = 90 cm-1. 
Experimental A Eg(T) data are defined in Fig. 25. Energies for SBə, 5B1, 54 and 
14, in cm! corresponding to AEg curves are: 


(26a) 


+ oO moO me» 


0 
0 
0 
0 
0 
0 


300 400 50 (26b) a Q 400 200 50 
300 400 10000 b 0 400 200 10000 
300 600 50 c 0 400 400 50 
300 600 10000 d 0 400 400 10000 
300 1000 50 


300 1000 10000 


<a|V [> =0. In the real case the transition probability Wa, might 
be slightly different from zero since the state |a > (|b >) might be a linear 
combination of a spin quintet (spin singlet) state and a small contribution 
of its counterpart due to spin-spin interactions. 
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a) 





0 50 100 150 200 TIK) 


b) 





0 50 100 150 200 T(K) 


Fig. 27. Temperature dependent quadrupole splittings for model (1). A=90 cm-1. 
Experimental AEg(T) data are defined in Fig. 25. Energies for ®Ba, 5B,, Aa, 
and 14, in cm-1 corresponding to A Eq curves are: 


(27a) a 0 150 150 50 (27b) a 0 100 300 50 
b 0 150 150 10000 b 0 100 300 10000 
c 0 200 150 50 c 0 150 300 50 
d 0 200 150 10000 d 0 150 300 10000 
e 0 300 150 50 
f 0 300 150 10000 


W ap being very small, however, means that the spin-lattice relaxation 
between |a> and |b > is very slow. With the spin-lattice relaxation 
time becoming larger than the precession time of the Fe? nucleus 
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two separate Mössbauer spectra would appear, one corresponding to 
the spin quintet |a> and the other to the spin singlet |b >, and the 
thermal average for evaluating AEg(T) should be only over the high- 
spin states. There are several examples of such high-spin low-spin tran- 
sitions in the literature (703) for compounds containing ferrous iron. 
In believing in slow spin-lattice relaxation between high- and low-spin 
states in deoxygenated myoglobin and hemoglobin one might try to 
explain Mössbauer spectra of deoxygenated hemoglobins, containing 
small-intensity lines (see Fig. 9a) by thermal mixtures of high- and low- 
spin quadrupole splittings. The AEQ value associated with the 14, 
state of about 1.3 mm/sec would be just in agreement with the “impurity” 
line positions of Fig. 9a. However, the intensity ratio, R, of low-spin 
to high-spin lines should be temperature dependent, in contrast to the 
experimental results for example of deoxygenated hemoglobin single 
chains (704), which sometimes under certain preparational conditions 
also show these “impurity” lines. Assuming the 1A state to lie in energy 
above the spin quintet ground state by 50 cm-1 (300 cm-1) the intensity 
ratio R would be unmeasurably small at T =4.2°K; at an elevated 
temperature of T =300°K the ratio R is estimated from energies Ex 
to be about 8% (2%). For these reasons we believe the energy separation 
EF (1A1) — E (5B2) to be at least 300 cm-t, thus the spin singlet 14, 
does not influence our JE Q(T) curves at all. 

To assess the reliability of the MO results further the energy terms 
Ex corresponding to curve b in Fig. 28a are taken to compute temper- 
ature-dependent magnetic susceptibilities z (T). The y-tensor ypq with 
P, g = x, y, z is defined by 


OMp 


Xpa = ƏHq |Ha ~ 0° (33) 


Mp represents the magnetization per mole. According to Eq. (33) 
we may write 


My = 2, xpa Ha ; (33a) 
q 


and the tensor elements ypq we derive from 


Ara = up A > <eallp +2 Slew > Lea |La +2 Ŝg lea > 
exp (— Eq’ | kT) — exp (— EalkT) 


(34) 
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Fig. 28. Temperature dependent quadrupole splittings for model (1). A=75 cm-1, 
Experimental A E Q(T) data are defined in Fig. 25. Energies for 5Ba, 5By, 5Ag and 
14, in cm“! corresponding to A Eg curves are: 


(28a) a 0 150 150 10000 (28b) a 0 200 100 10000 
b 0 150 175 10000 b 0 200 150 10000 
c 0 150 200 10000 


The spin states |ex > have already been defined by Eq. (16), and the 
Es are their corresponding energies. Symbol A stands for the Avogadro 
number. For cases of degeneracy we simply take 


lim exp (-Ea/kT) — exp (—Ea|RT) exp (—Ea/kT) 
Ea’”Ea Eq — Ea kT I 


The final powder susceptibility z (T) results from x(T)= (1/3) > xp», 
p 


where $ó =x, y, z. The solid y(T) curve in Fig. 30 corresponds to the 
A Eg (T) fit curve b in Fig. 28a. Closed circles in Fig. 30 are experimental 
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Fig. 29. Temperature dependent quadrupole splittings for model (3). A =75 cm-1. 

Experimental A Eq(T) data are defined in Fig. 25. Energies for 5B>s, 5B1, 541 and 

14, in cm-1 corresponding to A Eg curves are: 

(29a) a 0 120 160 10000 (29b) a O 180 100 10000 
b 0 120 180 10000 b 0 180 120 10000 


points taken from the susceptibility investigation of Mb and Hb by 
Nakano et al. (69). The y(T) curves corresponding to JE Q(T) curves a 
and b in Fig. 29a are nearly identical to the solid line of Fig. 30; the 
broken line, however, represents the situation that spin triplet states lie 
close to the spin quintet ground state with energy separations 
E (3E n) — E (B2) =200 cm-t and E(?EE) — E (5B) =300 cm-1. Thus, 
by comparing solid line, broken line, and experimental points in Fig. 30, 
we conclude that the energy sequence 5B>, 5B1, 5A ı gives a more realistic 
picture for deoxygenated heme compounds than the formerly derived 
ligand field result (52) with an energy sequence of 542, 3E (105). 
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0 50 100 150 200 T(K) 


Fig. 30. Temperatures dependent suceptibilities. Solid line corresponds to energies 
as specified for 4 Eg(T) curve b of Fig. 28a. Broken line corresponds to energetically 
low-lying spin triplet states (see text). Experimental points are taken from Ref. (69). 
x is given in emu per mole 


Summarizing the results, the present MO interpretation of experi- 
mental Mössbauer and susceptibility data of deoxygenated heme pro- 
teins is approximate, in the sense that in the three models under study 
(Fig. 23) we neglected the protein part of the molecule, several peripheral 
side chains and the probable non-planarity of the heme group, that we 
represented the histidine group by a water molecule, and further, that 
we assumed all interactomic distances to remain constant within the 
whole temperature range of 4.2°K<7T<200°K; moreover, the MO 
procedure itself is approximate as described in Refs. (72) and (76). These 
simplifications, however, are believed to be secondary in the present 
attempt to derive a gross description of the structural situation of the 
heme iron in deoxygenated heme compounds. From the work described 
it is therefore concluded that the MO interpretation of experimental 
Méssbauer and susceptibility data consistently agrees with the assump- 
tion that the heme iron in Mb and Hb is pentacoordinated and signif- 
icantly out of “plane” (by 0.4— 0.8 A, Fig. 23). The interpretation of 
recent X-ray data of Mb (706) and of Hb (82) partly support these find- 
ings, namely that the heme iron in Mb, Hb is penta-coordinated. 


VIII. Experimental Solution of the Phase Problem for Protein 
Structure Determination by Mössbauer Scattering 


The principle of phase determination by utilizing the interference be- 
tween electronic and nuclear resonance scattering and its application 
to proteins has been described in detail by Parak et al. (107). Here, we 
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briefly review the expectations and difficulties concerning the appli- 
cation of this new method to heme proteins. 

For the determination of the electron density map of a molecule the 
amplitudes and phases of the waves scattered by a single crystal are 
required for a number of Bragg reflections. Common X-ray techniques 
yield the product (Ayelpy): (Ane pz) —Ap, where Ay is the 
amplitude and øg the phase of the sum of all waves scattered by the 
electrons of the molecule into the Bragg reflection H. Thus, the scatter- 
ing amplitude Az is obtained, but the phase information is lost. The 
solution of this phase problem for protein structure determination is 
based on Perutz’ and Kendrew’s isomorphous replacement method 
(708-1717). In this procedure Bragg reflections have to be measured at 
least three times, first on a crystal of native molecules, and then on two 
crystals, in which reference scatterers (for example Hg atoms) have been 
substituted at well-defined positions. From the difference of the measured 
intensities one can calculate the relative phases without ambiguity. 

A supplementary method for the isomorphous replacement method is 
to provide an additional scatterer in the crystal. Applying the Mossbauer 
scattering technique (see II.5 and Fig. 2b} to single crystals of heme- 
proteins has the advantages that the Fe5’ scatterer is originally built into 
the protein, and that the nuclear scattering amplitude theoretically 
amounts to 6.215 times the scattering amplitude for the 80 electrons of 
one mercury atom. Therefore, Fe5? obviously is well suited as a reference 
scatterer in iron-containing proteins. Further, the wave-length of the 
Fe5? 14.4 keV y-radiation of 4 = 0.86 A is comparable to atomic spacings 
in the crystal, and the radiation scattered by the Fe®? nucleus can inter- 
fere with the radiation scattered by the electron shell of an atom (772, 
113). Additionally, it is important to note, that for a unique phase de- 
termination, the Fe5? nucleus is sufficient as a single reference scatterer 
compared to the twofold heavy atom substitution in the isomorphous 
replacement method. This can be understood in the following way. The 
phase of resonantly scattered y-radiation is shifted by 5 with respect to 
the incident y-radiation from the source S in Fig. 2b. Moving the source 
with a Doppler velocity v> vy relative to the heme protein single crystal 
destroys the nuclear resonance scattering which occurs at resonance ve- 
locity vy by Doppler effect. The nuclear scattering amplitude is then 
equal to zero, and the phase shift between incoming and the non-reso- 
nantly scattered wave also becomes zero. Decreasing the Doppler velocity 
v and approaching vr increases both the phase shift and the nuclear 
scattering amplitude. By arbitrary choice of two different Doppler velo- 
cities close to vr, two different nuclear scattering amplitudes with known 
absolute value and phase can be selected as an alternative for the two 
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heavy atom scattering amplitudes of the isomorphous replacement 
method (707, 774). The construction of the phase angle yx with the use 
of the two known nuclear scattering amplitudes. mye and Are and the 
primarily known scattering amplitude Fo(H), not including nuclear 
resonance scattering, then is analogous to the procedure for the X-ray 
structure analysis (Fig. 31). 


= £, 
IE, Hl GE 






IFo (H)! 


` 
x IFo (H)! 


— ` 


tH 1 


Fig. 31. Principle of the phase determination [Taken from Ref. (774)]. npe(#re) 
is the nuclear resonance scattering amplitude, Fo(H) is the scattering amplitude 
of the unit cell, not including nuclear resonance scattering, Foa(H) (Foa(H)) 
represents the scattering amplitude of the unit cell, including the nuclear resonance 
scattering. nye and FoA(H) correspond to Doppler velocity vr, and fife, Foa(H) 


to # (ë is very clos to vr). po(H) =Go(H) is the correct phase, which should be 
determined. 


For these reasons the Méssbauer scattering technique in principle is 
a suitable alternative to the isomorphous replacement method. Several 
experimental difficulties, however, as for example the preparation of 
strong Co5? sources (up to 500 mCi/mm2), the enrichment of heme- 
protein single crystals with Fe5? (49), and the experimentation at low 
temperatures (77°K) (775), have prevented the nuclear resonance 
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scattering method for protein structure determination from becoming 
successful so far; in particular the struggle continues for the determi- 
nation of the first phase value yy. 
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